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Computer  Programs  for  Scoring  Test  Data  with 
Item  Characteristic  Curve  Models 


Although  latent  trait  test  theory,  or  item  characteristic  curve  (ICC) 
theory,  has  been  developing  since  Lawley's  (1943)  paper  more  than  30  years  ago, 
applications  of  the  theory  have  appeared  only  recently.  However,  there  are 
indications  that  latent  trait  test  theory  is  beginning  to  reach  the  practitioner 
who  is  concerned  with  test  development  and  usage  in  applied  settings.  This  is 
evidenced , not  only  by  the  increasing  number  of  journal  articles  concerned  with 
latent  trait  test  theory  (e.g. , the  summer  1977  special  issue  of  the  Journal  of 
Educational  Measurement  on  applications  of  latent  trait  models)  and  in 
presentations  and  training  sessions  at  professional  meetings,  but  also  by  its 
application  in  adaptive  (Weiss,  1976)  or  tailored  (Lord,  1970)  testing. 

A potential  disadvantage  of  latent  trait  test  theory  is  that  its  use 
often  involves  complex  computational  procedures.  To  apply  ICC  models  to  the 
development  of  tests  and  their  scoring,  the  psychometrician  must  be  able  to 
estimate  the  ICC  parameters  of  the  items  in  the  test, and  then  use  them  in 
conjunction  with  the  response  data  of  a new  group  of  testees  in  order  to  estimate 
their  trait  scores  (e.g.,  ability  or  achievement  levels).  A number  of  compu- 
ter programs  are  available  for  estimating  ICC  item  parameters  (these  are 
summarized  in  Appendix  Table  A).  However,  there  appeared  to  be  no  general 
programs  available  for  scoring  test  data  with  ICC  models  when  item  parameter 
estimates  were  available  from  previous  data  sets.  This  report  describes 
several  programs  designed  to  meet  this  need. 

An  Introduction  to  Test  Scoring 

The  problem  of  test  scoring  can  be  conceptualized  as  the  process  of 
summarizing  a testee's  answers  to  a set  of  test  questions  into  a single  number 
in  such  a way  that  the  score  will  be  indicative  of  the  testee's  position  on  the 
trait  being  measured  by  the  test.  The  most  common  test  scoring  strategy  is  to 
add  the  number  of  correct  answers  and  to  transform  the  score  into  some  type  of 
standard  score  or  percentile  to  add  interpretability.  Historically,  the  number- 
correct  score  has  been  used  because  it  is  easy  to  calculate,  and  in  pre-computer 
days  this  was  an  essential  requirement  of  a test  scoring  procedure.  As  a 
general  procedure  for  scoring  tests  of  ability  and  achievement,  however,  the 
number-correct  score  has  several  deficiencies. 

Inadequacies  of  the  Number-Correct  Score 


One  major  problem  with  the  number-correct  score  is  that  it  is  possible 
for  the  same  number-correct  score  to  be  obtained  in  several  different  ways; 
that  is,  several  response  patterns  can  result  in  the  same  number-correct  score. 
If  the  items  in  a test  are  all  of  equal  difficulty  and  discrimination,  and 
therefore  are  essentially  replicates  of  each  other,  this  will  have  little 
effect  on  the  number-correct  score,  since  different  response  patterns  among 
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replicate  items  are  of  little  consequence.  But  it  is  a very  rare  test — and 
one  which  would  have  little  general  measurement  utility — which  would  have  items 
that  are  all  replicates  of  each  other  with  regard  to  difficulty  and  discrimi- 
nation. 

When  test  items  differ  with  respect  to  difficulty  or  discrimination,  they 
are  no  longer  replicates.  Under  these  circumstances,  different  patterns  of 
response  to  the  same  set  of  items  convey  different  information  with  regard  to 
a testee's  trait  level.  The  testee  who  correctly  answers  only  five  very 
difficult  items  in  a test  is  likely  of  higher  ability  than  the  testee  who 
correctly  answers  only  five  very  easy  items  in  the  same  test.  Although  the 
total  number-correct  score  is  the  same  for  these  two  testees,  their  trait 
level  estimates  derived  from  latent  trait  or  ICC  theory  will  differ.  An 
additional  unattractive  feature  of  the  number-correct  score  is  the  fact  that 
the  number  of  possible  scores  is  determined  by  the  number  of  items  in  the  test. 
Thus,  if  a test  consists  of  only  10  items,  only  10  unique  scores  are  possible. 
Although  this  may  be  sufficient  in  some  applications,  in  others  it  might  be 
desirable  to  obtain  a finer  gradation  of  scores. 

The  inadequacy  of  the  number-correct  score  as  a general  test-scoring 
procedure  is  most  obvious  when  considering  how  to  score  responses  of  testees 
who  have  been  administered  different  sets  of  items,  as  in  adaptive  or  tailored 
testing.  In  these  kinds  of  tests,  number-correct  scores  are  completely 
inappropriate,  since  different  testees  will  receive  items  of  different  diffi- 
culties and  discriminations  as  well  as  different  numbers  of  items  in  an  adap- 
tive test.  In  addition,  the  proportion  of  correct  responses  obtained  by  all 
testees  will  be  approximately  the  same  in  a well-designed  adaptive  test  (e.g., 
Weiss,  1975). 

ICC-Based  Scoring 

The  scoring  programs  described  in  this  report  use  considerably  more 
refined  approaches  than  a mere  adding  of  correct  answers  and  are  usable  for 
scoring  both  conventional  and  adaptive  test  data.  This  refinement  is  possible 
because  ICC  theory  makes  very  explicit  specifications  about  the  relationship 
between  performance  on  a test  item  and  the  testee's  position  on  the  trait,  0. 
This  relationship  is  referred  to  as  the  item  characteristic  curve  (ICC;  Lord  & 
Novick,  1968)  when  the  items  are  scored  into  two  categories  (correct  or 
incorrect)  or,  when  there  are  more  than  two  score  categories,  as  the  operating 
characteristic  function  (Samejima,  1969). 

In  the  context  of  latent  trait  test  theory,  scoring  may  be  conceptualized 
as  finding  the  value  of  0 (i.e.,  the  trait  being  measured)  most  "compatible," 
in  some  sense,  with  a given  pattern  of  responses  to  the  test  items,  given  the 
ICC  item  parameters  for  each  item  answered.  For  maximum  likelihood  scoring, 
the  score  associated  with  a given  response  vector  is  that  value  of  0 for  which 
the  likelihood  of  the  response  vector  is  maximum.  For  Bayesian  scoring,  the 
score  is  usually  either  the  value  of  0 that  minimizes  the  mean  squared  differ- 
ence between  the  estimated  0 and  the  "true"  8,  or  the  value  of  0 that  is  most 
probable  given  the  observed  responses. 

Maximum  likelihood  scoring.  The  details  of  the  maximum  likelihood  scoring 
procedure  are  presented  below.  However,  a conceptual  explanation  based  on  two 
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dichotomously  scored  test  items  will  serve  to  explicate  its  rationale. 

Figure  1 shows  response  probability  curves  for  two  test  items — Item  1, 
which  was  answered  correctly  (resulting  in  an  ICC  plot  of  the  probability  of 
a correct  response),  and  Item  2,  which  was  answered  incorrectly  (resulting  in 
a descending  plot  of  the  probability  of  an  incorrect  response,  or  1 minus  the 
ICC).  The  ICC  curves  for  the  two  items  are  described  by  three  parameters: 

(1)  difficulty,  b , which  is  the  location  of  the  ICC  on  the  trait  (9)  contin- 
uum at  the  point  of  maximum  slope  of  the  ICC  (£>=-. 5 for  Item  1 and  .75  for 
Item  2);  (2)  their  discrimination,  a,  which  is  proportional  to  the  slope  of 
the  ICC  at b (a=. 8 for  Item  1 and  1.4  for  Item  2);  and  (3)  "guessing,"  a , the 
lower  asymptote  of  the  probability  of  a correct  response  at  0=-“>  (c=.16  for 
Item  1 and  l-.78=.22  for  Item  2). 


Figure  1 

Response  Probability  Plots  for  a Correctly  Answered  Item  (Item  1) 
and  an  Incorrectly  Answered  Item  (Item  2) 


Item  1 [P , (B) ] 


Item  2 [Q2 (B)=1-P2(0) ] 


Trait  Level  (9) 


The  first  step  in  maximum  likelihood  scoring  consists  of  determining  the 
likelihood  of  the  response  pattern  (correct  response  to  Item  1 and  incorrect 
response  to  Item  2).  Assuming  local  independence,  which  means  that  responses 
to  the  test  items  have  nothing  in  common  except  their  relationship  to  the 
underlying  trait,  0,  the  likelihood  of  a response  pattern  at  any  value  of  0 
can  be  determined  by  multiplying  the  separate  probabilities  of  the  responses 
in  the  response  pattern  for  that  value  of  9.  The  value  of  0 for  which  the 
likelihood  is  maximum  is  the  maximum  likelihood  estimate  of  0. 

Conceptually,  this  can  be  illustrated  with  the  ICCs  in  Figure  1 by  using 
discrete  values  of  0,  such  as  those  shown  in  Table  1.  For  example,  at  9--1.0, 
the  probability  of  a correct  response  to  Item  1 (scored  as  1)  is  .442  and  the 
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probability  of  an  incorrect  response  to  Item  2 (scored  as  0)  is  .768;  multiply- 
ing these  values  gives  the  likelihood  of  the  [1,0]  response  pattern  as  .340. 

At  9=+1.0,  the  probability  of  a correct  response  [1]  to  Item  1 is  .903  and  the 
probability  of  an  incorrect  response  to  Item  2 is  .277;  the  likelihood  of  the 
[1,0]  response  pattern  is  therefore  .250.  Similarly,  at  9=0.0,  the  probability 
of  a correct  response  to  Item  1 is  .718  and  the  probability  of  an  incorrect 
response  to  Item  2 is  .668,  resulting  in  a likelihood  for  the  [1,0]  response 
pattern  of  .479.  This  process  of  computing  likelihoods  for  the  [1,0]  response 
pattern  can  be  repeated  for  a large  number  of  values  along  the  9 continuum. 


Table  1 

Probability  of  a Correct  Response  to 
Item  1 [P^ (0) ] and  Probability  of  an 

Incorrect  Response  to  Item  2 [^C®)]  f°r 

Selected  Values  of  0 (Item  1:  a=.8, 

£>=-. 5,  c=.16;  Item  2:  a=1.4,  £>=.75, 

<?=.22),  and  Values  of  the  Likelihood 
Function  [£(6)] 


0 

Pl(6) 

e2(0) 

L(6) 

-3.0 

.187 

.780 

.146 

-2.5 

.212 

.780 

.165 

-2.0 

.257 

.779 

.200 

-1.5 

.332 

.776 

.257 

-1.0 

.442 

.768 

.340 

.5 

.580 

.742 

.430 

0.0 

.718 

.668 

.479 

.5 

.828 

.503 

.416 

1.0 

.903 

.277 

.250 

1.5 

.948 

.112 

.106 

2.0 

.973 

.038 

.037 

2.5 

.986 

.012 

.012 

3.0 

.993 

.004 

.004 

The  result  of  computing  likelihoods  for  all  possible  values  of  9 based 
on  the  response  pattern  and  the  relevant  ICCs  can  be  a plot  of  the  likelihood 
values  as  a function  of  9.  This  plot,  shown  as  a solid  curve  in  Figure  2,  is 
called  a likelihood  function.  As  can  be  seen,  the  maximum  of  the  likelihood 
function  in  Figure  2 occurs  at  about  0=0.0  (actually  .01).  Thus,  9=.01  can 
be  considered  the  maximum  likelihood  estimate  of  0 associated  with  the  [1,0] 
response  pattern,  given  the  parameters  of  the  ICCs  for  the  items  generating 
that  response  pattern.  The  maximum  likelihood  0 estimate  is  thus  the  value 
of  0 which  maximizes  likelihood  of  the  given  response  pattern  for  items  with 
the  specified  ICCs. 

The  generalization  of  the  scoring  method  for  more  than  two  items  is 
straightforward.  For  each  value  of  0,  the  likelihood  would  be  determined  by 
multiplying  the  response  probabilities  for  the  appropriate  ICCs  (based  on 
the  specified  response  pattern)  across  all  items  that  have  been  answered. 
Thus,  for  n items,  n probabilities  would  be  multiplied  at  each  value  of  0 to 
obtain  the  likelihoods.  The  resultant  likelihood  values  for  all  values  of  0 


-5- 


could  be  plotted;  and  the  maximum  of  the  likelihood  function  would  be  used  to 
identify  the  value  of  0 that  gives  the  observed  response  pattern  the  greatest 
probability  of  occurrence. 

Figure  2 

Likelihood  Function  and  Bayesian  Posterior  Density 
Function  for  the  [1,0]  Response  Pattern 


Maximum  likelihood  scores  are  intuitively  appealing;  at  the  same  time, 
they  have  a number  of  optimal  statistical  characteristics,  at  least  asymptotic- 
ally (i.e.,  when  large  numbers  of  items  are  administered).  Of  special 
relevance  is  the  fact  that  as  the  number  of  items  in  the  response  pattern 
increases,  it  can  be  shown  (Kendall  & Stuart,  1961)  that  maximum  likelihood 
estimates  have  minimum  variance;  and  the  reciprocal  of  that  variance  is  known 
as  the  information  function  of  0.  As  a consequence,  different  scores  (i.e., 

0 estimates)  car,  have  different  degrees  of  accuracy  as  estimators  of  0 (Birnbaum, 
1968;  Samejima,  1969). 

Bayesian  scoring.  Although  the  numerical  details  of  maximum  likelihood 
and  Bayesian  scoring  are  substantially  different,  the  two  methods  are 
conceptually  very  similar.  Bayesian  scores  are  based  on  the  likelihood 
function  modified  by  the  prior  probability  density  function  of  6.  The  prior 
probability  density  function  describes  the  assumed  distribution  of  0 in  the 
population  of  individuals  to  be  tested. 

To  illustrate,  call  L(0)  the  likelihood  of  the  response  pattern  for  a 
given  0 value.  Now  call  f(Q)  the  prior  probability  density  associated  with 
that  value  of  0.  The  modified  likelihood,  which  may  be  called  n(0)  is  then 

p(0)”/*(9)£(0)  / yi.f(0)M0)]de.  [1] 
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Equation  1 is  called  the  posterior  probability  density  function.  Just  as  the 
maximum  likelihood  score  is  the  value  of  0 for  which  £(0)  is  maximum,  one  kind 
of  Bayesian  score  is  the  value  of  0 for  which  p(0)  is  maximum.  Such  scores 
are  called  Bayes  modal  estimates  by  Samejima  (1969)  because  they  are  based  on 
the  mode  of  the  posterior  density  function.  A different  type  of  Bayesian 
estimate  is  based  on  the  mean  of  the  posterior  density  function.  Owen's 
(1975)  Bayesian  scoring  procedure,  which  will  be  described  in  detail  below, 
is  an  example  of  this  approach.  In  his  procedure,  the  prior  probability 
densities  are  provided  by  a normal  density  function.  Other  Bayesian  scoring 
procedures  are  also  available  (Sympson,  1977). 

These  concepts  can  be  illustrated  using  the  likelihoods  associated  with 
the  [1 ,0]  response  pattern  discussed  earlier.  Table  2 shows  for  several 
values  of  0 the  probability  of  a correct  response  to  the  first  item  [Pj(0)]; 
the  probability  of  an  incorrect  response  to  the  second  item  the 

likelihood  of  the  response  pattern  [£(0)]  (these  first  three  columns  correspond 
to  the  data  in  Table  1);  the  prior  probability  densities  [/(©)],  which  in 
this  case  are  ordinates  of  a normal  distribution  with  mean  of  zero  and  standard 
deviation  of  1;  and  the  posterior  density  function  [p (0) ] , computed  using 
Equation  1.  For  these  data 

yi/(0)i(e)]d9=.348.  [2] 

The  resulting  posterior  density  function,  [p(0)J,  is  shown  as  the  dashed  curve 
in  Figure  2. 


Table  2 

Response  Probabilities  [P^(0),  ( ©) J » Likelihoods  [L(0)], 

Weights  [U(0)],  and  Posterior  Density  Function  [p(0) ] for  a 
Two-Item  Response  Pattern 
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The  mode  of  the  posterior  density  function  in  Figure  2 is  located  near 
0=0,  so  the  Bayesian  modal  estimate  and  the  maximum  likelihood  estimate  are 
about  the  same  for  this  data.  The  Bayesian  8 estimate  based  on  the  mean  of 
the  p( 8)  distribution  is  -.12^  This  0 estimate  does  not  coincide  with  the 
maximum  likelihood  estimate  (0=.O1);  as  will  be  further  shown  below,  estimates 
of  0 obtained  from  different  ICC  scoring  methods  do  not  generally  agree. 

Differences  Among  Scoring  Methods 

The  programs  described  in  this  report  are  capable  of  scoring  test  data 
using  most  of  the  ICC  response  models  available.  The  selection  among  models 
should  not  be  arbitrary,  especially  when  individual  decisions  are  to  be  made  on 
the  basis  of  test  scores.  Dichotomous  data  can  be  scored  by  means  of  the  one-, 
two-,  and  three-parameter  ICC  models,  using  either  a normal  or  logistic  ogive 
ICC.  Thus,  given  the  decision  with  regard  to  the  number  of  parameters  that 
describe  the  ICC,  there  still  remains  the  problem  of  choosing  between  the 
normal  or  logistic  ogive  response  models  for  scoring  purposes.  Unfortunately , 
there  are  as  yet  no  firm  guidelines  for  choosing  between  these  two  response 
models.  Same j ima  (1969)  has  shown  that  the  normal  and  logistic  ogive  models 
differ  with  respect  to  their  scoring  "philosophies,"  but  the  practical  impli- 
cations of  these  differences  remain  to  be  investigated. 

To  illustrate  the  differences  among  the  models  and  different  ICC  scoring 
procedures,  all  response  patterns  for  a five-item  test  were  scored  by  maximum 
likelihood,  assuming  both  normal  and  logistic  ogive  ICCs,  and  by  Owen's  (1975) 
Bayesian  scoring  method.  Table  3 gives  the  item  parameters  assumed  for  the 
hypothetical  five-item  test.  For  all  items,  the  c (guessing)  parameter  was 
set  at  0.0,  indicating  that  a two-parameter  ICC  model  was  used.  Items  varied 
in  difficulty  (i>)  from  -2  to  +2  and  had  discriminations  of  1.00  or  1.50. 

Table  3 


Item  Paremeters  for  Five-Item  Test 


1 1 era 

a 

b 

C 

1 

A 

1.00 

-2.00 

.00 

2 

1.50 

-1.00 

.00 

3 

1.00 

0.00 

.00 

4 

1.50 

1.00 

.00 

5 

1.00 

2.00 

.00 

In  a five-item  test  in  which  each  item  is  scored  dichotomously , there  are 
25=32  different  response  patterns.  These  response  patterns  are  shown  in 
Table  4 along  with  the  scores  associated  with  them.  It  is  obvious  from  the 
data  in  Table  4 that  for  a given  response  pattern , the  scores  (all  of  which  are 
on  the  same  metric)  differed  somewhat.  This  indicates  that  the  scoring 
procedures  are  not  interchangeable. 

For  example,  consider  the  five  response  patterns  which  have  20%  correct, 
namely  Patterns  2,  3,  5,  9,  and  17.  Not  only  do  the  0 estimates  (scores)  for 
a given  response  pattern  differ  among  the  three  scoring  procedures,  but  there 
are  some  differences  in  the  ordering  of  the  0 estimates  derived  from  these 
response  patterns  within  each  procedure.  For  maximum  likelihood  scoring  using 
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a normal  ogive  ICC,  the  ordering  of  the  0 estimates  derived  from  the  five 
response  patterns  was  exactly  the  same  as  that  obtained  from  the  Bayesian 
scoring  procedure,  although  the  numerical  values  of  the  0 estimates  were 
uniforml v higher  for  the  Bayesian  procedure.  For  both  these  scoring  methods, 
there  was  a tendency  for  higher  ability  estimates  to  be  obtained  when  a more 
difficult  item  was  answered  correctly.  For  example,  the  lowest  0 estimate  was 
obtained  by  both  scoring  methods  when  the  easiest  item  (Item  1)  was  answered 
correctly  (Response  Pattern  17);  when  only  Item  2 was  answered  correctly 
(Response  Pattern  9) , the  0 estimates  from  both  the  Bayesian  and  maximum 
likelihood  normal  procedures  increased.  In  addition,  both  scoring  methods  took 
into  account  the  discriminations  of  the  items  involved.  For  example.  Response 


Scores  Given 

Table  4 

to  Each  Response  Pattern 

by  Three  Scoring  Methods 

Response 

Maximum 

Likelihood 

Pattern 

Normal 

Logistic 

Bayesian 

1.  00000 

o ok 

co  k 

-1.72 

2.  00001 

-.93 

O 

\D 

1 

-.64 

3.  00010 

-.61 

-1.19 

-.38 

4.  00011 

-.13 

-.46 

.11 

5.  00100 

-1.42 

-1.60 

-1.06 

6.  00101 

-.50 

-.84 

-.28 

7.  00110 

-.30 

-.46 

-.11 

8.  00111 

.13 

. 46 

.30 

9.  01000 

-1.24 

-1.19 

-.89 

10.  01001 

-.23 

-.46 

-.15 

11.  01010 

.03 

.00 

.00 

12.  01011 

.50 

.84 

.41 

13.  01100 

-.60 

-.46 

-.42 

14.  01101 

.23 

.46 

.17 

15.  OHIO 

.39 

.84 

.28 

16.  01111 

.93 

1.60 

.64 

17.  10000 

-1.63 

-1.60 

-1.16 

18.  10001 

-.39 

-.84 

-.24 

19.  10010 

-.17 

-.46 

-.06 

20.  10011 

.30 

.46 

.39 

21.  10100 

-.78 

-.84 

-.58 

22.  10101 

.03 

.00 

.11 

23.  10110 

.17 

.46 

.23 

24.  10111 

.61 

1.19 

.62 

25.  11000 

-.42 

-.46 

-.29 

26.  11001 

.60 

.46 

.51 

27.  11010 

.78 

.84 

.63 

28.  11011 

1.42 

1.60 

1.09 

29.  11100 

.42 

.46 

.31 

30.  11101 

1.24 

1.19 

.93 

31.  11110 

1.63 

1.60 

1.08 

32.  11111 

oo  k 

ook 

1.55 

k 

For  maximum  likelihood  scoring,  it  is  not  possible  to  score 
response  patterns  with  all  correct  or  incorrect  answers. 
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Pattem  2 (with  a correct  response  to  Item  5,  the  most  difficult  item)  was 
assigned  higher  scores  than  Pattern  5;  but  Pattern  2 was  assigned  lower  scores 
than  Pattern  3 (which  had  a correct  response  to  Item  A,  the  second  most  diffi- 
cult item),  since  in  Response  Pattern  3 a correct  answer  was  given  to  an  item 
(Item  4)  with  a higher  discrimination  than  that  of  Pattern  2 (Item  5). 

On  the  other  hand,  assuming  a logistic  ogive  ICC  for  the  maximum  likelihood 
scoring  procedure,  estimated  values  of  0 were  related  to  the  discriminations  of 
the  items  answered  correctly.  Those  response  patterns  for  which  the  discrimi- 
nations of  the  items  answered  correctly  were  the  same  were  assigned  the  same 
score,  namely  -1.60  for  Patterns  2,  5,  and  17  and  -1.19  for  Patterns  3 and  9. 

For  the  latter  two  response  patterns,  the  discriminations  of  the  items  answered 
correctly  were  1.50;  for  the  former  three  response  patterns,  they  were  1.00. 
Thus,  the  magnitude  of  the  scores  was  a function  of  the  item  discriminations, 
and  the  item  difficulties  did  not  affect  the  0 estimates. 

These  data  indicate  that  the  assumption  of  different  forms  of  the  ICC 
within  the  maximum  likelihood  scoring  procedure  will,  in  general,  result  in 
different  0 estimates.  Since  the  Bayesian  0 estimates  were  different  from 
both  the  maximum  likelihood  estimates,  these  three  ICC-based  scoring  procedures 
are  not  interchangeable.  However,  additional  research  is  required  to  further 
delineate  the  similarities  and  differences  among  the  0 estimates  derived  by 
different  ICC-based  scoring  procedures  and,  more  importantly,  to  assess  the 
implications  of  these  differences  in  practical  applications. 


General  Description  of  the  Programs 

This  report  describes  three  computer  programs  for  scoring  test  data  with 
ICC  models — LINDSCO,  ADADSCO,  and  LINPSCO.  Table  5 summarizes  the  major 
features  of  these  programs.  LINDSCO  (LINear  Dichotomous  SCOring)  is  designed 


Table  5 

Summary  of  Program  Capabilities 


Model  and 

Dichotomous 

Polychotomous 

Scoring 

Linear 

Adaptive 

(LINPSCO) 

Procedure 

(LINDSCO) 

(ADADSCO) 

Graded 

Nominal 

Logistic  Ogive 

Bayes iana 
Maximum 

NO 

NO 

NO 

NO 

Likelihood 

YES 

YES 

YES 

YES 

Normal  Ogive 

Bayes iana 
Maximum 

YES 

YES 

NO 

NO 

Likelihood 

YES 

YES 

YES 

NO 

aThe  Bayesian  scoring  procedure  is  based  on  Owen  (1975). 


to  be  used  for  scoring  test  data  for  conventional  (linear)  tests  in  which  all 
items  are  administered  to  each  testee.  It  requires  responses  to  be  dichotomous; 
that  is,  responses  are  scored  into  one  of  two  categories,  such  as  "correct" 
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I and  "Incorrect."  Omissions  are  permitted,  but  they  are  Ignored  In  the 

computations.  The  number  of  omitted  Items  Is  tallied  from  the  number  of 
Items  administered  and  reported  as  part  of  the  output  for  each  testee.  Either 
the  normal  or  logistic  ogive  response  model  can  be  used  with  ICCs  described  by 
one,  two,  or  three  parameters  for  maximum  likelihood  scoring.  Response  patterns 
may  also  be  scored  by  Owen's  (1975)  Bayesian  method  which  assumes  a normal 
ogive  ICC.  The  user  can  also  specify,  in  addition  to  a total  test  score, 
subscores  on  as  many  as  25  subscales. 

ADADSCO  (ADAptive  Dichotomous  SCOrlng)  is  similar  to  LINDSCO,  but  it  is 
designed  specifically  for  scoring  item  response  data  derived  from  adaptive  >. 

testing.  Since  in  adaptive  testing  each  respondent  answers  a different  set 
of  test  items,  the  program  must  locate  for  each  testee  the  item  parameter 
estimates  of  each  attempted  item;  LINDSCO,  in  contrast,  does  the  item 
search  only  once.  ADADSCO  also  differs  from  LINDSCO  in  that  it  has  no 
subscale  scoring  capabilities. 

LINPSCO  (LINear  Polychotomous  SCOring)  is  designed  to  score  data  from 
linear  (conventional)  tests  in  which  each  testee  is  administered  all  items, 
and  items  are  scored  into  more  than  two  categories.  Three  models  are 
available:  the  graded  normal  and  logistic  ogive  models  (Samejima,  1969), 
and  the  nominal  logistic  model  (Bock,  1972).  In  LINPSCO  only  maximum  likeli- 
hood scoring  is  available,  and  subscale  scoring  is  not  possible. 

All  three  programs  compute  both  test  information  and  response  pattern  infor- 
mation values  when  maximum  likelihood  scoring  is  used.  Response  pattern 
information  (Samejima,  1973)  provides  an  estimate  of  the  precision  of  measure- 
ment for  a specified  response  pattern  and  can  be  used  to  compare  the  quality  of 
trait  estimates  derived  from  specific  test  administration  and/or  scoring  pro- 
cedures (e.g.,  Bejar,  Weiss,  & Gialluca,  1977). 

NUMERICAL  PROCEDURES 


Dichotomous  Data 


Maximum  Likelihood 


The  numerical  procedure  for  maximum  likelihood  scoring  of  dichotomous 
data  consists  of  two  stages.  Ip  the  first  stage  an  initial  estimate  is  sought 
by  the  bisection  method.  Once  this  initial  estimate  is  obtained,  it  is  refined 
further  by  the  Newton-Raphson  method. 

The  bisection  routine  begins  in  the  interval  ±5.00.  If  the  sign  of  the 
first  derivative  of  the  likelihood  function  during  the  firs’:  iteration  is  the 
same  when  evaluated  at  5.00  and  at  -5.00,  a value  of  0.0  in  returned  as  the 
initial  estimate  of  0.  Otherwise,  five  additional  iterative  are  performed. 
After  the  sixth  iteration,  the  width  of  the  interval  has  been  reduced  to 
10/(26)  » 10/64  ■ .15.  The  midpoint  of  that  interval  is  the  initial  estimate 
which  is  then  refined  further  by  Newton-Raphson  iterations  of  the  form: 

Vl  - - V//") 


il 


a 


[3] 

j 

j 


1 


where 

9^^  is  the  new  estimate, 

8 is  the  estimate  from  the  last  iteration, 

m 

f is  the  first  derivative  of  the  log-likelihood  function  evaluated  at 

8 , and 
m 

f"  is  the  second  derivative  of  the  log-likelihood  function  evaluated 

at  6 . 
m 

This  iterative  process  is  continued  until  |§  . , - 9 [ < .005.  If  that 

' m+1  m ' 

criterion  has  not  been  met  at  the  end  of  50  iterations,  the  case  is  said  to  be 
nonconvergent. 

Formulas  for  derivatives.  Let  v*  {u  , g~ 1,  2,  ...  n}  be  a response 
vector  such  that  ^ 

| 1 if  the  item  is  answered  correctly 
u = V 

^ ! 0 if  the  item  is  answered  incorrectly. 

Note  that  for  scoring  purposes,  the  response  vector  does  not  include  rejected 
or  omitted  items.  The  probability  that  U =1  for  a given  value  of  9 and  item 

parameters  a„,  b , and  a is  given  by  ^ 

9 9 9 

V9)  = a9 + (1-V[1  + e"1-7a^(0 " Vr1  m 


for  the  logistic  ogive  model  and  by 


Pg(9)  = og  + (1 


„>-L-  f 

9 /2F  X 


Og  + (1-^)  4>  lag(e-bg)  1 


for  the  normal  ogive  model,  where  <J>  stands  for  the  standard  cumulative  normal 
distribution. 

The  log-likelihood  function  for  the  response  vector  is 


V 6)  -JlogP  (9>*fc  Qgwl'u9 


E [m  logP  (9)  + (1-m  >log0  (0)] 
*7  y y y y 


, -.V*  ►*'  • . • 

■ •* 


[6] 
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where  Q ( 0)  - 1-P  (9)  and  P (0)  is  given  by  either  Equation  4 or  5. 

y y y 

In  general,  the  first  and  second  derivatives  of  the  log-likelihood  function 
of  a response  vector  are  given  by 


3LW(0) 

ae~ 


i 

9 


and 


9 2Lv(Q) 

ae7 


e ■ 

9 


7 1 \/3«  <e)\2  / . \ /32G  (0)\-| 

+ a'V[(^6))(-V-) + (ysy)  (-*-)_ 


18] 


For  the  logistic  ogive  model,  after  simplification  and  letting 
x m 1.7a  (0-6 ),  these  expressions  are 

*j  y 


9 Ly(0) 
90 


- -1.7  E 

9 


a e 


1+e 


x -\ 


+1.7  E 

9 


u a e 
_g_g _ 

c + eJ 

g 


[9] 


and 


9 zLy(0) 
907 


-2.89  E 

9 


r 2 X 

ar  & 

" u a1 a ex  ' 
9 9 9 

9 

+ 2.89  E 

9 

. (l+eX)2  . 

J 

[10] 


For  the  normal  ogive  model,  letting  x = -a2(0-fc  )2/2,  the  corresponding 
expressions  are  & ^ 


V9) 

90 


“KV-*/ 

Vu'V*lVe‘V1 


(1-u^)  (2TT)‘1/2a-cff)a£7ex 

l-{e“+(l-e^)*[a  <0-6  ) ] } 


[11] 


and 
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92Ly(6) 

— J U 
90  g g 


[(2tt)-1/2(1  -cjaaex]z 
{[a  -(l-o  )Ha  (Q-b  )]F 

v y %j 


(2TT)"1/V(0-^)(l-Cff)eX 

■{V(1  -f)nag(bg-Q)]} 


[(2TT)"i/2(l -cja^e*]2  (2Tr)'1/2a*(6-feff)(l-off)ex  | 

\AcgHl-cgma^(Q-bg)  ]K  + l-{^+(l-^)-J»[a^(e-/)]}  ) 


[12] 


Computation  of  information.  With  maximum  likelihood  scoring,  two  measures 
of  information  are  computed  for  each  response  pattern.  One  is  response 
pattern  information  (Samejima,  1973)  denoted  by  ?(§);  the  other  is  test  infor- 
mation (Birnbaum,  1968;  Samejima,  1969)  denoted  by  /(§).  Test  information  is 
defined  as  the  expected  value  of  the  second  derivative  of  the  log-likelihood 
function,  i.e.. 


^ | 92logl/y(0)  )■ 

E j 90*  | 


[13] 


Response  pattern  information,  on  the  other  hand,  is  defined  by 


A A 


1(0) 


92logLy(0)- 

a©2 


[iM 


that  is,  the  "observed,"  as  opposed  to  expected,  value  of  the  second  derivative 
of  the  log-likelihood  function  evaluated  at  §. 


These  two  measures  of  information  will  be  the  same  for  models  in  which 
there  is  a sufficient  statistic  for  the  response  vector.  In  particular, 
this  is  true  in  the  one-  and  two-parameter  logistic  ogive  models.  It  is 
also  true  for  the  "zero"  parameter  normal  ogive  model,  i.e.,  when  the  items 
are  parallel.  The  value  of  ?(0)  for  a given  response  pattern  is  simply  the 
value  of  the  second  derivative  of  the  log-likelihood  function  at  the  last 
iteration,  i.e.,  evaluated  at  the  estimated  value  of  0. 


/(§)  is  computed  by 


1(0) 


{p:(0)>2 

1 — 2— — 

9 P (0){1.O-P  (0)} 

g g 


[15] 
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where  P (0)  is  given,  in  general,  by  Equation  4 for  the  logistic  ogive  model  and 

y 

by  Equation  5 for  the  normal  ogive  model. 

For  the  normal  ogive  model. 


a_(l-c_) 


P'( §)  - -£L:-  [e~al[Q~h9]212}  ; 


and  for  the  logistic  ogive  model. 


ye> 


. , .1.7 a (0-b  ) 

l.7^(l-gg)g  gK  g 

u^-V^Vi2 


[16] 
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Standard  error.  The  standard  error  of  measurement  associated  with  0 is 

computed  as  l//£(0) , that^is,  the  reciprocal  square  root  of  response  pattern 
information  evaluated  at  0. 


Bayesian  Scoring 

The  Bayesian  scoring  procedure  used  by  LINDSCO  and  ADADSCO  is  derived  from 
Owen's  (1975)  sequential  adaptive  testing  strategy.  However,  since  the 
present  application  assumes  that  the  test  items  have  already  been  administered, 
only  the  scoring  aspect  is  of  interest. 


The  procedure  makes  the  assumption  that  the  prior  distribution  of  0 is 
normal,  with  mean  ^=*0.0  and  variance  Oo=1.00,  where  subscript  0 denotes  the 
fact  that  no  items  have  yet  been  administered.  After  the  w11*1  item  is  adminis- 
tered, the  mean  and  variance  of  the  posterior  density  function  are  computed 

according  to  the  following  equations.  If  the  response  to  the  m + l1"*1  item  is 
correct , 


Following  an  incorrect  answer 


In  Equations  18  through  21 (from  Owen,  1975), 

4> (Z?)  is  the  normal  probability  density  function, 

<J >(D)  is  the  cumulative  normal  distribution  function. 


A * a + (1 -a  ) <&  ( -D ) . 

9 9 


[22] 


[23] 


After  the  last  item  has  been  administered,  the  posterior  mean  is  the  esti- 
mated 0 and  the  posterior  variance  is  a measure  of  the  error  associated  with  that 
estimate.  Because  the  posterior  distribution  after  every  item  is  administered 
is  approximated  by  a normal  distribution  in  this  procedure,  there  is  a certain 
amount  of  inaccuracy  in  the  estimate.  Moreover,  the  resulting  scores  are 
order  dependent  (Sympson,  1977),  i.e.,  if  a response  vector  were  to  be  scored 
after  rearranging  the  items,  the  resulting  0 estimate  would  be  slightly 
different. 


Computation  of  Expected  Proportion  of  Correct  Answers 

The  expected  proportion  of  correct  answers  (EXPTOT)  is  defined  as 


EXPTOT 


E P (0) /NI  , 
9 9 


[24] 


where  P (8)  is  computed  from  Equation  4 for  the  logistic  ogive  model  and  Equation  5 

1/ 

for  normal  ICCs.  NI  is  the  number  of  items  on  which  the  estimate  of  0 is 
based.  EXPTOT  is  simply  an  estimate  of  the  true  score  associated  with  0 
(Lord  & Novick,  1968,  p.  387). 
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Poluahotomous  Data 


LINPSCO  is  capable  of  scoring  polychotomous  data  when  item  parameters 
have  been  estimated  according  to  a graded  model  of  either  normal  or  logistic  ogive 
form  (Samejima,  1969)  or  according  to  Bock's  (1972)  nominal  logistic  model. 

For  the  graded  model,  the  numerical  procedure  consists  of  a bisection  stage 
of  six  iterations  followed  by  Newton-Raphson  iterations.  For  the  nominal 
logistic  model,  the  initial  estimate  obtained  from  the  bisection  stage  is 
refined  further  by  the  secant  method  rather  than  by  Newton-Raphson  iterations. 

In  each  case,  the  bisection  phase  begins  in  the  interval  ±5.00.  During 
the  first  iteration,  if  the  sign  of  the  first  derivative  of  the  log-likelihood 
function  is  the  same  when  evaluated  at  5.00  and  at  -5.00,  a value  of  0.0  is 
returned  as  the  initial  estimate.  Otherwise,  five  additional  iterations  are 
performed.  After  six  iterations,  the  width  of  the  interval  is  reduced  to 
10/(26)=10/64».15.  The  midpoint  of  that  interval  is  taken  as  the  initial 
estimate. 

The  Newton-Raphson  procedure  used  with  the  graded  models  refines  the 
initial  estimate  with  iterations  of  the  form  shown  in  Equation  3.  This 


iterative  procedure  is  continued  until 


- 9 | is  less  than  .005  or  the 


number  of  iterations  is  greater  than  50.  The  secant  procedure  is  similar 
to  Newton-Raphson  iterations,  except  that  f"  in  Equation  3 is  an  approximation 
to  the  second  derivative  of  the  log-likelihood  function  given  by 


(0  - e„  . ) 

m n-l 


Graded  Models 

Let  V m {x  , g~  1,  2,  ...  } be  a response  vector  exclusive  of  omitted  and 
rejected  items  such  that 

(1  = if  the  "best  response  was  given 

2 = if  the  second  "best"  response  was  given 


(m^-1  if  the  next  to  worst  response  was  given 
nig  if  the  worst  response  was  given. 

For  the  graded  logistic  ogive  model,  the  probability  that  x takes  one  of 
the  values  between  1 and  m is  given  in  general  by 


Uxn  , yxn-i  _i 

P„  (0)  - P„  - [1  + e 9)~l  - [1  + e 9 ] 1 


-17- 


where 


V = -a  D(Q-b  ) , 

\ 9 *? 

-i = -%D(Q-bx  * and 


fl  * 1.7  is  a scaling  factor. 


When  x =1, 

9 


_i 

Px  = [1+e  g]  1 . 

9 


When  x = m , 

3 9 


—X  —1 

P - 1 - [1  + e xg  L]  1. 
x 

9 


For  the  graded  normal  ogive  model,  the  probability  that  x takes  one  of 

y 

the  values  between  1 and  m is  given  in  general  by 

y 


V9’ 


■Hp 

•'y*  -1 


(2ir)~  I 


■'*['%]  -’[V1]  , 


where 


w * a (0-b  ) 

Jx  g x ' 
9 9 


yx  -1  = ao(9-bx  -1} 

9*9 


When  x =1, 

9 


Px  “ (2*>' 

9 


y*9  fnit 
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When  xa  m mg* 


i fV1  et2/2dt  • 

J — OO 

The  log-likelihood  function  for  a given  response  vector  is  given  by 


P - 1 - (2tt) 

9 J — OO 


[35] 


Ly(0)  - log  n Px 


*0 


E rx  [log  P ] , 


[361 


where 


l 1 if  the  x th  response  category  is  chosen 

r * \ 0 


I 0 


'g  i u otherw ' se 
The  general  first  derivative  of  the  log-likelihood  function  is 


3L  (0) 

2 . E E r L 

39  *X* 


[37] 


Same j iraa  (1969)  referes  to  L*  as  the  basic  function.  Since  Lx^"(3P*^/3e)(P*^ 


3 L (0)  3P*</30 

y -EE  r_  a- 


•3e  " & K‘xg  ?x- 


[38] 


9 


The  general  second  derivative  of  the  log-likelihood  function  is  given  by 


32l  (0)  , a2  V302  , 

2- — = E E r [-(L  ) + p 1 


[39] 


Specifically,  for  the  graded  logistic  ogive  model, 

3L.(0) 


- - * S V-7{1-^  - 

y 


30  9 xn  9 ~g  g 


[40] 
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32L  (0) 


A m 

Z Z 2.89 a2  j2  [-{(Lx  )2  + UQ*  Q*  >' 

y "o  » x_  L ' g a g g 


i 


-i  - »(Pi  -i  \ -i'!]  , 

g 9 9 J ’ 


where 

P*  = [l  + e~lmla9^~bxg*]~l 

px  -i = [1  + «~1*7a&(0“Vl)r1  , 

9 

Po  = 0 • 

P*  =1  * 

9 

Q*  = 1 - P*  , and 

x„  x 

0 9 

Q*  = 1 - P . 
x -1  x -1 

For  the  graded  normal  ogive  model,  letting  2^ 
corresponding  expressions  are 


•[a2(0  - 6 )2]  / 2,  the 

* 


3£y(0) 


rx  a 3x  V1  \ 

—2~2  [e  9 - e 9 ])/  P„  (6) 

m f xg 


The  second  derivative  is  given  by 


92L.(0) 


s k\  [-*•/]♦ 


~~2  { ( 0-i  ) eXg  - (8  - bx 
m xg  X9 


eV>]  ' 
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zXq- 1 

When  x -1,  e -0,  and  P is  given  by  Equation  34;  when  x , e 9*0 
9 0 9 

and  P^  is  given  by  Equation  35. 


Nominal  Logistic  Model 

For  the  nominal  logistic  model,  the  probability  of  x , given  9,  is  given  by 

9 


P (0)  = P 
x x 

9 9 


- / W,  ? e(<%9+Bs>  , 


[501 


S=1 


where  ag  and  8g  are  the  slope  and  intercept  parameter  for  the  s*"*1  response 
category. 

The  secant  method  requires  only  the  first  derivative  of  the  log-likelihood 
function.  That  derivative  is 


m„ 


— 9 . . ot.<,9+6c 

9L  (0)  A "x  (ax  ^)e 


v 


= E 


9 9 


90  9 ™9  (as0  + 8s) 

E e 
s= 1 


Computation  of  Information 


[51] 


Response  pattern  information  is  computed  as  the  value  of  the  second 
derivative  at  the  last  iteration.  For  the  nominal  logistic  model,  that  value  is 
an  approximation.  Test  information  is  computed  from  the  general  formula  given 
by  Samejima  (1969), 


1(0)  -EE  (9P  /90)2P 

g x x x 

9 9 9 


[52] 


This  expression  involves  only  the  first  derivative  of  the  response  model. 
The  appropriate  expressions  are  listed  below. 


For  the  graded  normal  ogive  model, 

9P 
x 

~ 90 


S' . i [>  . ,V]  , 


/2rr 


where  Z - -[a2(0-i  )2]  / 2 . 

A * X9 


[53] 


For  the  graded  logistic  ogive  model, 


9P 


x 


-55s-  ■ 1.7  a [P*  (1-P*  ) 
99  gxx 


Px  -1  (1'P*  -1)] 

9 9 


[54] 


where  P* 
x 
9 


[i«-i-7a  w-Vr1. 


For  the  nominal  logistic  model. 


9 P 


\e 


<°k0e  + 9xa>  na  <%9  + > 

9 9 TT  e 9 9 (a  -<*s) 

s=l  ' g 


99 


j?  e(0V  + V 


V* 

s-1 
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PSP  OF  THE  PROGRAMS 


Input 


For  each  of  the  programs,  three  types  of  input  are  required: 

1.  The  Program 7 Parameter's,  which  consist  of  specifications  as  to  the 
number  of  items  in  the  pool,  the  options  chosen,  the  scoring  key, 
and  so  forth. 

2.  The  Item  Pool,  which  contains  the  item  pirameter  estimates  on  as 
many  as  600  items  for  LINDSCO  and  ADADSCO,  and  100  items  for  LINPSCO. 

3.  The  Test  Response  Data  consists  of  tcstee  name  and  identification 
number  and  each  testee's  item  responses.  For  LINDSCO,  item  responses 
need  not  be  dichotomized  beforehand;  for  ADADSCO,  they  must  be  dichot- 
omized unless  a key  is  provided  as  part  of  the  item  pool.  For 
ADADSCO,  the  number  of  items  attempted  and  the  identification 

number  of  each  item  attempted  must  also  be  supplied  as  part  of 
the  test  response  data.  For  LINPSCO,  the  test  response  data  must 
be  supplied  in  such  a way  that  the  first  category  corresponds  to 
the  "best"  response,  while  the  last  category  corresponds  to  the  "worst" 
response,  based  on  previously  obtained  item  parameterization  data. 

Testee  response  data  containing  all  correct  or  incorrect  answers  cannot 
be  scored  by  maximum  likelihood.  If  such  a response  pattern  is  found,  a 
message  is  printed,  and  the  estimated  9 is  set  to  10.00  if  all  responses  are 
correct  and  to  -10.00  if  all  responses  are  incorrect.  The  information  is 
set  to  0.0  in  both  cases.  Response  patterns  with  all  answers  correct  or 
incorrect  present  no  problem  for  Bayesian  scoring,  and  they  are  processed 
normally;  however,  a message  is  still  printed.  Appendix  B gives  examples  of 
the  use  of  each  of  these  programs. 


Table  6 

Input  Program  Parameters  for  LINDSCO,  ADADSCO.and  LINPSCO:  Card  Set 
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Program  Parameters 

Table  6 describes  the  input  program  parameters  for  all  three  programs,  using 
Card  Set  1 (all  numeric  information  is  right  justified).  After  Card  Set  1,  the 
program  parameter  and  input  for  each  of  the  three  programs  differs,  as  indicated 
below. 


LINDSCO  (Card  Set  2-10). 
Card  Set  2 (8A10). 


Card  Set  3 (1615). 


Card  Set  4 (8011). 


Card  Set  5 (8011). 


Card  Set  6 (8A10). 


Card  Set  7 (8A1C) . 


Card  8 (15). 


Card  Set  9 (215). 
(Omit  if  the  number 
of  subscales  is  0. ) 


Card  Set  10  (1615). 
(Omit  if  the  number 
of  subscales  is  0. ) 


The  variable  format  for  the  item  pool  is  punched 
on  this  card,  using  I-fields  (see  Item  Pool  below). 

Punch  in  five-column  fields  the  item  identification 
number  of  the  items  in  the  test  in  the  same  order 
in  which  they  appear  in  the  test.  Continue  on  as 
many  cards  as  necessary. 

A "1"  in  a given  column  is  punched  to  omit  a 
specified  item  from  all  computations,  e.g.,  if  the 
10th  item  is  to  be  omitted,  punch  "1"  in  column  10; 
if  the  100th  item  is  to  be  omitted,  punch  "1"  in 
column  20  of  the  second  card.  Continue  on  as  many 
cards  as  necessary. 

This  card  contains  the  scoring  key  for  the  test.  In 
general,  the  nth  column  contains  the  key  for  the  nth 
item,  as  in  Card  Set  4.  Continue  on  as  many  cards 
as  necessary. 

Variable  format  for  reading  the  subject  information 
and  test  response  data  (see  Test  Response  Data  below 
for  field  type  specifications). 

The  description  of  the  run  is  written  on  three  cards. 
The  three  cards  must  be  included  even  if  they 
are  blank. 

Punch  the  number  of  sub scales  to  be  scored  in 
columns  1-5,  maximum  is  25.  If  no  subscales  are 
to  be  scored,  punch  "0"  in  column  5;  in  that  case, 
this  is  the  last  card  set. 

For  each  scale,  punch  the  following  information: 
Columns  1-5:  Number  of  items  In  subscale  (maximum 
is  60)  . 

Columns  6-10:  Scale  number.  Repeat  for  each 

subscale  beginning  on  a new  card. 

Punch  in  five-column  fields  the  item  identification 
number  of  the  items  in  the  subscales.  Continue  on 
as  many  cards  as  necessary.  Repeat  for  each 
subscale,  beginning  on  a new  card  for  each  subscale. 
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ADADSCO  (Card  Set  2-5). 
Card  Set  2 (8A10) . 


Card  Set  3 (1615). 


Card  Set  4 (8A10). 


Card  Set  5 (8A10). 


Variable  format  for  item  pool,  using  I-fields. 

It  must  be  contained  on  one  card  (see  Item  Pool 
below) . 

Columns  1-5:  The  number  of  items  to  be  omitted, 

i.e.,  excluded  from  the  computations. 

If  none,  punch  "0"  in  column  5. 

Columns  6-10  and  subsequent  five-column  fields: 

The  item  identification  numbers  of 
items  to  be  omitted.  Continue  on  as 
many  cards  as  necessary.  If  more  than 
one  card  is  necessary,  begin  punching 
on  the  second  card  in  columns  1-5. 

Variable  input  format  for  reading  subject  information 
and  test  response  data.  It  must  be  contained  on  one 
card  (see  Test  Response  Data  below  for  field  type 
specifications) . 

Description  of  the  run  is  written  on  three  cards. 

These  cards  are  required,  even  if  they  are  left  blank. 


LINPSCO  (Card  Set  2-7). 
Card  Set  2 (8A10). 

Card  Set  3 (8011'). 

Card  Set  4 (1615). 

Card  Set  5 (8011). 


Card  Set  6 (8A10). 


Variable  format  for  the  item  pool.  It  must  be 
contained  on  one  card  (see  Item  Pool  below  for 
field  type  specifications). 

Punch  in  the  column  the  number  of  response 
categories  minus  1 for  the  n1^1  item.  Continue 
on  as  many  cards  as  necessary. 

Punch  in  five-column  fields  the  item  identification 
numbers  of  the  items  in  the  test.  The  numbers  must 
appear  in  the  same  order  as  the  items  appear  in 
the  test.  Continue  on  as  many  cards  as  necessary. 

The  information  on  this  card  is  used  to  omit 
specified  items  from  the  computations.  To  omit 
the  nth  item,  punch  a "1"  in  the  nth  column  of 
this  card;  otherwise,  punch  "O.”  If  no  items  are 
to  be  omitted,  punch  as  many  zeros  as  there  are 
items  in  the  test.  Continue  on  as  many  cards  as 
necessary. 

Variable  format  for  subject  information  and  test 
response  data.  It  must  be  contained  on  one  card 
(see  Test  Response  Data  below  for  field  type 
specifications) . 

Description  of  the  run  is  written  on  three  cards. 


Card  Set  7 (8A10) . 


1 
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Item  Pool 

LINDSCO  and  ADADSCO.  To  score  the  response  data,  a file  containing  the 
item  pool  item  parameter  estimates  must  be  prepared  beforehand  and  placed  in 
a file  called  IPOOL.  The  file  consists  of  a line  for  each  item  in  the 
pool  with  the  following  information: 

1.  A unique  item  number; 

2.  Estimate  of  the  a parameter; 

3.  Estimate  of  the  b parameter; 

4.  Estimate  of  the  a parameter; 

5.  Correct  alternative  for  this  item,  i.e.,  the  keyed  response. 

For  LINDSCO,  only  Items  1 through  4 must  be  supplied;  for  ADADSCO,  Item  5 must 
be  supplied  also,  although  it  could  be  a "dummy"  key  (e.g.,  a blank),  since  the 
data  may  already  be  scored  (see  columns  44-45  for  Card  Set  1). 

The  exact  format  of  this  information  is  not  critical,  since  it  is  read  with 
a user-specified  variable  format.  However,  the  following  limitations  must  be 
observed:  (1)  the  information  must  be  read  in  the  above  order;  (2)  the  item 

number  must  be  read  in  integer  mode;  (3)  the  item  parameter  estimates  must  be 
read  in  floating  point;  and  (4)  the  key,  if  ADADSCO  is  being  used,  must  be  read 
in  integer  mode. 

A typical  format  for  LINDSCO  could  be 

(10X,I4,3F10. 2)  . 

For  ADADSCO,  a typical  format  might  be 

(10X, 14 ,3F10. 2, 12)  . 

All  three  parameter  estimates  must  be  read  even  if  the  user  is  using  a 
one-  or  two-parameter  model.  This  presents  no  difficulties,  however,  since 
in  the  case  of,  say,  a two-parameter  model,  the  third  parameter  is  0 for  all 
items.  This  may  be  accomplished  by  reading  blanks  or  zeros,  or  by  editing 
item  parameter  estimates  (see  below). 

The  number  of  items  in  the  pool  may  range  from  the  number  of  items  in 
the  test,  M,  to  600.  If  the  item  pool  for  LINDSCO  consists  of  only  the  items 
being  scored  in  the  test,  then  0PT2  should  be  set  to  1.  This  indicates  to 
the  program  that  items  do  not  have  to  be  searched.  On  the  other  hand,  if 
the  pool  consists  of  items  in  addition  to  those  used  in  the  present  test, 
then  0PT2  should  be  set  to  0.  This  instructs  the  program  to  search  for  the 
item  and  to  retrieve  the  corresponding  item  parameters.  For  both  LINDSCO  and 
ADADSCO,  if  at  least  one  of  the  items  being  called  for  is  not  found  in  the  pool, 
the  program  prints  a message;  and  the  unavailable  item  is  treated  as  an 
omitted  item. 

Editing  of  item  parameter  estimates . LINDSCO  and  ADADSCO  have  several 
options  to  edit  item  parameter  estimates.  If  0PT3-1,  the  program  checks  that 
the  item  parameter  estimates  are  within  certain  bounds.  For  the  discrimination 
(a)  parameter,  the  program  checks  to  see  if  the  estimate  exceeds  AMAX;  if  it 
does,  it  is  set  to  AMAX.  For  the  difficulty  (b)  parameter,  if  the  estimate 
is  below  BMIN,  it  is  set  to  BMIN;  if  it  is  above  BMAX,  it  is  set  to  BMAX.  For 
the  "guessing"  (a)  parameter,  the  program  checks  to  see  if  the  estimate  exceeds 
CMAX;  if  it  does,  it  is  set  to  CMAX.  If  the  user  wants  to  edit  only  one  or  two 
parameters,  the  limits  of  the  other  parameters  should  be  chosen  so  that  the 
editing  has  no  effect. 
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A more  radical  form  of  editing  is  also  possible.  If  OPT3=2 , then  in 
addition  to  the  editing  caused  by  0PT3=1,  the  program  sets  all  a parameter 
estimates  to  CMAX.  If  CMAX=0.0,  this  implies  that  a two-parameter  model  is  in 
effect.  If  OPT3=3,  then  in  addition  to  the  editing  caused  by  0PT3=1  and 
OPT3=2,  the  program  sets  all  a parameter  estimates  to  AMAX. 

LINPSCO.  For  polychotomous  scoring,  the  item  pool  consists  of  the  following 
information  for  graded  normal  and  logistic  ogive  models: 

1.  A unique  item  identification  number; 

2.  The  "discrimination"  parameter,  which  is  common  to  all  response 
categories; 

3.  m -1  "difficulty"  parameters,  where  is  the  number  of  response 
categories  in  the  item.  Since  m can  be  at  most  10,  there  would 

a 

be  at  most  9 difficulty  parameters. 

The  exact  format  for  reading  this  information  is  not  crucial,  since  it  is  read 
by  a user-supplied  format  statement.  However,  the  following  restrictions  must 
be  observed:  (1)  the  identification  number  is  read  first,  in  integer  mode; 

(2)  next,  the  estimated  discrimination  parameter  is  read  in  floating  point 
mode;  (3)  the  ntg- 1 "difficulty"  parameters  are  read  next,  with  the  difficulty 

of  the  best  alternative  followed  by  the  second  best  alternative,  and  so  forth. 

Since  the  program  allows  the  number  of  categories  to  differ  from  item  to 
item,  the  format  should  be  specified  so  that  it  can  read  the  information 
for  the  item  with  the  most  response  categories.  For  example,  if  in  a given 
test,  the  maximum  number  of  response  categories  is  seven,  then  there  should 
be  at  most  six  difficulty  parameters.  The  format  for  such  pools  might  be 
as  follows: 

(I4,6X,F5. 2/10X.6F5. 2)  . 

In  this  format  the  item  identification  number  is  read  in  the  14  field;  the 
discrimination  parameter  is  read  next  in  format  F5.2;  and  the  six  difficulty 
values  are  read  from  the  next  card,  beginning  in  column  11. 

For  the  nominal  logistic  model,  the  item  information  is  read  in  the  following 
order: 

1.  A unique  item  identification  number; 

2.  nig  "slope"  parameters;  and 

3.  rrig  intercept  parameters. 

r 

Differing  from  the  graded  models,  in  the  nominal  model  there  is  a pair  of 
parameters  (a  slope  and  an  intercept)  associated  with  each  response  category. 

Since  the  response  categories  are  not  ordered  in  the  nominal  model,  the  order 
in  which  the  parameters  are  read  is  unimportant.  However,  the  ordinal  position 
in  which  the  parameters  appear  in  the  pool  must  correspond  with  the  integer, 
associated  with  that  response  category.  As  in  the  graded  models,  the  format 
should  be  able  to  read  the  information  for  the  item  with  most  response  categories. 
For  example,  if  the  maximum  number  of  response  categories  is  five,  the  format 
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could  be 

(I4,16X,5F5.2,5X,5F5.2)  . 

In  this  format,  the  item  identification  number  is  read  in  the  14  field;  next, 
the  five  slope  parameters  are  read  in  5F5.2;  and  finally,  the  five  intercept 
parameters  are  read  in  the  last  set  of  5F5.2  fields. 

Test  Response  Data 

Data  for  all  testees  must  be  on  a file  called  DATA.  The  structure  of  this 
file  differs  slightly  for  each  of  the  programs.  In  all  cases,  however,  the 
last  record  of  DATA  must  be  an  end-of-record  marker. 

LINDSCO.  This  program  requires  that  for  each  individual  the  following 
information  be  provided  on  DATA: 

1.  Name, 

2.  Identification  number,  and 

3.  Responses  to  the  test  items. 

The  exact  format  of  this  information  is  not  critical,  since  it  is  read  with  a 
user-supplied  variable  format;  but  the  information  must  appear  in  the  above 
order.  Two  words  are  used  for  testee  name;  thus,  name  should  be  read  with  two 
alphanumeric  words,  e.g.,  2A10.  This  allows  for  up  to  20  characters. 

The  testee  identification  is  read  with  an  alphanumeric  field  of  at  least  1 
column,  e.g.,  Al,  A9.  Test  item  responses  are  read  with  an  integer  format, 
e.g.,  2011. 

The  test  data  may  be  raw  item  responses  (l.e.,  the  number  of  the  alternatives 
chosen)  or  scored  (i.e.,  0 for  incorrect  and  1 for  correct).  However,  in  either 
case,  a scoring  key  must  be  provided  (see  Card  Set  5 for  LINDSCO).  The  key  will 
contain  the  number  of  the  correct  alternative  if  raw  data  are  read.  If  the 
data  are  already  scored,  a "dummy"  key  full  of  "l's"  must  be  provided. 

Omitted  items  are  indicated  by  the  integer  IOMIT  (see  columns  46  and  47  of 
Card  Set  1 for  LINDSCO).  For  raw  data,  this  will  normally  be  an  integer 
greater  than  the  number  of  alternatives.  Similarly,  for  scored  (0-1)  data 
IOMIT  must  be  an  integer  greater  than  1. 

ADADSCO.  The  program  requires  that  the  following  information  be  provided 
on  DATA  for  each  individual: 

1 . Name , 

2.  Identification  number, 

3.  Number  of  items  answered  by  the  testee  (i.e.,  number  of  items  attempted), 

4.  Item  identification  numbers  of  items  attempted,  and 

5.  Responses  to  the  test  items. 

This  information  is  read  in  the  above  order  with  a user-supplied  variable  format; 
thus,  the  exact  format  is  not  critical.  However,  the  following  limitations  must 
be  observed.  Even  though  the  number  of  items  administered  usually  varies  across 
individuals  in  an  adaptive  test,  this  program  assumes  that  the  data  record  for 
each  testee  is  formally  the  same  (i.e.,  that  there  is  the  same  number  of  data 
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lines  per  testee  and  that  these  lines  contain  similar  information).  Thus, 
if  the  maximum  number  of  items  taken  by  anyone  is  MMAX  (see  Card  Set  1 for 
ADADSCO) , but  any  particular  testee  takes  M items,  where  M<MMAX,  then  that 
testee' s record  should  be  "padded"  to  MMAX  items.  This  can  be  accomplished  by 
leaving  an  appropriate  number  of  blank  fields.  The  name  is  dimensioned  for  two 
words  so  the  format  should  allow  for  two  words,  e.g.,  2A10.  The  identification 
number  is  read  with  an  alphanumeric  format,  e.g.,  A8.  The  number  of  items 
is  read  in  integer  mode.  The  item  identification  numbers  and  item  responses 
are  also  read  in  integer  mode.  Note  that  in  reading  the  item  identification 
numbers  and  the  item  responses,  the  format  should  read  MMAX  of  each,  even 
if  some  of  these  will  be  blank  for  a given  individual. 

As  an  example  assume  that  MMAX  was  25;  then  the  variable  format  could  be 
(2A10,Al0,I2/20I4/5I4/25Il) . 

In  this  format,  the  name,  testee  identification,  and  M are  read  from  the  first 
card;  the  item  identification  numbers  are  read  from  the  next  two  cards;  and 
finally,  the  item  responses  are  read  from  the  fourth  card.  Note  that  for 
testees  attempting  20  items  or  less,  the  third  card  will  be  blank. 

The  item  responses  may  be  scored  or  raw  data.  For  scored  data,  the 
responses  have  been  reduced  to  three  categories:  correct,  incorrect, 
and  omitted.  In  this  case,  IFLAG  should  be  set  to  the  integer  corresponding  to 
the  correct  code,  and  IOMIT  should  be  set  to  the  code  for  omitted  responses. 

Note  that  if  IFLAOO,  the  program  ignores  the  key  read  as  part  of  the  item  pool. 
For  raw  data,  the  key  will  have  been  read  as  part  of  the  item  pool;  IFLAC  must 
therefore  be  set  to  0.  IOMIT  will  still  be  operational,  however;  and  it  must 
be  set  to  an  integer  other  than  the  highest  numbered  response  alternative . 

LINPSCO.  The  DATA  file  is  similar  to  LINDSCO's  with  the  exception  that 
the  item  responses  must  include  the  response  category  chosen  by  the  testee  for 
a given  item.  For  graded  models,  the  convention  that  the  best  response  category 
be  coded  "1,"  second  best  "2,"  and  so  forth,  must  be  obeyed.  For  the  nominal 
logistic  model,  this  convention  does  not  apply;  but  care  must  be  taken  so  that 
a category's  response  code  matches  the  ordinal  position  of  that  category  in 
the  IPOOL  file.  For  either  graded  or  nominal  data,  the  code  for  omitted  responses 
should  be  an  integer  greater  than  the  maximum  number  of  response  categories. 


Output 

Four  kinds  of  output  are  produced  by  each  program:  program  parameters, 
item  parameters,  computational  messages,  and  testee  data. 

Program  Parameters 

The  output  consists  of  the  information  in  Card  Sot  1,  the  description  of 
the  run,  and  the  variable  formats  for  reading  the  item  pool  and  the  testee 's 
raw  data. 

Item  Parameters 

LINPSCO  and  LINPSCO.  The  output  consists  of  item  identification  number, 
scoring  key,  rejection  key  (i.e.,  whether  or  not  the  item  was  included  in  the 
computations),  and  the  item  parameter  estimates.  If  the  estimates  have  been 
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eclited,  the  edited  values  will  be  printed.  An  option  (see  column  10  of  Card 
Set  1)  permits  all  of  this  information  to  be  punched  as  well.  If  subscale 
scoring  has  been  requested,  the  item  identification  number  of  the  items  in  each 
subscale  will  be  printed. 

ADADSCO.  The  user  has  the  option,  but  only  for  the  first  10  testees,  to 
print  the  following:  testee's  name  and  identification  number;  and  for  each  item 
attempted,  the  item  identification  number,  the  response  to  that  item,  and  the 
item  parameter  estimates. 

Computational  Messages 

The  program  will  print  a testee's  name  and  Identification  number  if  (1) 
a response  pattern  is  found  with  all  items  correct  or  incorrect,  excluding 
omitted  or  rejected  items;  (2)  a zero  score  has  been  obtained;  or  (3)  it 
was  not  possible  to  achieve  convergence  in  scoring  the  testee's  responses.  For 
polychotomous  data,  a perfect  or  zero  vector  occurs  if  the  testee  responds 
with  the  best  or  worst  response  categories  in  all  attempted  items,  exclusive  of 
omitted  or  rejected  items.  If  an  item  is  not  found  in  the  pool  or  has  extreme 
parameter  estimates,  an  informative  message  is  printed. 

The  number  of  testees  read  and  the  number  of  convergence  failures  are  also 
printed.  If  Bayesian  scoring  has  been  requested,  the  number  of  nonconvergent 
cases  will  be  zero. 

Testee  Data 

LINDSCO.  For  each  testee,  the  following  information  is  written  on  a file 
called  TAPE3: 

1 . Name ; 

2.  Testee  identification  number; 

3.  Scale  number,  or  in  the  case  of  total  score,  a "T"; 

4.  Proportion  of  items  answered  correctly; 

5.  Maximum  likelihood  or  Bayesian  estimate  of  0; 

6.  The  response  pattern  information  for  maximum  likelihood  scoring  or 
the  posterior  variance  of  0 for  Bayesian  scoring; 

7.  The  number  of  items  used  in  the  estimation  of  0,  excluding  items 
rejected,  omitted,  or  not  found; 

8.  The  test  information  associated  with  the  estimated  0 (for  Bayesian 
scoring,  the  information  is  computed  using  the  normal  ogive  model); 

9.  The  true  score  corresponding  to  the  estimated  0; 

10.  For  maximum  likelihood  scoring, 

a.  The  number  of  Newton-Raphson  iterations  needed  to  achieve 
convergence  and 

b.  The  standard  error  of  0. 

The  format  used  for  writing  this  information  for  total  scores  is 

(X, 2A10,A9,*T* , F5. 2 ,2F7 . 2 , 14 , 2F7. 2 , 14 ,F7. 2)  . 

The  subscale  results  are  written  with 

(X, 2A10 , A9 , 12 , F5 . 2 , 2F7 . 2 , 14 , 2F7 . 2 , 14 , F7 . 2 ) . 


ADADSCO.  The  same  information  is  written  as  that  for  LINDSCO  with  the 
exception  of  the  scale  number.  The  format  is 
(X,2A10,A9,F5. 2 ,2F7. 2,I4,2F7. 2 ,I4,F7. 2)  . 

LINPSCO.  For  LINPSCO,  the  following  information  is  written: 

1 . Name ; 

2.  Testee  identification  number; 

3.  Proportion  of  "best"  responses; 

4.  Maximum  likelihood  estimate  of  9; 

5.  The  response  pattern  information; 

6.  The  number  of  items  used  in  the  estimation  of  0 excluding  items 
rejected,  omitted,  or  not  found; 

7.  The  number  of  iterations  needed  to  achieve  convergence. 

8.  The  test  information  associated  with  the  estimated  0; 

9.  Estimated  standard  error  of  measurement. 


AVAILABILITY 


FORTRAN  source  code  listings  of  the  three  programs  are  in  Appendix  C 
(LINDSCO),  Appendix  D (ADADSCO),  and  Appendix  E (LINPSCO).  Copies  of  the 
FORTRAN  source  code  are  available  on  cards  or  tape  at  nominal  cost  from 

Psychometric  Methods  Program 
Department  of  Psychology 
University  of  Minnesota 
75  East  River  Ro?-* 

Minneapolis,  Minnesota  55455 

Telephone:  612-376-7378 

Potential  users  of  these  programs  should  note  that  the  programs  were  written  for 
Control  Data  Corporation  CYBER  series  computers.  Because  of  the  large  word 
size  of  the  CYBER  computers,  accurate  computation  on  other  computers  may  require 
the  use  of  double-precision  arithmetic.  Minimal  additional  modifications  required 
may  include  (1)  modification  of  A10  fields  to  smaller  sizes  used  by  other 
computers  and  (2)  modification  of  FORTRAN  statements  unique  to  the  CYBER  series 
computers. 
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Item  Parameter  Estimation  Programs 
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Appendix  B: 

Examples  of  Program  Use 


The  following  examples  serve  to  illustrate  the  use  of  each  of  the  three 
programs.  These  results  should  also  be  useful  in  testing  the  accuracy  of  the 
results  of  the  programs  in  different  installations. 

LINDSCO 

The  IPOOL  file  for  these  examples  was 


1 

1.000 

-2.000 

.25 

2 

1.500 

-1.000 

.25 

3 

1.000 

0.000 

.25 

4 

1.500 

1.000 

.25 

5 

1.000 

2.000 

.25 

The  DATA  file  is  also  shown  below.  The  first  field  contains  the  names;  the 
second,  the  subject  identification;  and  the  third,  the  response  patterns. 


Name 

I.D. 

Responses 

A00 

1 

ooooo 

A01 

2 

oooui 

A02 

3 

00010 

A03 

4 

oooii 

A04 

5 

ooioo 

A05 

6 

00101 

A06 

7 

00110 

A07 

8 

00111 

A06 

9 

oiooo 

A09 

10 

oiool 

All) 

11 

01010 

All 

12 

Oloil 

A12 

13 

0Ud0 

A13 

14 

Oliol 

A14 

15 

OHIO 

A15 

16 

ouu 

Alb 

17 

10000 

A17 

18 

10001 

A18 

19 

10010 

A19 

20 

100U 

A20 

21 

10100 

A21 

22 

10101 

A22 

23 

10110 

A23 

24 

10111 

A24 

25 

llooo 

A25 

26 

liool 

A2b 
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Example  1.  This  example  illustrates  the  use  of  the  normal  ogive  model 
(0PT4-1)  for  a five-item  test  (m*5)  with  a pool  containing  five  items  (INUP*5). 
The  example  also  illustrates  the  use  of  parameter  editing  (0PT3-1)  in 
which  BMAX“BMIN“0 . 0 , which  in  effect  sets  all  b parameter  estimates  to  0.0. 
AMAX-2.00,  which  means  that  if  there  were  a parameter  estimates  greater  than 
2,  they  would  be  set  to  2.00.  CMAX«.10,  which  means  that  any  a parameter 

estimates  greater  than  .10  will  be  edited  to  .10.  This  example  also  illustrates 
the  use  of  subscales.  The  program  parameter  cards  for  this  example  were 


5 5 111  0.00  1.00  2.00  0,00  0.00  .10  4444 

(9X, I1»3F10.3> 

12  3 4 5 

00000 
11111 

(2A4'A2'10X'5ll> 

RUNS  BASED  ON  ALL  POSSIBLE  RESPONSE  VECTORS  FOR  A FIVE  ITEM  TEST 


1 

3 1 

2 3 5 


The  output  corresponding  to  this  example  is  shown  on  the  following 
pages. 
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Example  2.  This  example  is  identical  to  Example  1 except  that  the 
Bayesian  scoring  routine  was  used  (OPT4=3)  instead  of  the  maximum  likelihood 
normal  ogive.  Only  the  scoring  results  are  shovn. 
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Example  3.  This  example  illustrates  the  use  of  the  maximum  likelihood 
logistic  scoring  routine  (e.g.,  OPT4-2)  without  subscale  scoring.  Only  the 
scoring  results  are  shown. 
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a. 

Z 

C* 

Z 

H 

u 

0) 

u 

u 

o 

CJ 

c 

o 


u 

o 

a 

o 

»- 

a. 

TD 

v 

u 

u 

o> 

c 

X 

w 


w 

c 

o 

•H 


03 

0) 


0 

|4 

0> 

X* 

1 

3 

2 


aj 

c 


g 


(1) 

3 

CO 

03 

0) 


r 


0 
u 
u 
w 

J- 

T3 

C 

03 

aj 

tn 

OJ 

0) 

4-1 

1 


cn 

UJ 


* * * * 

..... 

. . 

♦ * ♦ 

♦ * ♦ 

* * * * 

♦ * 

* * * 

* ♦ * * 

* • 

* * * 

AoO 

1 

T 

0 

-10. OU 

0 

5 

0 

0 

0 

Aul 

i 

T 

•20 

-1 . 0 j 

1.20 

5 

.90 

.20 

2 

• 91 

AU2 

J 

T 

•20 

-99.99 

“99.99 

5 

-99.99 

-99.99 

99 

-99.99 

Au3 

4 

T 

• 40 

-.29 

3.b9 

5 

3.b8 

.42 

3 

.53 

AU4 

J 

T 

• 20 

-1.0b 

1.20 

5 

.90 

.20 

2 

• 91 

Aub 

o 

T 

• 40 

— • 4o 

3.32 

5 

2.95 

.35 

1 

.55 

Au6 

7 

T 

• 40 

-.29 

3.59 

5 

3.68 

.42 

3 

• 53 

Au7 

o 

T 

• u0 

.07 

4.49 

5 

4.46 

.58 

2 

.47 

Auti 

9 

T 

•20 

-99.99 

-99.99 

5 

-99.99 

-99.99 

99 

-99*99 

AU9 

10 

T 

• 40 

-.29 

3.59 

5 

3.68 

.42 

3 

.53 

AiO 

11 

T 

• 40 

-.11 

4.01 

5 

4.26 

.50 

2 

.50 

Ail 

li 

T 

• 60 

• 2b 

4.22 

5 

4.26 

• 66 

2 

.49 

AX2 

13 

T 

• 40 

-.29 

3.59 

5 

3.68 

.42 

3 

.53 

Ai3 

14 

T 

• 60 

.07 

4.49 

5 

4.46 

.58 

2 

.47 

Ai4 

lb 

T 

• b0 

.2b 

4.22 

5 

4.26 

.66 

2 

.49 

Aib 

lo 

T 

• t)0 

.71 

2.74 

5 

2.72 

.83 

2 

.60 

Aib 

17 

T 

•20 

-1.0b 

1.20 

5 

.90 

.20 

2 

.91 

Al7 

lo 

T 

• 40 

— • 4o 

3.32 

5 

2.95 

.35 

1 

.55 

A*ti 

19 

T 

• 40 

-.29 

3.59 

5 

3.68 

.42 

3 

.53 

Ai9 

20 

T 

• o0 

.07 

4.49 

5 

4.46 

.58 

2 

.47 

A *0 

2i 

T 

• 40 

-•4o 

3.32 

5 

2.95 

.35 

1 

• 55 

Ail 

2i 

T 

• bO 

-.10 

4.53 

5 

4.28 

.50 

2 

.47 

Ai2 

23 

T 

• 60 

.07 

4.49 

b 

4.46 

.58 

2 

.47 

Ai3 

24 

T 

• 60 

• 4b 

3.70 

b 

3.70 

.74 

2 

.52 

Ai4 

2b 

T 

• 40 

-.29 

3.59 

5 

3.68 

.42 

3 

.53 

Aib 

2o 

T 

• 60 

.07 

4.49 

5 

4.46 

.58 

2 

.47 

Aib 

27 

T 

• 60 

.2b 

4.22 

b 

4.26 

.66 

2 

.49 

Ai7 

2o 

T 

• 60 

.71 

2.74 

5 

2.72 

.83 

2 

• 60 

Ai8 

29 

T 

• o0 

.07 

4.49 

b 

4.46 

.58 

2 

.47 

Ai9 

30 

T 

• «0 

.4b 

3.70 

5 

3.70 

.74 

2 

.52 

AjO 

31 

T 

• o0 

.71 

2.74 

5 

2.72 

.83 

2 

.60 

A J 1 

3i 

T 

1.00 

10.00 

0 

b 

0 

0 

0 

ADADSCO 


example 

consisted 

of  10  items: 

1 

1.00 

-2.00 

.25  1 

2 

1.25 

-1.50 

.25  1 

3 

1.50 

-1.00 

.25  1 

4 

1.75 

-0.50 

.25  1 

5 

1.00 

0.00 

.25  1 

6 

1.25 

0.50 

.25  1 

7 

1.50 

1.00 

• 25  1 

8 

1.75 

1.50 

.25  1 

9 

1.00 

2.00 

.25  1 

10 

1.25 

2.50 

.25  1 

The  data  for  the  16  subjects  used  in  the  example  are  shown  below: 


loiiioiiSo 
l2345b?8  910 

2 Q 4 

1010 

2  4 6 8 

3 C 5 

11111 

12  3 4 5 

4 D 5 

00000 

678  910 
b £ 8 

11000000 
12366789 

6 F 2 

10 

3 7 

7 6 6 

011111 

12  4 6 8 9 
6 H 9 

101010101 

123456789 
9 I 7 

1100110 
13456  810 
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11  K 8 
10110110 

2 3 4 6 7 0 910 

12  L 9 

lOlllllll 

1234578  910 

13  M 5 

01001 

1 3 5 7 9 

14  N 6 

000001 

4 5 6 7 8 9 

15  P 7 

1011110 

1 2 4 5 6 7 8 

16  Q 6 

011011 

2457  810 


The  program  control  cards  for  this  example  were 


10  10  1012  0.00  1.00  2.00  0.00  0.00  .10  1 3 

(8X.I2.3F 10. 2»12) 

0 

(A2.1X.2A2. X2»/10Il./l0l2) 

description 


In  this  example  the  maximum  number  of  items  attempted  by  anyone  was 
10  (MMAX=10).  Although  the  code  for  omitted  items  was  3 (IOMIT=3),  IFLAG=1, 
which  means  the  key  to  each  item  was  read  from  IPOOL;  however,  in  this  case  it 
was  1 for  all  items.  0PT1=1  means  that  item  information  for  the  first  10 
subjects  will  be  printed.  Editing  of  item  parameters  was  requested  (OPT3=l). 
The  scoring  algorithm  was  maximum  likelihood  logistic. 

The  entire  output  for  this  example  is  shown  on  the  following  pages. 


AQADSCO 


ADAPTIVE  jICHOTOhuS  bCORI  ,6  #ITH  THRel  PARAMETER  MOOELS 


u;iIVEi  SIT*  OF  MINNESOTA 
MPLS.  Mt,\N.  5545b 


XNUP  - 10 

tiMAX  - 10 

lOMIt  = 3 

IFLAo  = 1 

«PT1  = 1 

i,PT2  = 0 

OPT3  = 1 

uPT4  = 2 

lb  = 0 

TSS  - l.nu 

a max  = t . oo 

UMAX  - 0 

oMlN  - 0 

CMAX  = .10 

VARIABLE  TORMAT  For  POOL=<8X.I2»3F10.2.I2) 
vARUoLE  IORNiAT  FOR  DRTA=(A?,1X.2A2.12'/10I1,/1012) 
OESCkIPTIuN 


1 

1 1 

i .do 

0 

.10 

2 

0 

1.25 

0 

• 10 

3 

1 

1.50 

0 

.10 

4 

1 

1.75 

0 

.10 

5 

1 

1.00 

0 

.10 

6 

0 

1.25 

0 

.10 

7 

1 

1.50 

0 

.10 

8 

1 

1.75 

0 

.10 

9 

1 

1.00 

0 

.10 

10 

0 

1.25 

0 

• 10 

2 

8 

2 

1 

1.25 

0 

• 10 

4 

0 

1.75 

0 

• 10 

6 

1 

1.25 

0 

• 10 

8 

0 

1.75 

0 

.10 

3 

C 

1 

1 

1.00 

0 

• 10 

2 

1 

1.25 

0 

.10 

3 

1 

1.50 

0 

• 10 

4 

1 

1.75 

0 

• 10 

5 

1 

1.00 

0 

.10 

bUBJECT  C 

ID: 

4 

U 

6 

0 

1.25 

0 

• 10 

7 

0 

1.50 

0 

• 10 

8 

0 

1.75 

0 

.10 

9 

0 

1.00 

0 

• 10 

10 

0 

1.25 

0 

• 10 

HAS  AtL  ANSWERS  RIGHT 


..  , ,1  . 


OIZirXLHicniPionitM  Testee  Name 
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jubjlct  a 

1U=  4 HAS  NO  right  ANSMEKS 

5 

e 

1 

1 

1.00 

0 

>10 

2 

1 

1.2b 

0 

• 10 

3 

0 

l.bO 

0 

• 10 

b 

0 

1.00 

0 

• 10 

6 

n 

1.2b 

0 

.10 

7 

0 

l.bO 

0 

• 10 

8 

0 

1.7b 

0 

• M 

9 

0 

1.00 

0 

‘•10 

b 

F 

3 

1 

l.bO 

0 

.10 

7 

0 

l.bO 

0 

• 10 

7 

0 

1 

0 

1.00 

0 

• 10 

2 

1 

1.2b 

0 

• 1(1 

4 

I 

1.7b 

0 

’•10 

b 

1 

1.2b 

0 

• 10 

A 

1 

1.7b 

0 

• 10 

9 

1 

1.00 

0 

• 10 

8 

H 

1 

1 

1.00 

0 

• 10 

2 

0 

1.2b 

0 

• 10 

3 

1 

l.bO 

0 

• 10 

4 

0 

1.7b 

0 

• 10 

s 

1 

1.00 

0 

• 1. 

o 

0 

1.2b 

0 

• 10 

7 

1 

l.bO 

0 

• ln 

a 

0 

l • 7b 

0 

• 10 

9 

1 

1.00 

0 

• 10 

9 

1 

1 

1 

1.00 

0 

• lu 

3 

J 

l.bO 

0 

• 10 

4 

0 

1.7b 

0 

• 10 

b 

0 

1.00 

0 

• 1C 

6 

1 

1.2b 

0 

• 10 

a 

1 

1.7b 

0 

• 10 

iO 

0 

1.2b 

0 

.10 

iv 

J 

4 

1 

1.7b 

0 

•10 

8 

t 

1.7b 

0 

• 10 

9 

0 

1.00 

0 

• 10 

wASCb 

«LAD= 

16  C ASF  b 

HOT  COnVi.rGcd=  C 

c 

o 

ion  Correct 

« 

2 

0 

u 

**- 

It 

u 

u 

0 

c 

x: 

<x> 

C 

u 

a> 

0) 

■ 

S JC 

4J 

62  ^ 

o a. 

62 

u 

£ 

62 

ft 

6. 

If 

c 

■ 

T3 

§ 

Jd  W. 

3 0 

«i  c 
a.  c 

<6* 

a. 

** 

£ 

e w 

si 

C ft 

* 

Q)  Of 

Hi 

£ 

0 

18 

K 62 

£ S 
« *6* 

If 

1 

u 

62  0 
e ft 

u 

If 

i 

H 25 

a. 

1 £ 

£ ~ 

z 

H < 

bJ 

Z 

» * t 4 

1 

► • • • 1 

.70 

■ * • » ( 
.34 

> • • • 

9.6b 

* • 

10 

• • • • 

9.65 

* • t 

.71 

a • 

3 

2 

.50 

-•22 

5.06 

4 

4.66 

.43 

2 

3 

1.00 

10.00 

0 

5 

0 

0 

0 

4 

0 

-10.00 

0 

5 

0 

0 

0 

5 

• 2b 

-.78 

3.15 

8 

2.69 

.24 

1 

b 

.50 

-.09 

2.57 

2 

2.57 

.50 

2 

7 

.83 

• 79 

3.13 

6 

3.12 

.87 

2 

8 

.56 

-.09 

0.85 

9 

9.55 

.50 

2 

9 

.57 

•05 

8.03 

7 

8.03 

.58 

2 

lu 

.67 

.43 

3.18 

3 

3.20 

.77 

2 

11 

.63 

• 13 

9.62 

8 

9.62 

.62 

2 

lc 

.89 

.89 

4.02 

9 

4.02 

.89 

3 

13 

.40 

-.29 

3.59 

5 

3.68 

.42 

3 

14 

.17 

-1.08 

1.40 

6 

.85 

.18 

2 

lb 

.71 

.25 

7.55 

7 

7.52 

.68 

2 

lo 

.67 

• 23 

7.03 

6 

7.03 

.67 

2 

jru*  Estimated  Standard 

Error  of  Measurement 


r 


-45- 


5 


LINPSCO 

Following  are  sample  runs  from  LINPSCO  using  graded  models  and  the 
nominal  logistic  model. 

Graded  models.  The  IPOOL  file  for  these  examples  was 


i 

x • D 

3.0 

2 « u 

4 • b 

w.O 

tm 

x • b 

2.0 

1*0 

3.U 

r, 

X • w 

' 

x • b 

1.0 

0 . t 

l.b 

o . 0 

*4 

i • j 

0.0 

-1.0 

0.0 

- 1 . s 

X • J 

-1.0 

-c  • 0 

-1.5 

-w.n 

u 

x • b 

-? . 0 

-J.  u 

-3.1) 

The  DATA  file,  including  subject  identif ication  and  item  responses,  was  as 
follows.  Note  that  in  coding  the  item  responses,  a "1"  indicated  the  "best" 
response  and  a "3"  indicated  the  "poorest,"  as  specified  by  the  item 
difficulty  parameters  in  IPOOL. 


1 

jw  J.xi  1 

2 

EoPll 1 

3 

„■>«:  It  IP 

4 

liPitl 

5 

o 

7 

J ^21 

b 

....PEcS 

9 

111 J33 

10 

l-llEi? 

JLl 

E v>  btiP 

1?. 

?j3i) 1 

PtDlil 

i.4 

Pxlill 

The  following  is  an  example  of  the  logistic  graded  model  (0PT4-1)  with  a 
1.7  scaling  factor.  In  this  example  the  b parameters  for  the  items  were  taken 
from  columns  3-4  of  the  IPOOL  file.  The  option  and  format  cards  for  this 
example  were 

o o^Ol  1 | . 7(J  4 

ill  #lX»f-i.i#2(  1X»F4.1>  ) 

<22^22 

1 ?.  3 4 *>  6 

oOujOJ 
( 2/t  7 * Al » o 1 1 ) 

uXmv.PLE  Rj,j  0'  THE  LOGISTIC  GRADED  MODEL — USES  THE  FIRST  PAIR  OF 
RA^aVE T -Rj  FRO',,  ITEM  PQwL 


B 


r 
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The  output  from  this  run  was  as  follows: 


tlNPgCO 


uINEAH  POLYCHOTOMUS  SCORING  rflTH  TWO  PARAMETER  MODELS 


st/mtu’s?  0f»s?°GM“ 

UNIVERSITY  OF  MINNESOTA 

mPls.  minn.  55455 


INUP  = 

mmax  = 

lOMir  r 
oPTl  = 
0PT2  = 
0PT4  = 
MAXCaT  = 


6 

6 

4 

0 

1 

1 

2 


U =1.7 

VARIaULE  FORMAT  FOR  POOL  = ( I l * 1X*F3« 1 . 2 < lX»Fu . 1 ) I 
VARIABLE  FORMAT  FOR  DATA  = (2a7»A1.6U> 

tXAMPLE  RUN  OF  THE  LOGISTIC  oPADEO  MODEL— -USES  THE  FIRST  PAIR  OF  B 
r-ARAmCTCHS  FROM  ITEM  POOL 


I (EM  ID  = 1 

A:  1.50 

a:  3.oo  2.oo 


b:  2.00  i.oo 


n5?  10  i.50  3 
b:  i.oo  o 


ITEM  ID  = 4 

A:  1.50 

a:  o -i.oo 


x,l?  10  r.5o  b 

b:  -i.oo  -2.oo 


HEM  ID  = 6 

a:  1.50 

a:  -2.00  -3.00 


REJECTION  = 0 


REJECTION  = 0 


REJECTION  = 0 


REJECTION  = 0 


REJECTION  = 0 


REJECTION  = 0 


CmSLS  READ  = 14  CASEb  NOT  CONVEHGEn  = 


• /*"  ‘ 


X 

.33 

-•SO 

4.72 

6 

2 

4 • pS 

.46 

ft. 

.So 

• so 

4.72 

6 

2 

4 . 25 

• 4b 

%i 

• 3 3 

05 

3.69 

6 

3 

4.28 

.62 

•* 

• 6 i 

1*62 

2.63 

6 

3 

4.11 

• o 2 

••do 

4.21 

6 

3 

4.  Hi 

.49 

J 

-•JO 

5.16 

6 

2 

4.41 

.44 

/ 

• .JwJ 

•si 

4.60 

6 

2 

4. 25 

.46 

o 

-•Si 

4.«rt 

o 

2 

4*25 

.46 

J 

. 6-J 

-•00 

.96 

6 

3 

4.41 

1.02 

1 u 

• So 

• 52 

2 . 06 

o 

2 

4.25 

.61 

1* 

0 

“i  • SI 

4.06 

6 

2 

4.11 

.45 

1a. 

.5b 

— • 00 

4 . ?0 

6 

3 

4.41 

.49 

1^ 

. r».j 

1*  Jl 

4.06 

o 

2 

4.11 

.45 

1h 

.00 

2 • o9 

3.12 

6 

2 

2 . 67 

.57 

Following  is  an  example  of  use  of  the  normal  ogive  graded  model  (OPT4-2) 
using  the  same  DATA  and  IPOOL  as  the  previous  example.  In  this  example  the 
b parameters  for  the  items  were  taken  from  columns  5 and  6 of  the  IPQOI.  file. 
Input  control  cards  for  this  example  were 


o 6201  2 4 

a I f 1 x . f 3 . 1 # 1 1 X »K4 . 1 1 lx  »F4 • 1 ) 

224222 

1 2 3 4 5 6 

b 0 b u 0 0 

12A7»A1*oI1) 

uXA’^PLl  Rj.i  OF  THE  MORTAL  OGIVE  GRADED  MODEL— USES  SECOND  PAIR  OF  B 
H Ana ME  f Ei<S  FRO*  ITEV  POOL 
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Output  was  as  follows: 


LINPSCo 


LINEAR  POLYCHOTOMuS  SCORING  WITH  TWO  PARAMETER  models 


KKHWKTO  fJf«oS?OG*AM 

UNIVERSITY  OF  MINNESOTA 

mPls.  minn.  55455 


m = 1 

iOMIT  = 4 

OPT1  = 0 

OPT2  = 1 

0PT4  = 2 

MAXQaT  = 2 

VARIaULE  FORMA'!  FOR  POOL  =(Il»lX»F3*l*llX»F4.itlXtF4,l) 

VARlAbLE  FORMAT  FOR  UATA=t2A7.Al»6ll) 

EXAMPLE  RUN  OF  The  normal  OGIVE  GRADED  MODEL--USES  SECOND  PAIR  OF  B 
parameters  From  item  pool 


HEM  ID  = 1 REJECTION  = 0 

a:  1.50 

b:  4.50  j.oo 


HEM  10  = cn  t REJECTION  = 0 

A • 1 • DU 

b:  3.00  1.50 


ITEM  ID  = „ 3 REJECTION  - 0 

A;  1.50 

B:  1.50  0 


ITEM  ID  = 4 REJECTION  = 0 

a:  1.50 

u • n _ i . c.n 
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X 

• 3 j 

-•75 

2.97 

6 

2 

3*36 

.58 

Cm 

.50 

. 75 

2.47 

6 

2 

3.36 

.58 

w 

.32 

-•04 

P.O0 

6 

2 

3.42 

.35 

. 6 / 

1.35 

0.05 

6 

2 

3.39 

.35 

'j 

.32 

-•00 

11 .98 

6 

2 

3.42 

.29 

U 

0 

-.00 

11.58 

6 

2 

3.42 

.29 

/ 

.35 

1.09 

c* . 50 

b 

2 

3.39 

.32 

o 

0 

-1.09 

°.50 

G 

2 

3.39 

.32 

J 

. 5o 

-•00 

12.80 

6 

1 

3.42 

.28 

lu 

. 5u 

.75 

12.29 

t> 

2 

3.36 

.29 

1 X 

0 

-2*25 

5.10 

b 

2 

3*35 

.44 

1C. 

. 5o 

-•00 

3.19 

U 

2 

3.42 

. 5o 

1^ 

.5i.i 

2*25 

5.11 

6 

2 

3*35 

.44 

1 + 

.33 

5*92 

2.0  5 

G 

2 

2*25 

.70 

Nomi nal  logistic  mode l . Following  is  an  example  of  use  of  the  nominal 
logistic  model  (OPT4=3).  The  DATA  file  was 


10001  9946454i34llU22442211121114llll4lllll211111Ulll 

10002  301 0w.b34b4111441 14 14 141 14411 11 14illll  14 1 1 14111414111 

10003  38  9Hn00b424l411«illl44l2l4l2lll4llll4411illl4llli241 
iOuOb  381 0424^042422244442421 114244242 12224212 121212222241 

10004  j8  94b43444l4l4444l41444444l411l444111l4ll414U44411 
IOuOb  j81  0452t>b4il  1441  jl214bbll2l  1223  li I3l 1 1 1 1 S55555555551 
I0u04  j810-b445443bll 111514111 14411 151 il 1411411 11111114111 
IOuOb  38  9^33494232455143233143524155555555555555555555555 

10009  38l0^34444i42l24x52232l344i3i242i4124221j 1432i222554 

10010  381051509444444111444511 141 11141 14141111 111111144141 

10ull  J01O329454ilil  11x1511411 13  xlillllln  1111  lllllllH  111 
10ul2  3810-b492442323n34333ll422225l3x325ll22l44511242l5l 

10013  3810540444414444x4411144144411 1441411441 114241511111 
lOulb  38  9324bb4l445l2l3254441l4illHlili21241i34411112l41 

10014  38lO5903b421l24lx42144l44441«;ll2i4i21141il4111l54411 

100 15  3010312954211141x124411 i 1411 lllixl 1 1 11 1 l i 1 1111 ll 1141 
10ol4  38l03952b5244l32ll3142111425l211llll4111i41111l3l45l 
lOoltt  3810324844114411X41324111414114^44431441114114115141 
10U19  38103282544 1441 lX442411124o3l 1344241 124 112131 11 141 14 
10020  2&10h 11434444444141441 1144444414X4444141 114144144444 
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Input  control  cards  were 


43  22300  3 5 

(I4#4<1X.F9.6> ./»4X.4(lX*F9.6> ) 

0333333-5333333^333333333333333*333333333333 
3417  3251  3422  3421  3277  3408  3010  3405  3213  3079  3002  3220  3404  3430  3021  3221 
3076  323b  3029  3078  3239  3023 
IOOOOOIOOUOOQOOUOOOOOO 
(2A7.A1.43I1) 

EXAMPLE  RUN  OF  THE  NOMINAL  LOGISTIC  MODEL 


And  the  output  from  this  run  was 


LINPSCO 


LINEmR  POLYCHOIOMuS  SCORING  WITH  TWO  PARAMETER  MODELS 


SHOWS?1 8 ^c,8M0SRAH 

UNIVERSITY  OF  MINNESOTA 

mpls.  minn.  554b5 


INUP  = 43 

M MAX  = 22 

lOMIl  = 5 

opt 1 = n 

0PT2  = 0 

0PT4  = 3 

MAXCmT  = 3 

VARIABLE  FORMAT  For  POOL  =(I4*4(1X»F9.6) #/.4*»4(lX*F9.b) ) 

VARIABLE  FORMAT  FOR  uATA=(2A7t Al»43ll) 
example  run  of  the  nominal  Logistic  model 
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1 l EM 

ID 

= „ u 

X: 

0 0 

bs 

0 0 

I l£M 

ID 

= 32ol 

a: 

1.04  .01 

b: 

.91  -.71 

"VZ 

ID 

^.bH^-.Ob 

H J 

.96  -.61 

“S? 

ID 

” .9?***"  ‘ . 46 

o: 

1.32  .20 

1 1 EM 

ID 

= 32.7 

a: 

.94  .13 

b: 

1.75  -.20 

nw 

ID 

= 3406 

1.27  .17 

b: 

.72  -.50 

11  EM 

ID 

s „ 0 

X: 

0 0 

o: 

0 0 

I 1 EM 

ID 

= 34ob 

a; 

1.22  -.16 

b: 

.70  -.47 

l'H 

10 

= 32.3 

.94  .12 

b: 

1 .62  -.23 

1 tLM 

ID 

= 30  79 

a: 

1.30  -.14 

q: 

1 • 75  — * 66 

1 JEM 

10 

= 30 

a; 

1.1b  -.01 

b: 

.16  -.bl 

‘■it? 

1J 

1.3^-. 04 

g: 

1.66  -.09 

“W 

11 

= 34  0 4 

.99  -.07 

b: 

1*02  — . 66 

hem 

ID 

= 34^u 

a; 

1.24  .08 

b; 

2 . 36  — . Oo 

ns? 

10 

= 30.- 1 

2.04  .39 

b: 

3.71  .13 

“i? 

ID 

1.31“' ‘"‘.13 

b: 

2.12  • 28 

I 1 EM 

ID 

- 30/6 

a; 

1.29  .23 

b; 

3.02  -.45 

1 i EM 

a: 

1U 

I.232j2.ib 

b: 

1.77  -.24 

*'S? 

ID 

= .6^.16 

b; 

1.31  -1.38 

1 fEM 

ID 

= 30/8 

a: 

1.69  .13 

u: 

3.63  -.59 

“sr 

ID 

r.332j\o4 

b: 

2.39  -.11 

ns? 

ID 

1.43°*"-. 07 

g: 

2.24  -.33 

CmSLS  PLAU  = 20 

* ■*  tr  v 

REJECT  ION^  1 
b 0 


REJECTION  = 
-1.1b  .li 
-1.43  1.24 


WU??4  = 

-1.59  i.ib 


-?^^CTi?S4= 

-2.29  .77 


= 0 


REJECTION  : 
-1.01  -.Ob 
-l.HO  .41 


REJECTION  — 
-1.19  -.24 

-1.4o  1.24 


NtJjjCl  10NQ=  1 

U 0 


REJECliON  - 
-1.0b  -.01 

-1.27  1.04 


-?^CI1?iV  0 

-2 . *4 1 1.01 


Rc-JlCTIOI.  = 0 

-1.20  -.04 

-2.0u  .97 


REJLCUON  = 
-l.bO  .3b 
-1.04  1.99 


. 14 


-1.73 


RloLCTION  = 0 

-1.0b  .14 

-1.1b  .79 


REJECTION  = 
-1.44  .11 

-2.64  1.13 


-3.79  -.Ob 


pljection  = 

— 1 . 1 3 — .3b 

-2.bb  .09 


Rt JeCT  iON  = 
-.69  -.02 

-2.16  -.39 


rejection  = 
-.92  -.12 
-1.(14  .31 


-?^CTI?!ir 

-1 .69  l.bb 


REJlCT ion  = 
— 1 • 3 o — • 4o 
-3.13  .10 


—3 . bo  1 • 3 1 


^^CT*?y,: 

-1.64  -.37 


CASES  NOT  CONVERGED  = 
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•H 
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0) 
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T3 
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Vi 
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« 
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4-1 
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to 

•H 

5 
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»-• 

03 

to 

o 

T3 

4-1 

•H 

0) 

6 

Vl 

•H 

C 

•H 

(4-1 

rH 

4J 

<u 

a 

4J 

to 

■U 

o 

a) 

4-1 

4J 

4J 

(0 

4-4 

c 

to 

M 

M 

E 

C/3 

0) 

c 

•H 

a. 

M 

TD 

o 

•J 

(4-4 

U-4 

o 

T3 

M 

•H 

a> 

o 

o 

(4-1 

0) 

■U 

e 

03 

c 

4-4 

<U 

Vi 

c 

u 

Vl 

M 

3 

Q) 

O 

6 

o 

Q) 

0) 

E 

U 

a 

•H 

a 

X 

4-1 

•H 

05 

o 

X 

03 

g 

B 

03 

4J 

CD 

v» 

(0 

a> 

3 

O 

0) 

03 

H 

a. 

X 

Z 

z 

H 

W 

6 99404 

5 

.50 

-*63 

9.16 

20 

3 

9.12 

6103034 

5 

.65 

-.02 

6.24 

20 

3 

6*21 

6 9hQuO 

5 

.55 

-•33 

7.94 

20 

2 

7.60 

6104245 

0 

.15 

-1.31 

12.19 

20 

2 

12.11 

8 99543 

4 

.35 

-1*10 

11.43 

20 

2 

11.44 

8104526 

5 

.59 

-.47 

7.11 

17 

2 

7.11 

8105544 

5 

.71 

-.07 

5.17 

17 

3 

5*16 

8 99334 

9 

.19 

-1*73 

9.55 

16 

2 

9.55 

8105344 

4 

.32 

“1*40 

11.69 

19 

2 

11.89 

6105130 

9 

.47 

-.69 

9 . C 2 

19 

3 

9.00 

8105294 

5 

♦ 93 

3.12 

.59 

19 

3 

• 59 

8103849 

2 

.lb 

-1.90 

10.66 

19 

3 

10.85 

8105404 

4 

.45 

-.72 

° • 60 

20 

3 

9.59 

8 95245 

5 

.56 

-*5l 

7.98 

18 

3 

7.97 

B1059b3 

5 

.40 

-1  • 05 

11.34 

20 

2 

11.24 

8105129 

5 

.70 

• 00 

6.17 

20 

3 

6*14 

8105952 

6 

.47 

-.85 

10.09 

19 

3 

10.09 

6105248 

4 

.55 

-.49 

8.39 

20 

3 

8.39 

8105282 

5 

.40 

-1.03 

11 .18 

20 

3 

11.18 

5109114 

3 

.3u 

-•98 

10.94 

20 

3 

10.95 

.33 

.40 

.35 

.29 

.30 

.37 

.44 

.32 

.29 

.33 

1.30 

.30 


APPENDIX  C 

LINDSCO  FORTRAN  PROGRAM  LISTING 


m 


PROGRAM  LINDSCO  (INPUT, OUT  PUT,  OAT  A, IPOOL, T APE1=DA T ft , T AOf 2 = I aOOL , T A 
1PE  3, PUNCH! 

0 1  HE  NS  ION  ITEH(600I,  *1600),  BT600T,  CI600I,  INAO(JJC>,  KFY(300>, 

1 1 °E J(  300  I * IFORN(6l , IRAN(300),  INADS(300>,  IRESP(3jO>,  ISAOIFO.25 
2»,  AON (300),  BOH ( 3001,  COHI300I,  ISAOS(300»,  0ESr(?4»,  NAMF(L),  I F 
30»H1(8>,  NISS(25,2» 

I N 1 1 Gt R 0PT1,0PT2,0PT  3,  OPT  4 

PtAL  ITOT 

N=NC=0 


RcAD  OPTIONS  ANO  PROGRAM  PARAMETER  FROM  INPUT  FILE  DATA  IS  TAP£2 


I POOL  = 2 

READ  SO,  lNUP,M,OPT l,OPT2,OPT3,OPT4, TS,TSS, AN AX , AM 1 N, UN A X, CM AX , I OH 
1IT 

pEAO  S3,  (IFORNtm  , 1=1,91 

p E AD  (IPOOL, IFORM1T  ( I T £ H( 1 1 , A ( I » , a ( I) , C ( IT , 1=1 , I NU“I 


•START  c.  AU1NG  THE  SPECIFIC  OATA  (SPECIFIC  FO°  THE  PUNT  FifjH  TH-  pPUT 

• p.ao  is  the  item  io#s  administered 

•IP-.  J IS  THr  PtJECTEO  ITEM  ID  S 

•Kt*  IS  THE  KfcV  FOP  TH-  ITEMS  IN  INAO 

•NS-.  C is  TH.  NUMBER  OF  SJPSCALES  THAT  MILL  BE  GIVEN  TO  THE  -'JOG-AM 
•MISS  will  MAVt  THE  *UMPEP  OF  I T c MS  IN  EACH  SUBSCALE  TOGEIM-->  WITH  (H- 
*,-<t  OF  THE  SUbSCALES 

•IS  10  IS  THE  ITrM  ID  S IN  EAC“  SUBSCALE 


!- AO  51,  CINAD(IT,I=1,M» 

Jt  40  52,  ( IRE  J ( I ) , 1 = 1 ,M> 

P AD  52,  (REVm,I=l,N| 

Rc  AO  53,  (IFOPMdl  , 1 = 1,  3| 

p£AO  55,  OESC 
PRINT  49 

PRINT  54,  INUP,  M,  IOMI!  , 0PT1,0PT2,0PT  J.  3PT4,  TS,TSE  , A-<a*,  THAX  , phIn,  C 
1HAX, IE  0RM1, IFCRM,UESC 

C THOSE  ITEM  10  S THAT  APt  IN  IP£J  ARfc  SET  TO  7c  PO  2t-?0  ITE“  10  S 

C WILL  BE  SKIPPED  DURING  THE  COMPUTATIONS 

CO  1 1=1, M 

1 IF  (1»EJ(II.EQ.1>  INADT  I)  = D 

C HEAD  SUBSCALES 

READ  48,  NSSC 
IF  T‘SSC.EQ.0)  GO  TO  5 

-£AO  5b,  (NISS ( I,:l ,NISS (I ,21 , I=1,NSSC» 

C -EAD  SSC  INOEX 

JO  2 Jj=l,60 

!)0  2 11  = 1,25 

2 ISADt JJ,I1)=0 
30  3 1=1, NSSC 
NI=N1SS(I,1I 

-'EAD  51,  (ISAD(J,IT  ,J=1,NIT 

3 CONTINUE 

00  4 11=1, NSSC 
NI=NISS( II, IT 

PRINT  58,  II,(ISAD(JJ,II),JJ=1,NIT 

4 JONTINUE 
POINT  49 

5 JONTINUE 

IF  (OPT2.E0.1)  GO  TO  9 


1 

? 

T 

4 

5 
b 

7 

8 
9 

10 

11 

12 

IT 

14 

1 5 
IF 

17 

18 
IP 

23 
21 
2? 

2 ’ 

24 
2 ~ 
?5 
?7 
23 
2B 
3' 
31 
3? 

33 

34 

35 

, , 
3* 
3° 
4 1 

41 

42 

43 

44 

45 

45 

4 T 
43 

40 

51 

51 

52 

53 

54 

55 
55 
57 

5 * 
c 9 
63 
61 
e? 

63 
f 4 
65 
6F 


* * 57 

•IN  THl  N-  x I no  LOO'-'  THE  TTEM  PARAMETERS  CORRESPONDING  TO  THE  TTFHS  * 68 

•IN  luao  ARt  PMfIVH)  FnOH  A.n.C,  AND  LOADER  INTO  AOM,  BOH,  COM  Pf$B.  * 60 

•THE  cMKlt  > IN  THE  AOM, RDM, COM  API  ZfROEO  TON  1 »<■  CASE  OF  ZERO  ITFM  * 70 

•It)  IH  TM-  INAC  * 71 

* * 7? 

* *73 

DO  8 J = 1,M  75 

IF  UNAnjJ)  .NT  .0)  00  TO  6 75 

ADm(  J|=>«UM(J)=C0M(  J)  =0  77 

GO  Tr  8 70 

6 C 0 N T J MJr  79 

IFOUNQ=0  AO 

DO  7 i=l,INUP  81 

IF  (ilvAOIJI.N.-  .IlrKlII  GO  TO  7 AT 

ifo’)nd=i  at 

AnAIJI=A(II  A A 

•TO -A  ( J»  =0  C 11  85 

CD(1TJI=CCT»  A 6 

GO  111  A 8 7 

7 CONTINUE  8 8 

IF  I I'-OUND.tC.O)  IfIA(J(J»=0  89 

8 CONTINUF  90 

GO  Tr;  u qi 

C THt  NEXT  Si  CT  10*1  IS  I'SFO  it  OPTION  i IS  ON,  IT  MILL  t AKF  THE  FT  9? 

C f»AHAMHEFS  T-7CN  THt-  POOL  HIT  MOIIT  MAKING  lISc  OF  INAO  98 

9 00  l"  1=1, M 94 

AOM(i)=A(I)  95 

HDM(I»=Rm  95 

CDM(I)=C(T>  97 

10  CONliNUt  98 

0 IT  001  I ON  3 IS  ON  TH*-  PARAMETERS  A.B.C  ARE  CONSTRAINED  HltHIN  9 99 

C OF  AMA*. ATTN, UMAX, OMIN.CHA*  100 

11  IF  (ORTJ.F  0.0  GO  TO  13  101 

00  lc  1=1, M 10? 

IF  ( 1 N AD  ( I ) , t 0 , 0 ) ',0  TO  1?  103 

IF  ( ( ADMTII ,GT. AMAX|  ,09.  (OPT3.FO.3H  A 0 M ( II = A MA  X 104 

IF  ( “DM  TIT  . L T • MM  IN  I HDMU>=9MIN  105 

IF  (rfPMTII .GT.BMAX)  BOM ( I ) =BMAX  105 

IF  KC08lll.fiT.CMAI)  .OH.  (OPT3.EQ.2H  C()M(1)=CMAX  1C7 

1?  CONTINUE  1 C 8 

13  IF  (nPTl.NF.lt  GO  TO  IE  109 

DO  14  1=1. M 110 

PUNCH  57,  INADT  I » ,<FY  < 1 1 .IPf  J«  It  , AOMTI  » ,H()M(T»  ,CD**CI>  111 

14  CONIlNUr  11? 

15  PRIM  57,  (IMOm  .KEY  ( It,  IPtJdl  ,AOM<II  ,90M<I>  ,COM<I>  ,I  = 1,N>  113 

PPINI  49  114 

* * 116 

♦ * U7 

•RE  AO  A SUBJECT  FPOM  TAP:!  CALCULATE  THt-  T A , LOOP  BACK  TO  5 CTC.  * 118 

* * 119 

• * 120 

16  HE  AO  (1.1FORHI  NAME, ID, (IRAHtll , I = 1,M»  1 22 

00  17  J J= 1 , M 123 

17  I Rt  SP ( J J ) = 3 124 

IO=UNIT (It  125 

C CHECK  THt  END  OF  FILT  ON  DATA  FILE  125 

IF  (10. LE  .31  GO  TO  18  127 

PRINI  59,  10  128 

18  IF  (IO.EQ.O)  (-0  TO  97  129 

N = N+ 1 130 

1TOT=1.0  131 

C SET  THE  ITFM  in  jo  Z-PO  FQP  THE  0MMISS10NS  (THAT  IS  READ  IN  FPO  13? 
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C 133 

C ANO  Ski  THE  RESPONSE  VECTOR  TO  1 IF  THE  ANSMER  IS  CORRECT  134 

c 13*; 

DO  21  1 = 1, M 136 

1 N AOS (11= INAO (II  137 

IF  (IPAHin.EO.  IOMIT)  19,20  13* 

19  I NADS ( I ) =0  139 

GO  TO  21  140 

20  IF  <<IRAN(II .EQ.KFYtlM .AND. (IRE J(I) .EQ.0>>  IRESP(I>=1  141 

ITOT=ITOT*FLOAI (IRtSP(I) ) 14? 

21  CONTINUE  143 

I T I NAOS=  0 144 

00  22  KK=l,M  14t' 

22  IF  (INADS(KK)  . fc  0 . 0 1 ITX’jAOS=IT  I NADS* 1 146 

KL=M- I T I NAOS  147 

ITOl =ITOT /FLOAT  I CKl ) ) 14* 

IF  (1TOT.EQ.O.T  23,24  149 

23  PRINI  60,  NAME, ID  1*0 

IF  (0PT4.t3.3J  GO  TO  26  151 

T=-1G.0  15’ 

SFOR«=0.0  153 

TINFO=0.0  154 

I KL  =*L  155 

I T Ec=  0 156 

iXPTOT*Q.O  157 

HRITt  (3,63)  NAME  , ID,  HOT,  T, ST  ORM, IKL, T INFO , E XPTOT , ITt R 15* 

GO  TO  32  159 

26  IF  (1101.60.1.)  25,26  162 

25  PRIM  61,  NAME, ID  161 

IF  (0PT4.Ed.3l  GO  TO  26  l*’ 

r = l J.O  16  3 

SFOPf'  = 0.0  1E6 

t I NFn=  o . 0 165 

IKl=*L  166 

ITEP=C  167 

E*PTOI=0.0  16* 

HR  I It  (3,63)  NAMt.IO, I10T,1,SF0RM,IKL,TINF0,E*PT0T,ITFR  169 

GO  TC  32  17 r 

171 

• * 17? 

• • 173 

• NON  I HE  DATA  IS  READY  TO  MAKE  THE  CALLS  TO  IMF  APPROPRIATE  ROUTINE  174 

*rn  EST  IMA  I t THE  THETA. OPTION  4 HILL  OETERMIN-'  THE  METHOD  BY  WHICH  « 17* 

•IMF  THETA  ESTIMATF  HILL  BE  F OUNO  * 176 

• » 177 

• * 17* 

m 

2*  I KL  =*L  IS* 

IF  (CPT4-?)  27,28,29  1*1 

27  CALL  MAXLNO  ( I PESP  , I N A3  S , M , H , ADM  , <30M,C  OM  , 6 3 , . T05,  T , SF  ORM  , I F A IL  , I I N 1*2 

IF  0«t  XPTOT, ITEK.SEM)  1*3 

GO  TO  30  1*4 

2*  CALL  MAXLK  ( M,  INADS  , I»E  SP,  M , AEjM,  BOM  , CDM,  5 0 , . 0 05  , I F A IL  , SF  ORM,  T,  T I NF  185 

1C. tXPTOT,  ITER, SEMI  1*4 

GO  TO  30  1*7 

29  I = T S 18* 

SF0RM=ISS  139 

SfM=L.O  19C 

CALL  MAYfcS  (M,INAOS,lKE5P,M,A9M,HOM,COM,T,SFOPM.TINFO,FXPTOT,ITtRI  191 

GO  TO  31  19? 

33  CONTINUE  193 

IF  ( IF  AIL  .EQ. C)  60  TO  31  194 

PRI*II  62,  NAME,  ID  19c 

C SFORM  ANO  T ARE  Si T TO  -99.99  IN  MAXLK  IF  NOT  CONV  196 

NC=NC*1  197 


**»K  I TM!  SUCCFSFUl  RtSILTS  TOTHF  FILE  TAPE  3 

* 


* 


31  WFIP  ( 3 , c 3 T f.AMt,ID,IloT,T,SFORN,IKL,TTNFO,cXPIOT,IT'R,SEM 

Z IF  THEM  I?  NO  SU8SCALE  CALCULATIONS  LOO°  °ACK  TO  PtAO  A SU8JEC 

C CIHEPWUF  CONTINUE  NT  TH  THr  CP**PUIATTONS  FOR  SU3SCALES 

32  If  (NSSC.'-O.Ol  GO  TO  16 
10  uf  I=1,NSSC 

N I = ’)I  SS  ( 1 , 1) 

I10T=u.P 
'JO  3a  J=1,N 
ISA  os ( J I =INADS( JT 
CO  34  <=1,NI 

IF  « lSA1(<,I)-INAnS( J»)  34,33,34 

33  ITOT-ITOI+FLOAT (IV-SP(J) I 

GO  TO  35 

34  CONTINUE 
ISADST J)=3 

35  CONTINUE 

I T IFAOS=0 

no  ?«  KK  = i,M 

36  IF  ( ISAOS(XK) ,Nt • » ) IT1 SAOS= 1 1 IS AOS* 1 
KL  = I I ISADS 

ITOTMTOI  /FLOAT  ( KL  > 

IF  ( I TOT .f  0. G . OT  37,  3H 
3 7 PRIM  66,  NAME  , IO,NISS(T  ,21 

IF  (PPTh.EQ.3T  GO  TO  40 

r = - 1 r . o o 

T INE('  = EXPT')T  = C.  0 
I T E p = u 
SFCHM  = 0 • C 

WRIV  (3,o4T  f AMF,in,NlSS(I,2l,IT0T,T,SF0RM,IKL,TINF0,EXPT0T,ITER 

GO  Tfl  4b 

3 H IF  ( 1 TOT  .5  0.1 .01  40 

30  PPlNI  67,  NAME,I0,N1SMI,2) 

IF  (pPT4.EQ.3T  GO  to  <3 
T = 1 :.c 

T INFI'  = EXFT0T  = C.  1 
ITEP-u 

SFO°M=0 . 0 

WRITE  (3,6-4  1 NAME  , 10,  MSS  ( 1,21  , I TOT,  T.SFOWM.IKL,  T INFO,  EXPTOT,  ITEI? 
GO  TO  46 

40  I <L-<1 

IF  (OPT  4 — ? I 41,42,43 

41  CALL  P AX  L NO  ( I PE  SP.  I S AOS  , Ft,  M,  ADM  , BOM,  CDM,  6 3 , . 01 , T , S FORM,  IF  A I L , T I NF 
lO.tXPTOT,  IT6  P, SEMT 

GO  TO  44 

42  CALL  MAXLK  ( M, I S A )S , IPE SP, M, ADM, BOM, COM, 5 0 • . 0 1, I F A I L, SFORM • T , T I NFO 
l,EX»TOT,lltR,StN) 

GO  TO  44 

43  I = TT 
SF06«=TSS 
SE  M = U . 0 

CALL  MAYtS  (M, ISADS, IkESP.M, AON, BOM, CON, T.SFOPM.T INFO. EXPTOT,ITFRT 

GO  TP  45 

44  IF  (1FAIL.EQ.0T  G-1  TO  46 

PRINI  66,  NAMF,IC,NISS(I,2) 

45  WRIIr  ( 3 , 64T  NAME , 1 0, NISS ( 1 , 21 . ITOT , T, SFO RM, IKL , T INFO, EXPTOT , ITER, 

1SEM 

46  CONTINUE 
GO  TO  16 
PFINI  65,  N, NC 


47 
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C 265 

bn  FORMAT  <151  266 

49  FORMAT  <1H1)  267 

50  FORMAT  < 214, X, All, 6F 5.2, 2X, 12)  268 

51  FORMAT  ( 1615)  269 

52  FORMAT  (6011)  270 

5 J FORMAT  ( 8 A10 ) 271 

54  FORMAT  <T50,*LINDSCO*,/,T50*==««=**. //////, T20, ’LINEAR  OICHOTOMUS  272 

1 SCORING  WITH  THREE  PARAMETER  MODELS*, ////, T4P*PSVCHOMETRIC  METHOD  273 

2S  PROGRAM*, /,T40*DEPARTMENT  OF  PSYCHOLOGY*, /, T40*UNIVERSITY  OF  MIN  274 

3NESOI A*,/,T40*MPLS.  MINN.  55455*, /,///, T20* INUP* , T27*=* 15,/ ,T2Q,* M 275 

4MAX*, T27*=*I5,/,T20,*IOMIT*,T27*=*,I5,/,T20*OPTl*,T27*=*I5,/,T20*O  276 

5PT2*,T27*  = *I5,/,T20*OPT  3*, T 27*=* 15, / , T 2 0* 0PT4*, T2 7*=* , 15 , /, T20*TS*  277 

6»T27*=*F5.2»/,T20*TSS*,T27*=*f'5«2»/»T20*AMAX*»T27*=*F5.2»/»T20,*DM  278 

7 AX*. T 27, *=*F5.2,/,T2  0*BMIN*,T27,*  = *F5.2,/,T20*CMAX*,T27,*=*F5.2,/,  27  9 

8 T 2 0 » * VAR  I ABLE  FORMAT  FOR  POOL  = *8A10,/,T20, ’VARIABLE  FORMAT  FOR  OAT  28C 

9A=*8A10,/,T20,8A10,/,T20,8A10,/,T20,8A10, /!  281 

55  FORMA!  ( 8A 10 ) 282 

5b  FORMAT  <215!  283 

57  FORMAT  < *1*,29X , • I TE  M ID  S*2 X, *KEYS*2X , ’REJECTIONS* ,4X ,* A*, 9X,*B* , 284 

19X,*C*,/,63</30X,I5,5X,I3,5X,I3,lX,F10.2,lX,F10.2,lX,F10.2!)  285 

58  FORMAT  < ///,40X, ’ITEMS  IN  SUBSCALE  NO= • 13 ,/, 1 0 < 20 16 , /) 1 286 

59  FORMAT  < 10X, ’PARITY  ERROR  ON  T APE  * , 9 OX , 1 2 I 287 

60  FORMAT  < 10X, ’SUBJECT  =*,2A10,*  ID  =*,A9,*  MAS  NO  ANSWERS  *,*CORPF  288 

1C T IN  TOTAL  SCALt*!  289 

61  FORMAT  <10X, ’SUBJECT  =*,2A10,*  IO  =*,A9,*  HAS  ALL  ANSWERS  *,*CORPE  290 

1CT  I*'  TOTAL  SCALt*!  291 

6?  FORMAT  <1QX, ’COMPUTATIONAL  PROBLEMS  WITH  SUBJECT  =*,2A10,*  IO  =*,A  29? 

19,*  IN  TOTAL  TE*)  293 

63  FORMAT  <X,2A10,A9,*  T * , F5 . 2 , 2F 7 . 2 , 1 4, 2F 7 . 2, 14 ,F7 . 2)  294 

64  FORMA!  )X,2A10,A9,I2,F5.2,2F7.2,I4,2F7.2,I4,F7.2I  295 

65  FORMAT  </10X,*CASES  REA0=*,I6,*  CAStS  NOT  CON VE RGED  = * , 15 ) 29F 

66  FORMAT  < 1 0 X , *SUB Jt C T =*,2A10,*  IO  =*,A9,*  HAS  NO  ANSWERS  *,*CORRF  297 

1CT  IN  SUBSCALF  *,I5)  298 

67  FORMAT  < 1 0 X , *SUB Jt C T =*,2A10,*  IO  =*,A9,*  HAS  ALL  ANSWERS  *,*CORRF  299 

1CT  IN  SUHSCALE  *,T6)  308 

68  FORMAT  )10X, ’MAXIMUM  LIXELIHOOO  ESTIMATION  DOES  NOT  CONVERGE*, ’FOR  3C1 

1 THE  SUBJtCT  = *,2A10,*  IO  = *,A9,*  ON  SUBSCALE  *,I5)  302 

END  303 

SUBROUTINE  BAYES  < M , I ! M , RE SP , N , A , B , C , B THE T , B V AR , T INFO , £ X PT OT , ITE R » 1 

INTFGER  PESP(M) , ITHTM!  2 

REAL  A (NT ,0<N) ,CCN)  3 

DO  1 1=1, M 4 

IF  CITMTII .EQ.0!  GO  TO  1 5 

CALL  BSCOR  < BTHt T , 3 VAR, B ( I) , A < I ! , C < I ) , RE SP ( I ) I 6 

1 CONTINUE  7 

CALL  NOSTAT  < M, IT M, A , H, C ,B THE T, T INFO,E XPT OT T 8 

ITER=0  9 

RETURN  in 

ENO  11 

SUBROUTINE  BSCOR  < 9 THET , BV AP , PIF , OIS.GUF S P, IRFSP)  1 

0= (OIF -B THE T ) /SORT <2. 0*  <1,0/OIS**2*BVAR) T 2 

£ RF  0=ERF  NP  <0 ) ‘ T 

E0S0=EXP(0**2I  4 

IF  (EOSQ.FQ.O.OI  RETURN  6 

E OSOI= 1 « O/EOSQ  6 

XKINV=0.6*<1.0-ERF3I  7 

XL INV=GUtS°*  <1.0-GJESPT*XKINV  8 

IF  (IXLINV.FQ.O.O!.OR.(XKINV.EQ.O.O)T  RETURN  9 

XL=1.0/XLINV  10 

IF  (IKESP.NE.IT  GO  TO  1 11 

S= 0.398942* <SQRT < BVAR) 7 SOR T (1.0* < 1 . 0 70 IS**2I /BV AR ) ! * < 1 . 0/XX INV! *E 0 12 

1SOI  13 

T =1.0-1. 772454*O*EOSQ*<1.0-ERFDI  14 

HTHtl=BTMFTF<1.0-GJESPI*XKINV*XL*S  15 

9VAC  = PVAR-(1 , 0 -GUt S PT  *XRINV*xl*?**2*<T-GUFSP*XL)  . 16 


; 


1 


C*  * * 

c*** 

C**» 

c**» 

C 


c*** 

1 

c 

2 
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y E T •J*-’  N 

MTHl  I =HTHE  T-0. 797  3A5*  MVAP/SfJPT  (i.a/0IS**2F9VAR))  *EOSQI*  (1.0/d.C* 
1 1 °F D) I 

PARTl  = 1.128379/(l.CMl.3/niS**2)*(l.  )/9VAP)  ) 

PART2  = 1.0/(CD''O*(i.')*Fpcc)  ) * * ’ 

PART3  = 0.66i»19e*O*: 9SCI*  ( ! . G *fc‘ RF  0) 

BVARsRVAKMl.OPART  1*PA“'T2*PART3I 

RETUFN 

ENJ 

PE  AL  EUNCI ION  EREN3  (X) 

OAIA  A 1 / 0 . 299 i 3 0 / 

DATA  A2/-0. 284597/ 

RATA  A3/1. 421414/ 

RAM  A4/-1  ..,52152/ 

0414  55/1.051405/ 

DATA  p/0.  <27541/ 

FPEN‘>=0.  0 

IF  (X.EQ.O.G)  PE  T'l.'N 
t S=31GN(1,  3 » X ) 
r = a 

IF  ( Y .LT  .5.0)  GO  TO  1 

F PFNF=F9 

pc  Tl)PN 

Y2=Y*Y 

T=l. J/(1.0*P*Y) 

AT=MAl*IA2PCA3*(A4  + A5*U*l)*T»»n*n 
tAr=Ai/r*F(Y2) 

EPFNP=(1.0-EAT)*fc5 

RF  MJFN 

END 

SIJMPOUTIN'  PAXLK  < M , I T M , RF  SP  , N , A , 9 ,C  , M 4X  , F “S.  IF  A I L , SOP  V , THE  T A, T INF 
10, E XPIOT  , KUM I T b » SE  M ) 

EXT'  RNAL  FOULCG, SD910G 
INTtGER  KESPTM) 

OIM'NSION  AIN) . BIN),  CIN),  ITMIM) 

USES  4AX1HUF  LIKELIHOOD  LOGISTIC  SCOPING  ALGORITHM  AND  RESPONSE 
MOO~l 

HIScCTION  IS  USED  TO  PP^vIOF  THf  INITIAL  GUESS  FOP  THE 
NE  W T ON-0 AFHSON  ML  T HOD 

CALL  BISECT  (FD0L03»PfcSu,A,n,C»M,ITM,5»GllESS) 

CALL  NEHTPAP  < F DDL OG , SOOLOG, «r S°,A,0,C,M, ITM, MAX , F PS, HUM IT S, GUESS , 
ITHe  TA.SOFV.IFAIL) 

IF  (IFAIL.eO.il  1.? 

NEWTON  PAPHSON  010  NOT  CONVERGE 

CALL  NWTEPP  (THETA, SDKV,SFM,TINFO,EXPTOT) 

KF  TURN 

CALL  LGSIAT  ( M,  I T M, A , H , 0 , T HF 1 A , T T NFO, t XPT OT ) 

SEM= 1.0 /SORT (A9S (SO° V) I 

PE  TURN 

ENU 

F IT NC • ION  FODLOG  ( »•'  S P,  I T M,  A , 3,  C,  M,  THE  T A> 

INTEGER  RESP(R) , FIGHT 
0 1 ML  NS  ION  AIM),  H ( M ) , C(M),  ITM(M) 

OATA  XHAX,XMIN/200. U.-20U. 0/ 

0 A I A 0 , P 1GH T/1,7,1/ 

CALCULATES  FIRST  DERIVATIVE  OF  LOG-L IKEL I HOOO  FUNCTION  OF  A 
RES °ONSt  VECTOR  FOR  THL  LOGISTIC  MOOEL 
SIIM=C  . 0 
00  1 1 = 1, M 

IF  UTM(I).EO.O)  GO  TO  1 
X = I)*A(I)*(THETA-H(T)  ) 

IF  (X.LT.XMIN)  X=XMIN 
IF  (v.GT.XMAX)  X=XMAX 
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24 

25 
1 
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x 

4 

5 
5 
7 
ft 
9 

m 

n 

12 

13 

14 

15 

16 
17 
IF 

19 

20 
21 

1 

? 

x 

4 

5 

6 
7 
A 
9 

in 

n 

12 

13 

14 

15 

16 
17 
tft 
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2n 

21 

22 

23 

1 
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,3 

4 
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7 
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1 


c*»» 

C**» 


1 


r * • * 
• * • 
» » * 

. , • 


C»** 


~ , , , 

1 

c*** 

2 

3 


EXF=EXPTXT  14 

ae=a<i»*exf  1? 

SUM=SUM-AE/CEXF»1.0)  15 
IF  (RESP(I).NE.RIbMTI  CO  TO  1 1? 
CE  = C<n»fcXF  14 
SUM=SUM»At/CE  19 
CONTINUE  23 
FDOLOG=-1.2*SUM  21 
RETURN  22 
ENO  23 
FUNCTION  SOOLOG  (RESP.ITH, A,B,C,H,THETA>  1 
INTEGtR  RESPIH) .RIGHT  ? 
DIMENSION  ITH(H).  A I HI  » BIN).  C(H)  3 
OATA  XMAX.XHIN/200. 0,-200.0/  6 
DATA  O, RIGHT/1. 7,1/  5 
CALCULATES  SECOND  DERIVATIVE  OF  LOG-LIKFL IHOOO  FUNCTION  F 
OF  A RESPONSE  VECTOR  FOR  THE  LOGISTIC  MODEL  T 
SUM=U.O  . 4 
DO  1 1=1, M 9 
IF  (ITM(I).EQ.OT  GO  TO  1 13 
X = 0*A(IT*ITHETA-B(IM  • 11 
IF  (X.LT.XNIN)  X = XM I N 12 
IF  (X.GT.VMAXT  X = X M A X IT 
tXF=EXPCXI  14 
AE=A<I)*FXF  lr 
SUM=Sl'M-A  TIT*A£ /(  II.  0*tXF)  Ml.OFf  XF)  » 

IF  (FtSPm.NE.  RIGHT)  GO  TO  1 17 
CE=0<IT*EXF  14 
SUM=SUM»A(I»*CTI)*AE/(CE*CfcT  19 
CONTINUE  ir 
300LCG=-£.49*SUH  21 
RETURN  22 


t NO 

SUBROUTINE  ‘IISECT  (F  l.RfSP.A,  1,C,M,  ITH,  NITER,  BN  ID) 
INT-GFP  NfcSPTH) 

0 I Mt  NS  ION  A ( M ) , B(MI,  C(M»,  ITN(MT 

calculates  apppoximate  root  of  fi  my  bisection; 

H I Sf  C I I NG  NlTtR  INUMBE  R OF  IT-PATIONST  TIMF3. 

HMin  is  bt jT  ciirklnt  guess  at  »odt  theta 


INITIALIZE  LEFT  BOUND  ANO  FI (BOUND)  ANO  RIGHT  BOUNO  F 1 (90UN0T 
BL=-S. 0 
BR=5.0 
BM I 0- Q • 0 


* 

0 

11 


1 L = F 1 (RES°,  ITH,A,H,C,H,9LT 

12 

TP  = Fi(RESc,ITM,A,9,C,H,riR) 

13 

TEST  FOP  NO  ROOT  IN  INT E RVAL --RE TURN 

IF  NO  SOLUTION 

14 

IF  ( < IL* IP) .GT . 0.31  RETURN 

15 

1 F 

NOR  CALCULATE  BISECTIONS  NITfR  TIRES 

17 

00  3 1 = 1 ,NITEP 

14 

TMIP=F1(PESP,ITH,A,B,C,H,BMI0T 

19 

IF  < ( TMIO’TL) .GT.Q.O)  GO  TO  1 

2T 

REPLACE  RIGHT  80UN0  WITH  OHIO 

£1 

HR  = -9MID 

£ ? 

GO  TO  2 

23 

REPLACE  LEFT  HOUND  WITH  BM ID 

24 

Tl =thio 

25 

BL  =PM ID 

2F 

F INO  NEH  MIDPOINT  3HID 

27 

BM I D= ( BL  APR 1 /2. 0 

24 

CONTINUE 

2° 

KE  TURN 

30 

ENO 

; « 

SUBROUTINE  NENTRAO  (F1,F2,RESP,A,B,C, 

M,  IT M,  NITER, FPS, NUM ITS, GUESS , 

i 

1 ME  Ti. SORV. IFA  IE  ) 

9 

ijt  v ' • --*•  ^ 
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* 


C*** 

c»*» 


c**« 

1 


Q * * * 

c*»* 

c 

C**' 


c*** 

c*** 


intent r resp(mi 

OIM.'SION  A(M),  B<«),  CTM),  ITM(M) 

CALCULATES  T HE  ROOT  06  FI  GTV'  N ITS  FIRST  DERIVATIVE  F2 
ANO  AN  INITIAL  GU- S3  USING  NE  h T ON-RA  PHSON  METHOD 

ihc.ta  is  Apot.  to  thl  rdct:  s:irv  is  F2 ( the  t a i 

NUMI I S=0 
I HE  I A = 6llt  SS 

LOU”  UNTIL  ER’XERS  OR  NUMBC1*  OF  ITERATIONS  flCCOME  S TOO  LARGF 

FO“V  = f KRFSP,  IT  N,  A,B,0,*",TFir  TAJ 

S0RV=F2(RESP,ITM,1,9,C,>*,TMFTa) 

t Ph  = F URV/SDF  V 

I H: I a=THETA-EKP 

NUlTI  3 = NUM IT  S*  1 

tXIT  LOOP  CR  IT  E.  Si  If’N 

IF  ( (NUM1TS.LT.M1 -P|  .AGO.  (ABS(LRF’)  .GT.EPSII  GO  TO  1 
LNU  LOOP.  TFS I F 0 J CON  VI  °GF  NO  AND  Sc  T TFAIL 
IFATL=0 

IF  (AOS  (fcR-O  .LT  .L*‘S)  * L T UF'  N 

NEWTON  RAPHSOf.  METHOD  UOt.S  N01  CONVERGE 
I F A t L =1 
RE  I URN 
E ND 

SUBROUTINE  NHTERK  ( T He  1 A ,SF  ORN,  Sf  M , T INF  0,  c XPT  DT) 

SEIS  eRROR  VALUES  FOR  THE  CASF  IN  WHICH  NEKTON  PAPHSON  FAILS 
TO  CONVERGE 
I Ht  TA=-99. 39 
SF  0P“=-99. BP 
SE  M:=-99.99 

t info =-99. 90 

EXPT0T=-99.99 

RE  TORN 

END 

SUBROUTINE  LGSTA1  ( M.  II  M,A,B»C»  THETA.!  INFO, "XPT OT  ) 

0 I Hc  MS  ION  ATM),  B(M),  C(M),  IT  HT  HI 
OATA  XMAX.XMIN/12. 0,-12. 0/ 

T INFO=0. 0 
E X P T 0 T = 0 .1) 

K OUST  = 0 

00  1 1=1, M 

IF  (ITM(I).EQ.C)  NO  TO  1 
K0UN1  =KOUNT*  1 

ARG'J=-1.  7*A(I)*fTMcTA-B(I)) 

IF  (ARGO.GT.XMAX)  ARGU=XMAX 
IF  (ARGU.LT.XMIN)  AF'GO=XMIN 
PaCUIMl.O-Cm  )*(l.u/(1.3*EXB(APGU)>  ) 

0=1.  C-P 

t ARG=tXP(-APGU) 

powlwtsEARG/Kl.UFcARG)  • (l.OFcAPG)  ) 

PPHlMt=PPRIME*(l.*i-C  (II  l*A(I)*1.7 
T INFO  = TINt-0*(PPPIM-  *PPRTME  ) / (p*0) 
t X PTOf =F XP  TOT ♦ p 
CONTINUE 

FX PTC T=EXF TOT/FLOAT (KOUNT) 

«F  TURN 
END 

SUBROUTINE  MAXLNO  ( Rf  SH  , I T M , M , N,  A,  B, C . MAX  , E PS,  THE T A , SORV, IF A IL, T I N 
1F0,IXPT0T,NUMITS,SiM) 
tXTFRNAL  F ONOGV, SONOG V 
INTEGcR  RESPTM) 

OIMENSION  ITM(N),  A ( N ) , BTNI,  C(NI 

OStS  MAXIMUM  LIKELIHOOD  NORMAL  OGIVF  SCORING  ALGORITHM  ANO 
RESPONSE  VECTOR 

BISECTION  IS  USED  TO  PROVIDE  THE  INITIAL  GUESS  FOR  THE 
NEKTON  RAPHSON  METHOO 


3 

4 

5 

A 

7 

9 

9 

10 

11 

12 

13 

14 

15 

16 
IT 
IB 

19 

20 
21 
2? 

23 

24 

25 
1 
2 

3 

4 

5 

6 
7 
B 
9 

IP 

1 


4 

B 

6 

7 

A 

9 

IP 

11 

12 

13 

14 

15 

16 
17 
lfl 

19 

20 
21 
22 

23 

24 
1 
2 

3 

4 

5 

6 
7 

A 

9 


l o f\j  k*  c;o  u 
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C 10 

CALL  PIStCT  (F0N0l.V,RtSP,A,B,0,MtITM,5,GUtSS)  11 

C 12 

CALL  NEM1CAP  (FONOGV.SONCGV.Pr SP.A.B.C.H, ITH,  MUX,  EPS,  NIJM  ITS  , GUESS  , 15 

1TH:.TA,SOkV,IFAILI  14 

IF  (lFAIL.EQ.il  1,2  15 

C**»  NEWTON  P A FHSPN  OI(J  NOT  CONVFOOE  1ft 

1 CALL  MWTtPK  (THtTA.SOFfV.SfM.TINFO.EXPTOn  17 

KETUWN  14 

C 19 

2 CALL  NOSTST  (H,  IIH,  A,4,i,,THF  M , T INFO , F XPT OT ) 20 

SO^V-AdS (SON  VI  21 

SFH=:.0/SOPT(SOPVI  2? 

VFTMWN  23 

-NO  24 

FIJN'llION  FONOftV  ( 4 - S P , 1 T M, A , n , C , M, THE T A I 1 

I N I r (>fc  K PE3P(M  .PI'.HT  2 

0 lit  NS  ION  A ( H I , 0(41,  C ( HI  , I TM  ( HI  3 

OAM  “I.PIGMT/3.14159?,  i / 4 

II  A f A XMAX,XMIN/7.C,-7.0/  5 

***  CALCULATES  FIi-ST  n-.oi/ATIVL  OF  L OG- L IKEL I HOOO  FUNCTION  OF  ft 

•**  A srsPONbt  VFCIOP  FOP  Tnt  NO°**AL  OGIVE  HOOEL  7 

4 

SL'Hzij.O  q 

k 0 1T“I  = 1.P/0QKT (2. 1*PI)  'll 

IK.  i J=1,M  11 

if  ([14(11. tQ. 01  CO  TO  " l’ 

T r 4 *=A(  I)»(TH.  TA-"(  III  13 

IF  ( I - MP.GT . *MAX ) TFMp=¥HAX  14 

IF  ( I - HP . LT . V-IU)  T(  HP=XHIN  15 

X = - ( I c HP»T  -HP)  /2.i,  1ft 

Ilf!*-  „ I = ;!UOT«I*  A ( 1)  * ( 1 . 3-C  ( I I I *CXP(X)  1T 

III  NC'M  = C ( 1 1 ♦(  1. 0-C(  I I I *C0FN|  (»c  4P)  14 

I‘  ('  i.SP(I).FP.hil(,HI  I 1,0  TO  l 19 

III  on*--  (1.1-0-  NOHI  20 

S114  = [,IJH+ ( 1 (4t  41 /O'- NOMI  21 

G9NllMJr  2? 

f ON ift  1/=5(IM  2 3 

/c  I lift  lv  24 

F NO  25 

FUNCIION  S9NOC V ( "E  F p , 1 T H,  A , » , C,  M,  THE  T A I 1 

INf  b'“  uFb»(M) ,KJ3H1  ? 

014 . N S 1 0 N AIM),  H(4I,  C ( H)  , IT4(h)  3 

DATA  P I , P I GH  T / 3 . 1 1 5 4 2 , f 4 

CAM  XMAX.X4IN/7.w,-7.Q/  5 

***  C A L O' 1 L 4 I E.  3 '•f.COND  Oi  P I V A T T Vr  IF  l OG- L I <l L I H 000  FUNCTIOI  f 

„•**  ft F .1  ktSPO  JSfc  VECT1P  FO  THf  -mi- H«L  OGIVE  HOOF L 7 

u * 

suh^l.o  o 

40JT'-I  = 1.  j/SOcT  (2.5*011  If 

00  * [=1,4  11 

IF  (riHdl.FQ.OI  (O  TO  2 12 

TFHJl=A(ll*(THtTA-9(III  17 

IF  ( I t 4Pl.bT.X4AX  I T t HP  1 = X 4 A V 14 

IF  (lr4Pl.LT.XHIM  T£HPI=XHIN  1= 

x = -f|.  H°1*TEHF1/2.1  1ft 

TF  HiV  = «*OGT*»I*  ( 1.0-0  ( I)  I *A(  II  •FXP(X)  17 

F I P*(UH=  1 E N32  * T E HP.  14 

StC‘il,H=Tl',02*A(II»T‘-HPl  19 

SO'  ,IIM=C(I)»(1.0-C(1)I*COFN(TFHP1)  20 

F 0-  Nf'4=S0r  NOH»S0E  MM  21 

IF  (ki_S°(Il.EQ.?lG-(TI  Go  TO  1 22 

FO--.t.  H=(l.,l-r,OrNOMI  * (l,.i-30r*  14)  23 

SO-*t(  “=- Il.n-STTFNOHJ  24 

1 b(}4=  l|M-  ( F IPFMI H/F  0."  NOMl  - (SEP NOH/SOENOM)  25 

2 CONTI  Mir  26 
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SONH'V  = SU“»  27 

RFT 2 1, 

EN3  2q 

bU  'iPO'  I T T N*-  MCS14I  (^*11M.fl,H,(;,THt'Tft,TI*iF0,CXtlT01  » 1 

DIM-  r _,ljKi  A(M),  H(*»,  C(M»,  IT*MP)  ? 

riaii  pi/,>.  im"p:/  3 

UAIAX1AX'X<1IN/7'.<U-7.0S  4 

TlNFO=i).0  5 

e x“rr'i=c.n 

K OU  M r 0 7 

0 3 

00  i 1=1, M q 

IF  (TTMm.FO.O)  !>0  TO  1 ir 

KOUM  =K0UNT*1  11 

ihMp=c(ii»(TH‘-ra-'(i)i  1? 

IF  (I EIP.bT. XMAX I T t M°= XMA  X 13 

IF  (iHhP.LT.XMlNI  T £ MF*=  X MI  N 14 

p = cui  ♦d.j-cm  )»c:gfn( tempi  ip 

0=1.  r-p  if, 

IFMP=-TfcM°»TF^P/i;.0  17 

PPRIMt  = < 1.  J/SORT  CE.a*PI)»*(1.0-CII»»*«(I)*EXP(TfMP)  18 

T I NF  ri=  T T NF  Of  (PPRIMf  * PPPIME  ) /<°*Q)  IP 

t XPT()I=fc  XPfOTFP  20 

1 CONTINUE  21 

f.  XPT01  = EXCT0T/K0UNT  22 

RETURN  23 

tNJ  ZU 


APPENDIX  D 

ADADSCO  FORTRAN  PROGRAM  LISTING 


PROGRAM  AOADSCO  ( I NPUT , OUT  PUT , OA T A , I POOL  , T A PF 1=0 A T A , T APE 2=1 POOL , T A 
1PE  3 ) 

DIMENSION  ITEM (600 1 t AI600),  B(600),  C(600>,  KEV(60U),  IREJ(«:i,  I 
IF  0RM2 ( 8) » IRAN  180)*  INADS(BO),  IRESP180),  ADM(80>,  BOMISOI,  CDM(80 
21,  bESC(24),  NAME ( 2 ) , IFORM1  (8) 

INTEGER  0PT1,0PT2,0PT3,0PTA 

REAL  ITOT 

N=NC-G 


C * 

C3FAO  OPTIONS  ANO  PROGRAM  PARAMETER  FROM  INPUT  FILE,  DATA  IS  ON  TAP-2  * 


I POOL  = 2 

RE  AO  2 0,  INUP,  MM AX,  OPT  1, OPT ? , OPT 3 , 0 PT4 , T S , T SS , A MA X , UMIN, BMA X , CMA X , 
1IFLAG, IOMIT 


READ  22,  CIFORM1CII ,1=1,8) 

C 1EORM1  IS  THE  VARIABLE  FORMAT  FOR  THE  ITEM  POOL 

REAP  ( IPOOL, IFORM1I  <I1FM(l) ,A(I) ,0(11 ,C(II , KEY (I) , 1=1, INUP) 


C * 

c * 

CS  T ART  READING  THE  SPECIFIC  DATA  (SPECIFIC  FOR  THE  PUN)  F POM  THE  INRUT  * 
C I N A D IS  THE  ITEM  IDAS  ADMINISTERED  * 

CIE'tJ  IS  THi  RtJECTED  ITEM  10  S * 

Cl RE SP  IS  I HE  RESPONSE  VECTOR  * 

C * 


C» 


READ  21,  MNUM, ( IRE J ( I ) , T=l, MNMMI 
KFAl  22,  ( IF0RM2 (I) , 1=1 , 3) 

IF0RM2  IS  THL  VARIABLE  FORMAT  FOR  THE  SUBJECT  DATA 

READ  24,  UtSC 

PRINT  23,  INUP, MMAX, I OMIT , IF L AG,  OPT1 , 0 PT 2 , OPT 3, OP T 4 , T S, TSS , AMA X , GM 
1AX,1MIN,CMAX,IF0RM1, IFORM2,OF  SC 


C * 

c * 

CREAfl  A SUHJECT  FROM  TAr-1  CALCULATE  THETA,  LOOP  BACK  TO  5 ETC.  * 

C * 

C * 

1 READ  ( 1 , I F OR  M2 ) 1 0,  N A MF  , M,  ( I Rc  SP  ( I)  , 1=  1 , MM  A X ) , ( I N A OS  ( I ) , I = 1 , MM  A X ) 
M=MINt)  (M,MMAX) 

I 0=UN I T (It 

C CHECK  THE  END  OF  FILE  ON  DATA  FILE 

C 

IF  (IO.LE.O!  GO  TO  2 
PRINT  25,  10 

2 IF  (10. EG. 3)  GO  TO  19 
N = N ♦ 1 

ITOT=3.0 

0 St  T TH-  ITtM  1JT  TO  ZFBO  FOP  THE  OMMISSIONS  (THAT  IS  RF  AO  IN  FRO 

C 

C AND  SET  THL  RESPONSE  VECTOR  TO  1 IF  THE  ANSWER  IS  CORRECT 

DO  5 1 = 1, M 

C SET  THE  ITEM  IDS  IN  1REEJ  TO  ZERO 

IF  (MNUM.EQ.Ot  GO  TO  A 
DO  3 IJ=1,MNUM 

IF  (iNAOSin  .Nf ..  IRFJ(IJ)  ) (,0  TO  3 
INADS(I)=0 
GO  TO  5 
5 CONTINUE 

4 CONTINUE 

IF  IIRFSPT It .F Q. 10“TT ) TNA  OS  ( T t = r 


1 

2 

3 

4 

5 

6 

7 

8 
a 

10 
11 
1? 
13 
1 A 
15 
IF 
17 

1 A 

IP 

20 

21 

22 

23 

2A 

25 

2F 

27 

2« 

29 

30 

31 
3' 
35 

3A 

35 

36 

37 
33 
30 
40 
A 1 
A2 
A3 
AA 
A5 
4F 
A 7 
A3 
40 
c o 

51 
5? 
c 3 

54 

55 
5G 
57 
53 
50 
6 3 
FI 
62 
F3 
6A 
65 
k 
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IF  UNARM!)  « c (3. 0 1 GO  JO  5 

67 

CALL  SEARCH  < INUP,  INAOS  ( I)  , AOMUJ  ,BOM<n  , CIMC  I)  ,KXt  , ITFM, 

A ,.l , C,  XE  X 

63 

1,  10) 

69 

c 

IF  lnt.  flag  iflag  IS  NOT  ?fro  IT  IS 

T M < F N TO  RE  THE  OJM««Y 

KEY 

70 

c 

IF  II  is  7ePC  INC'I  XEX  IS  pear  from 

THE  POOL  AN  LEFT  IN 

XKE 

71 

IF  (IFLAGaNc . 0 ) KKE=1FLAG 

7? 

leiZ-0 

73 

IF  (IhFSP(I)  ,F0.K<:  ) I«FS=1 

74 

I We  XI)  = I WF  S 

7= 

5 

CONTINUE 

7F 

77 

Q 

• 

73 

c 

• 

79 

CIN 

THr  Nr*(  00  LOOP  I HF  I Tt  M PAW A Mr  TE  RS  COPRESPONOING  TO  THE 

ITFHS  * 

80 

CIN 

IN  AO  rtt>L  WFTR1VF0  F-")M  A,P,C,  AND  LOAOcO  INTO  AOM, BOM , COM 

RESP.  • 

81 

CTHF 

ENTRI»S  IN  THE  A0M,HDM,COM  APE  7EP0E0 

FOR  THr  CASE  OF  Li 

RO 

ITEM  * 

82 

CIO 

IN  THE  1NAD 

• 

83 

C 

» 

84 

c 

• 

89 

c*** 

# ♦ 

86 

A 

If  OP  1 ION  3 IS  ON  THE  PARAMETERS 

A,H,C  APE  CONSIPAIN>() 

WITHIN  b 

87 

C 

OF  AMAy,AMIN,HM4X,RMIN,CMAV 

83 

IF  (UFT3.EQ.0)  GO  TO  7 

89 

00  F 1=1 ,M 

90 

IF  (TNAOS(I) .LO.O)  GO  10  b 

91 

IF  ((ADM(I).GT.AMAX).OW.<OPT3.EU.3)> 

A OH ( II =AM»x 

92 

IF  (HUMI1) .LT.BMIN)  B0MII)=9MIN 

93 

IF  (HTMII) .GT.HMAX)  BOM(I>=RMAX 

94 

IF  ( 1 COM ( II .GT .CP  AX ) .OR.  (OPT  3.F0.2) ) 

COM( 1) =CMAX 

96 

6 

CONI i NUc 

96 

C 

OPTION  1 WILL  P^jNI  IMF  SPrOIFIC 

nA I A IF  ITS  ON 

97 

C 

98 

7 

CO NT TNUF 

9Q 

IF  C (QPTl.FO.ru  .op. (0P1  l.GT . 10)  1 GO 

TO  3 

100 

qo  f 1 = 0PT 1*1 

101 

° P I N 1 26,  10, NAMf  , ( INAOSI T II  . rJESR( I II , AOM( II 1, 00M( ri 1 , 

COM! IT), 11= 

102 

11,1) 

103 

3 

CONTINUE 

104 

IT  I NA 0S=  G 

105 

no  o <<=i,h 

106 

1T0T=IT0I fIPESP(KK) 

107 

9 

IF  (INAOS(KX) .EQ.CI  iriNAOS=lI INAOS+ 

1 

103 

KL=m-ITTNAOS 

109 

I TUT  = 1T0I/FL0AT  < (XL  > ) 

110 

if  (dot. to. o)  10,11 

111 

10 

PVIM  Z7 , NAME, IT 

112 

IF  (OPT 4,60,3)  GO  TO  13 

113 

r=- io.o 

114 

SFOt'M  = 0.  0 

115 

r INFO=0.  c 

116 

EXPTOT  = 0 .(1 

117 

SEM=C.O 

113 

I T h W=0 

119 

I X L = XL 

120 

GO  TO  13 

121 

11 

IF  ( HOT  .EO.  1.  0)  12,13 

122 

12 

PPINI  31,  NAME, 10 

123 

IF  C0PT4.C0.3)  GO  TO  13 

424 

r=ia.o 

125 

SFO=>M=0. 0 

126 

r info=o.u 

127 

E X°T0T  = 3 .0 

128 

SE  M=0 . 0 

129 

ITtR=0 

130 

X XL  = XL 

131 

go  rn  i« 

132 

c • 
c * 

C NOW  THE  D*TA  IS  READY  TO  MAKE  THE  CALLS  TO  THF  APPROPRIATE  ROUTINE 
CTO  ESTIMATE  THE  THETA.  OPTIN  4 MILL  DETERMINE  THE  METHOO  BY  WHICH  • 
CTHt  T HE  I A ESTIMATE  MILL  BE  fOUNO  • 
C * 
C * 


I KL=Kl 

IF  (0PT4-2)  14,15.16 

CALL  HAXLNO  ( IRESP , I NADS, H, M , AOH , BOM, CON,  50, . 1 J 5,  T , SF ORM , IF  A IL  , T I N 
lFO.tXPTOT, ITER, SEMI 
GO  TO  IT 

CALL  MAXLK  (M, I NADS, IRESP, M, AOH, BOH, COM,  5 0, . 0 "5  , IF  A IL  , ScOCM , T,  T INF 
10, F.XPTOT, ITER, SEMI 

t;  0 TO  IT 
T = I S 

SF  0rH=T3S 

CALL  BAYES  (M,  INAOS,  IRESP,  M,  AOM,  BOH, COM,  T.SFOP'i,  T INFO,  "XPTOT  » 

I TtR=0 

Sf  m = 0 . 0 
GO  TO  1* 

IF  (IEAIL.EQ.0)  GO  TO  18 
pf-INI  28,  NAME, ID 

SF9RM  ANO  T ARE  SET  TO  -99. 90  IN  MAXLK  IF  NOT  CO  NY 

NC  = NC ♦ 1 


CWRITF  THL  SUCCESFUL  RESULTS  TCTHE  FILE  T APE  3 


W <•’  1 Tt  (3,29)  NAME,  ID,  ITOT,  T.SFOKM,  IKL,  TINFO.CXPTOT,  IT:  2,  SFM 
GO  TO  1 

“PINT  30,  N,  NC 


FCIMAT  (214, X, 411. 6F5. 2, 12, 121 
FORMAT  (1615) 

FORMA!  ( A A10 I 

FORMA)  ITSa,*ADAU3CO*,/,T50*=======*,//////,l  .'’0  , * AOAP)  I Vc  DICHOTOM 

IDS  SCORING  WITH  THREE  PAPAMF  TF  P MO OF  LS’,////» FrU’PSYCH  )MF tpic  M<TH 
2303  PROGRAM*,/,  T40»DEPAPTME  NT  OF  PSY  CHOLOGY  * , /,  I 4 0*'JN  I YE  PS  I T Y OF  m 
31  N Nr  SOT  A*,/,  T4  0 • MDL  S • MINN.  5S45 5 * , / , / // , T2 Q • INUP * , T 2 T * = * 15 , /, T 2 L , 
u*mmAX’,T27’=*I5,/»T20,*IOMIT*,T27*=*,I5»/,T2G’1FLAG*,T2T’=*,10,/,T 
62 )*OOT1*,T27*  = ’I5,/, T20’OPT2*,T27*  = *I5,/, T2  ’OPT  3*, )2/’  = *IS,/,T20’ 
6ORT4*,T2T*=*,I5,/,T20*TS*,T2T*=*F6.2,/,T20*TSj*,T2T*=*F5.2,/,T20*A 
TMAX*,T2T*  = *F5.2,/,  T 2 0,  • BMA  X*,  T 2T,  ♦ = • F5. 2, /,  T2n’PMIN»,t27,*  = *Fr-.2,/ 
“»T2&’CHAX*,T27,*=’F6.Z,/»T20, ’VARIABLE  FORMAT  FOR  PUOL=*HAli,/,T2D 
9 » * VAw I ABLF  FORMAT  FOR  0 A T A = * A A1 0 ,/ , T 20 , 8 A 1 0 , / , T 2 C , A A 1 u , / , T? I , 6 A 1 0 , 
F / ) 

FORMAT  ( SA10 ) 

FORMAT  (10K, ’PARITY  ERROR  ON  T APE* ,9  0X , I 2 ) 

F ORMAT  (10X,  A9,  2A10,/,  CX,  I4,?X,  I2,2X,  JF10.  2)  ) 

FORMAT  (10X,’  SUBJECT  ’,2AlO,*IO=  ’,A9,’HAS  NO  RIGHT  ANSWERS*) 
FOPMAT  (10X,»  MAXIMUM  LIKELIHOOD  ESTIMATION  DOES  NOT  CJNVfPGF*,*  F 
10R  THE  SUBJECT  = ♦ , 2A10 , * 10=  ’,A9) 

F ORMAT  (X.2A10,  A9,F8.2,?F7.2,I4,2F7.2,  14,  FT.  2) 

FORMAT  (lOX,*  CAStS  KEAO=*,I5,’  CASES  NOT  CONVF  pGi 0 = ’, 161 
FORMAT  (lOX,*  SUBJECT  *,2A10,’ID=  ’,A9,’HAS  ALL  ANSWERS  RIGHT*) 

F ND 

SUBROUTINE  SEARCH  ( I NUP  , 10  , A ,0,  C , KE  Y , I TM,  AO,0P,CP,<'.Y  J,  I0NU1) 

I)  I Mt  NS  I ON  AP(INUP),  HP(TNUP).  CP(TNUP),  KEY-(T*)P).  ITKT(III’) 


ooo  oc?r>  m o O o 
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INTlI^R  FLAG  3 

DO  1 i=l,TNUP  4 

If-  ( 10.  Nt  . IT M ( I n GO  TO  1 6 

a=ap(ti  6 

H=9P(I)  7 

C=OP(I)  6 

«EY=KEYP(I)  *3 

CALL  PCHtCK  (£,H,(  , ICMIONUNI  11 

rfiurn  li 

1 CONTINUt  12 

PHI.. I 2,  1J.IDNUM  13 

C 14 

ID=J  is 

k>-  pirn  16 

17 

16 

19 

FORMAT  (10X,*  I1M  =*,Iu,*  IS  NOT  IN  THE  POOL  F OR  SUBJECT  10  = *,A9  20 

if  21 

E NO  22 

SUPPOUTTK-  =>CPLC<  I A,”,  0,10,  IPNOH)  1 

•**  CHtCX  WHt  IN:  P 0*  MT  1 1 r (•  P ARAMEPRS  AR_  VAl  in  2 

*•*  IF  NUl,  ERROR  MESSAGE  10  PPIMFO  3 

4 

IF  (A.LE.O.QI  °R1;T  1,  IO.A.inNUH  6 

IF  ((-6.U.bT.P».OP.(B.GT.6.0M  °F  I NT  ?,  IP.O.IJNUM  6 

IF  lU.C.U.n  .OF  .(C.r.l  .1.  of  » PPINT  3,  IO.C.IDNUM  7 

RETURN  4 

q 

10 

11 

1 FORMAT  (10X,*ITEM  =*,I4,*  HAS  THf  INVALlO  A PARAMETtR  OF  *,F6.2,*  12 

1 A MUST  fc  GPEAli  7 THA4  0 . "* , /1TX , * f PROP  FOUND  *,*for  THC  SUBJECT  13 

2 WITH  10  =*,A9I  14 

2 FORMAI  (1PX,*TIEM  =*,I4,*  HAS  Th l EXTREME  fl  PARAMETER  OF  *,F6.2,5X  16 

l,*tc:ROR  FOUND  FOF  THE  SUBJECT  WITH  ID  =»,A9I  16 

$ fopmai  <icx,*itfm  =*,i4,*  has  thf  invalid  c parameter  of  *fF6.2,*  17 

1 C MUST  HE  RFTHE-N  0.0  AND  1 . r * , / 1 C X , * - R ROP  FOUND  *,*FOR  THE  S'JBJ  16 

2tCT  WITH  10  = *, A9I  19 

EMU  20 

SUBROUTINE  BAYES  ( M , 1 1 M , RF S» , N, A , " , C , U T HE T , MV AR , T 1 N FO , E X PTO T I 1 

INTFl.fR  RFSP(M)  ,ITM(M)  2 

PFAL  1(NI,HIN),C(M  3 

DO  1 1=1, W 4 

IF  UTM(U.EO.O)  (.0  TO  1 5 

CALL  HSCO°  (HTHET ,HVAp,  MI » , A( Tl ,C(I ) ,RFSP( I) > 6 

1 CONTINUE  7 

CALL  NOS  I A T ( M,  I T M , A , H,  C , H THt  T , T 1 NF 0 , E X3T  CT  I 6 

p F T UK  N q 

END  10 

SUHPOUTINE  OSCOR  ( 0 THE 1 , BV AP , LUF , 013, GU 'S P, IRESP » 1 

0=  (DIF -HI  H£T  | XSQF'I  (2.0*(1.0/OTS**2*DVAR»)  2 

ERFD=ERFNO(0)  3 

EQS0  = tXO  (D**2)  4 

IF  (H)SO.fa.O.OI  RFTUPN  6 

E QS(Jl  = l , O/EDSO  6 

XKINV=0.5*  (1.0-ERF  3)  7 

XLINV  = GUESP» (1. O-GUESP*  *XKINV  ft 

IF  ( (YLINV.E  0,0,0) . OR . (XKINV.f  0 . 0, 01  I RETURN  9 

XL=1.0/XLINV  10 

IF  (TRESP.NE.il  GO  10  1 11 

S=0 .3989  42* (SORT (HVAR) /SORT (1,0* ( 1 . 0 /D 13**2 1 /9V AR I I • ( 1 . 0 /XK INV) *E0  12 

1S0I  13 

r = 1 .0-1. /72454*D*rOSO*( l.O-FPFO)  14 

a THE  ( =9 THF  T*  ( 1 . O-G'JE  SPI  *XK INV**L*S  16 

BVAK=DVAR- (1.0-GULSP»*1ICINV*XL*S»*2*(T-OUESP*YLI  16 


A 
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I 


RE  TURN 

1 BTHET=BTHET-0. 7 97885* (BVAR/SQRT ( 1 . O/OIS” 2 + BVAR I I *£0501* II. 0/C  1. J* 

ItRFOII 

PART  1=1. 128379/11.0* (1. 0/OIS**21* Cl. O/BVARI ) 

PART2=1.0/(EDSQ*I1.0*ERF01 

PART3=0.564190*D*EOSQM1.0*ERFOI 

BV AR=BVARM1. O-PART 1*PART2*PAPT31 

RETURN 

ENO 

REAL  FUNCTION  ERFNP  (X> 

DATA  Al/0. 254830/ 

DATA  A2/-0. 284497/ 

OATA  A3/1. 421414/ 

OATA  A4/-1. 453152/ 

OATA  A5/1. 061405/ 

OATA  P/0.327591/ 

ERFNP=0. 0 

IF  (X.EQ.0.01  RETURN 
ES=SIbNC 1. 0. XI 

y = abmxj 

IF  (T.LT.6.01  GO  TO  1 
ERF  NP=ES 
RE  TURN 
1 T2=Y*T 

T=l.U/tl.0*P*Y» 

AT  = ( I Al*  |A2* (A3*(A4*A5*T  * * T > *T1*TI*T1 

t Af=AT/EXP|Y2l 

EPFNP=(1.0-EAT>*ES 

RE  TURN 

ENO 

SUHPCUTINE  N AXLK  I 9, I T M , RE SP , N, A , B , C , H A X , ' I l • SDR V* THE T A , T INF 

iO.fcXOfOT.NUHITS.StNI 
t X TFRNAL  FOOLOG.SOQLOG 
INTEGER  RE  SP (Ml 

D I 9t  NS ION  AIN),  BIX),  CCN),  IT9IM) 

C***  USt 5 MAX1HU9  LIKELIHOOD  LOGISTIC  SCOPING  ALGORITHM  AND  Rf S PONCE 
C*»»  HOOrl 

C**»  RI SECTION  IS  USED  TO  PROVIDt  THE  INITIAL  GUESS  FOP  THE 
C***  NF  H T ON-R  APMSON  HE  MOO 

CALL  BISECT  |FOnL03,RES“,A,R,r,M,ITM,5,GUF.SS> 

CALL  NENTRAP  (FDOLOG.SDOLOG.RF  SP, A,R,C,H, ITH, HAX.FPS, NU NITS, GUESS, 
1 f Ht T A , SORV, IF A IL ) 

IF  (IFAIL.EQ.il  1.? 

Nf  H TUN  RAFHSON  010  NOT  CONVERT, c 
1 call  NNT t RR  ITME I A, SDRV.SFH, T INFO.EXRTOT) 

RETURN 

CALL  LGSTAT  CH, 1TM,A,B,C,THFTA,TInF0,£XpT0TI 
SF 9=1.0 /SORT  I ANSI  SOP  VI  I 
RETURN 
k ND 

FUNCTION  FOOLOG  ( »E S P , I T N, A , H, C , H, THE T A I 
INTEGtR  rfspchi .right 
DIHfNSION  AIMI,  Bin,  C ( HI  , IT9I9I 
OATA  XNAX.XNIN/200. 0,-200.0/ 

DATA  0, RIGHT/1. 7,1/ 

C**'  CALCULATES  FIRST  DERIVATIVE  OF  L OG-L IKEL I MOOO  FUNCTION  OF  A 
C»*»  PFSPPNSE  VECTOR  FOR  THE  LOGISTIC  HOOEL 
SUHiU.O 
00  1 1=1,9 

IF  IIINI1I .EQ.OI  GO  TO  1 
X = 0*» (II • (THE TA-BI 1 1 I 
IF  (Y.LT.X9INI  X=  X9 1 N 
IF  (X.GT.X9AXI  X = X 9 A X 


17 

18 

19 

20 
21 
22 

23 

24 

25 
1 
2 

3 

4 

5 

6 
r 
8 
9 

10 

11 

12 

13 

14 

15 

16 
1 7 
18 

19 

20 
21 

1 

* 

4 

5 

E 

7 

A 

9 

n 

11 

12 
13 
1 4 

15 

16 
17 

1 A 

19 

2', 

21 

22 

23 


4 

A 

Ft 

7 

A 

Q 

10 
11 
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t Xr  = *•  XP  (X  1 14 

At=A(I)*EXF  15 

SUN=SUH-AF  /(tXF+1,9)  15 

IF  (RESPm.NE. RIGHT)  G"  TO  1 l* 

CE  =C ( 1 1 ♦£  XF  in 

sum=sum*  Ar  /ce  1^ 

1 OO NT  1 NUF  2P 

FDUL0(.*-1.7»SUM  21 

Kr  TURN  2^5 

t NO  2 3 

F UNO  I ION  SOHLOG  ( <ESP, ITM, A, H.C, M, THlTAI  1 

INI  l,rS  F'E  >P<MI  ,RISHT  2 

Oil  •'<  j 1 0 N ITM(M),  4 ( M ) * e(M),  C(M)  » 

04ra  XMAX,XMIN/20iT.  0,-230.0/  I. 

o 4 r a it, right/: .7,1/  5 

0***  CAt.CULATtS  SECOND  DtRIVATIVf  OF  LOG-LIK'LIHOOO  FUNCTION  5 

C***  OF  4 RESPCNSC  VECTOR  FOR  THE  LOGISTIC  MOOEL  7 

SUM=l.O  a 

00  i 1=1, M 9 

IF  (MHdl.EU.OI  GO  TO  1 IT 

x=i)*am  *uHti  A-em  ) n 

IF  (V.LT.XNIM  X=XMIN  12 

IF  IX.GT.X9AX)  X = XB AX  13 

t XF  =- XP(XI  14 

AF=A(I)*tXF  15 

SUM  = V(!H-AII)  *At/»  Il.OFfcXFIMl.fmXF))  IF 

IF  T R t.  S ° T 1 1 .Nf  .RIGHT)  GO  TO  1 17 

CE=C(I>*fc*F  18 

SUM=SUM+ ACII  *C<  II»AE/(Cr*CE>  19 

1 CONTINUE  ?0 

S0DLUG=-2. 99*SUH  21 

RETURN  22 

E NU  23 

S UHROUTI NF  RISECI  Tf  l,Rf SP , A , 9 , C , M , I T M , N I TrR , ON  I 0 » 1 

INTEGER  RESP(M|  2 

DIM*  NS  ION  A ( h I , B(H),  C ( Ml  , ITM(M)  3 

C***  CALCULATES  APPROXIMATE  ROOT  OF  FI  9r  BISECTION:  4 

C***  BISECTING  NITER  (NJMBfcR  OF  ITERATIONS)  IIM^S.  5 

C***  BMIl)  IS  B'.ST  CURRENT  GUESS  AT  ROOT  THtTA  5 

C 7 

C***  INITIALIZE  LEFT  HOUND  A NO  FI ( BOUND I ANO  RIGHT  BOUNO  FI (BOUND)  8 

BL  = -b . 0 9 

8K=5.0  IP 

BMI(1=G.O  11 

T L = F 1 CRES°,ITM,A,H,C,M,OL)  12 

TR  = F1  (RESR,ITM,A,8,C,M,0R)  13 

C***  TEST  FOP  NO  ROOT  IN  INT F RV AL --RE TURN  IF  NO  SOLUTION  14 

IF  HTL*TP)  .GT.O.OI  RETURN  15 

0 IF 

C***  NOW  CALCULATE  BISECTIONS  NITFR  TIMFS  17 

00  S 1=1, NITER  18 

I MI0=F1 IkF  5P,ITM,»,B,C,M,BMID)  19 

IF  ((IMIO’TH.GT.O.O)  GO  TO  1 20 

c***  R E PL 4CE  PIGHT  BOUNO  HI1H  BhIO  21 

BR=TMIO  22 

GO  TO  2 23 

C***  REPLACE  LEFT  BOUNO  WITH  BMID  24 

1 TL=TMIO  25 

BL  = BM 1 0 26 

C***  FIND  NEW  MIDPOINT  9MIO  27 

2 BMIO=(BLFBR)/2.0  28 

3 CONTINUE  29 

RETURN  30 

t NO  31 

SUBROUTINE  NFWTRAP  (F1,F2,RFS°,A,9,C,M* ITM, NITER, EPS, NU HITS, GUESS,  1 

ltHtTA.SORV.IFAIL)  ? 


k. 


O C i vi 


r 


INI  -l.  M S P I M I 

HIm-hsion  AIM),  Hill,  CIM,  II  1 ( *11 
C*»«  CALCULATE  S THt  KO'>T  OF  FI  GIV  N 1 1 3 FIRST  'V9IVATIV-  F2 
q.».  (SMI)  AN  INITIAL  Gll‘  SS  USTNG  Nl W TON - 5 A PH SO N MFTHOP 
c»**  I H r T A 13  APPR.  TO  THE  ROOT  5 SO^V  is  Ft(THfTA) 

NUMI I 3 = 0 
I Ht T J = GUt  * 3 

C***  LOOP  UNTIL  *'  k F « F p . » Ok  NUKOto  "F  IT-P1TIONJ  OFOONtS  TCO  L A FGr 
1 F OR V - F 1 ( V JP, IT H,C, «, THr TA) 

SP‘>/-F2(kl  3P,  ]TN,A,.«,c,  1,TH*  T A ) 
t c tzy  'tpy/sOP  V 
I M*.  T„  = THt  I A-t  P.N 
.0“  ; I S = N 1 1 M I T 3 ♦ 1 
2*»*  .nil  LOOP  C-  IT  F k I 'I 

If  ( | •lij'lp  S.L  I .NIT  _~)  . A‘  a.  I Ai<«.  (•  F’t -I  ,GT  .F  DS)  | r,0  TO  ; 

C*"  t NO  LOOP.  I * I I H,->  CONV  RGLl.C  A \ 0 Sf  T 1 F A I L 

I F A I L — 0 

IF  (a  IS  It  -R>  ,L  f ,L">»  F r TIPN 


0**»  Nf  U • N RAFH'O'i  HE  I Hon  miS  NOT  ClNI/i.»>GE 
r F A : I = 1 

• L I -*'N 
l NO 

-,H  I-  >.  ft*.-  • WT  : kh  ( I H-  I • ,0f  0 .i  V,  1 INFO,  I f >T  nT  I 

",»*♦  , • I OJI,  V-l'ir.G  FGk  I LAS  ir.  WHICH  NiWTJ  -'AJM_'_o<  ‘ATI') 

o» • * ip  ••  y s r , : 

i h - • ..  = - ■>  i . E 1 

.F  > <->.  )i 

Si  » 

I IV  = -OH.  h ’ 

t < ; ' | = - H l . 1 C» 

> • I ' • ' 

- * 0 

Jill  '!!*.  LGSTAI  ( O,  IT  % A , T , ] h TA.TI'FO,  X pT  II 
'•  : M I ' f.  A ( M)  , < T I , C I Ml  , I ’ HI  • 1 

oa  i »m;  x . *"  r / 1 u . 1 . - i ' . >./ 

II'.»  *0.0 

tx*'  1=1. 

< TIM  =0 

C 

i)  o . !*;,* 

IF  |,|i(il  ,o,oi  t c : 

KUO  I =K  0C‘  I ♦ 1 

«c  ."--l.rUI'M  I"-  T A-  » ( III 
jr  ( Gil  ,M  . VKAX  I UGU*ykAx 
If  I - GU  . 1.  T . X •‘INI  ARG'IsX'*!-, 

••=<:  i * ) ♦ i i-c  m ) 5 ( : . / 1 1 . n * ..i))» 

i-  i . -p 

* A-S.  = t XP  ( - A ■’  N"  ) 

Pi'-;*  AMi/Kl.Gl  ARIl’C.r*  A l l) 
i * . * - f . =PP'  IHt  » ( 1 . ’-0  ( 1 I * »A  ( I ) • ! . 7 

I I N-  = T 1 1 F .)♦  (►  PR1  I-  * P-).  i ».-  | / ()  * II 

; x I'  I I = " X p T o T * P 

i c o n r i *.i|F 

► x > i i =f  x- r ot /nu  . r 
c * r 1 1.  n 

t'JII 

Slj'p  ol  UN  k t X L N L ( - t ' H.  IT  H,  •*.  N,  t.  1,0,  1AX  ,f  ’S«TM*  T A , S OF  V,  l F A ,t  , I I *! 
IF  0,  X J I O I , Nl  I K I f S * ;■"«) 
r X I ‘ nal  f UNOGV,  SO'ICGV/ 

I M ■ (..  p kr  JP  ( 'll 

nil  ’SION  ITHINI,  MNI.  MM,  0(N* 

I J • , p . *‘AX  l**i|M  LKlLIHOOl  N()PMH  iT.Wt  SCOPING  ALGORITHM  AN!' 

• * t v-  ,*'(•. ‘.Sr  </ " C T C R 

**•  II  I.IIIN  I*  l.jFN  10  f'V  "VI H T*II  T I AL  G'l-SS  rf!c  TH' 

MWIIN  P A F H'lfll  HtT-«on 


k 


A 


X 

4 

c 

f> 

7 

II 

q 

in 

ii 

l? 

IT 

14 

I G 
IN 
17 

II 
10 
2r 
21 
2? 
2T 
74 
?F 
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4 

r. 

7 

Q 

1 1 
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U 

r 

C. 

7 

s 
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4 *) 
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n 
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ic 
i7 
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C 10 

CALL  hlStCT  (lDNO,.V,fitSt*,A»‘’.0»M,ITM,5,r-URS'3»  11 

C 12 

CALL  KEHIRAP  <F0Nn3V,S0‘  CG  V,  RF  SP  , A , H , C , M,  IT N, HAX , rP3, NUN  IT S , CUE SS  , 13 


1 THE  ro.SPFfV,  IFAIL)  14 

IF  (lFAIL.ET.il  1,2  15 

C***  NERTUN  RAFHSOi.  OIO  NOT  CONVERGE  If, 

1 CALL  FWTFpr  ( T«E1  A,  <50RW,SEM,  TINF0,E*',T0n  17 

RETURN  13 

C 19 

2 CALL  (OS  I AT  ( M,  1TP«A,H,  °*  THFT  A,TINFO»EX°1  ()T  1 2(1 

sof  v=  abs  nn»v>  21 

3tM=:.n/SORT(SORVI  2? 

RFTU-’N  2 3 

f NO  24 

FUNGI  ION  FONOGV  < Rc  S P , l T M,  A , B , C,  M , TH<- T A I 1 

INfrirP  fr^PlMl , WIGHT  ? 

DIMENSION  A ( Ml  , BIN),  C(M),  ITMCMI  .3 

D A I A PI, K1GHI/3. 141592,1/  4 

UA1A  XNAXtXMIt4/7.0,-7.0/  5 

C***  CALCULATES  FIRST  0 RIVATIVE  OF  L OG-L  I*t  L I HO  TO  FUNCTION  OF  F, 

C***  A RESPONSE  VICTOR  FOR  I HE  NORMAL  OGIVt  NO  Dc  L 7 

C 3 

5 U M = „ , 0 9 

RO'JTt->l  = l.e/SQ0T  (2.  J*PI)  10 

00  2 1=1, M 11 

IF  (1  mill  .£0.01  GO  TO  ? 1? 

r t NP=U ( 1 1 * (THE T A-B<  I ))  13 

IF  (IrMP.GT.XNAX)  T t HP-  X HA  X 14 

IF  (leMP.LT.X«lNt  TFM°=XPIN  15 

X=-(  ItHP’TEfiPl /2.1  15 

ONHffA I = ROOT0 1 * A 1 1 1 • ( 1 . 0 -C ( 1 1 ) *EX  P ( X)  17 

Of  NON=C I II  Ml.  O-CII 1 1 »COFNlTENPI  18 

IF  (klSP  111  .EQ.RIGH  I)  0(1  TO  1 19 

OF N0H=- ( l. O-OFNOM)  20 

1 SUM=SUNF  II1NMRAT  /O'.NOMl  21 

2 CONTINUE  22 

FONO(.V=SUM  23 

RETURN  24 

t N 0 25 

Ft) NCI  ION  SONOC.V  I fc F S P , I T M,  A , H , C,  M,  THE  T A 1 1 

1 N T c G t R K t 3P I Ml  , FIGHT  2 

0 I MENS  ION  A ( M ) , 0 ( H I , C ( H)  , ITNIMI  3 

U A 1 A PI, RTGHT/3. 141592,  1/  4 

0 A I 6 XMAX,XMIN77.0,-7.0/  c 

CALCULATti  StCONU  OtRIVATIVF  OF  LOG-L IKr L IHOOO  FUNCTION 
OF  A RESPONSE  VECTOR  FOR  THE  NOPHAL  OGIVE  HOQFL 


SUM  = |j.O 

ROOTPI=1.07SQ°T (2. 1*PII 

00  2 1*1, N 

IF  <1TN<  11 .EQ.01  GO  TO  ? 

TEMP1  = A< II *1 TMETA-Hf I)  I 
IF  (IENP1.GT.XMAX)  TEHP1=XHAX 
IF  (ItMPl.LT.XMINI  TEHP1=XHIN 
X = -TF  HP1 • TEMPI / 2. 0 

T E MP2=R00f  PI*Cl,0-C(I))*A(I)*fXP(X) 

F IRNUH=TEHP2*TEHP2 
SEC NUM=TEMP2*A III ’TEMPI 
S0ENOM*CIIIMi.0-Cm  l*rOFN»TFNPll 
F OF num=soenom*sde  nth 

IF  (ctSP(I).EO.RIGHTI  GO  TO  1 
FDt NUH=  U.D-SOENOHl  • C 1 . 0 -S OE NOMI 
SOI  NOM  = - (l.O-SOENOU 

SU  M -SUM- I F IP  NUH/F  OE  NOHl - ISECNIIH/SOENOH) 
CONTINUE 
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SPVHjV* 'll"  27 

*•  •?  f 'lw N 21 

lMI)  ?<4 

s')  <-i'HTii«r  r.o'ijar  jm, ii ",a ih:  ta,t info, r xmon  i 

Dll  ’ISI3N  Mm,  t'Cli  tt(H«  IT1I1I  ? 

041'  D14X,  XII)  /7,.,-7.sl/  5 

Til..  PI/  1.  (. 

ITUF"  = d.C  r- 

I X r>  r 0 I * J . J 0 

KTll|:1  7 

C 1 

no  . 1=1,"  <? 

IF  (110(1)  .M3, 01  (-1  TO  « 17 

x rii‘:i  =<0(jx  T*i  11 

l'Mr.(l)*(TH114-r>(Il)  1? 

IF  (I  03  .l.T  . V I4X  ) TrMMMway  IX 

rs  (i  • m^.i  r.  x n n » r f y*t  ?4  i<. 

psr:m*ti.3-cmircFMTi  "bi  ir 

i!=:.  -®  i* 

j.v,  - 1 ; MF  *Tf  mp/,'.  : i7 

w n ] -•  = ( i . 1/Sf  » 1 (,  . 1 * = l » 1*  ( 1 . )-  7 ( i 11  ( T ) *F  ("  (TM,0I  1“ 

r r .t  = t i f>  )+  ( » pi>ip i «•- j / ( 

r y .»)  1 *f»t.  10 1 c 1 

1 i.  O .1  i .0  21 

i . » * i m = - *>  i , i /koiimt 

i i-  • ?.? 

1 ?u 


APPENDIX  E 

LINPSCO  FORTRAN  PROGRAM  LISTING 


l 


C 

2 

c 


c 

c 

c 

c 

c 

r* 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c; 


pf  ogl  a 'i  li  sf’cco  ( ;.\pui  , out  put,  oat  a,  ibool,  tapei=oata,tape  2=ipoul,t  a 
i»E  .1,  PUNCH) 


DEFINITION  OF  VAOIARIFS  USED  IN  LlJPSCO 
I = T H t 1 A * THF  ARIL  I T T ESTIMATE 

ITEM:  AN  APkAY  O'-  IHl  ITEM?;  IN  THE  ITER  POOL 

IN'JPt  Tf-F  NUMREF  IF  ITEMS  IN  THF  ITEM  ^OGL 

At  AN  ARRAY  OF  !H-  OISORI MINATION  PAF4MF  IF  -MS)  ASSOCIATED  WITH 
tACH  IIFM 

Pt  AN  ARRAY  "F  1 Hf  OIFFJCULTY  nAp  AME  TFRS  ASSOCIATED  WITH  EACH  ITEM 
MAXCuTI  NO.  OF  k-fpQNS  CATEGORIES  MINUS  1 FOP  THt.  ITEM  IN  THE 
I T r M POOL  WHICH  HAS  T Hr  MAXIMUM  NO.  OF  RESPONSE 
CATEGORIES 

INAnj  flK  4PKAY  OF  THr  ITEMS  WHICH  HAVE  BEEN  COM  I N I S T t.  RE  0 TO  ALL 
SUBJECTS 

Ml  HE  NUM«EP  OF  ITEMS  IN  THt  TEST 

AOM  A NO  QOMt  THt  ARRAYS  OF  PARAMETERS  ASS3CIATru  WITH  TH-  I TIMS 
IN  INC  0 

INAOSI  AN  AP"-AY  OF  THt  ITEMS  U5FO  TO  CALCULATE  THcIA  FOP  A GIVEN 
SUBJECT  T INA0S=  1N.A0- (PE  JFCT  ION  A NO  OMISSIONS)  WHERF  10  S 
OF  REJECTED  OP  OMITTEO  I Tt  MS  APE  SET  TO  ZERO) 

IRFSPt  THE  PlSPONSE  VtF TOP  OF  THE  CUPptNI  SUHJFCT 
NCATl  AN  APkA  Y OF  THr.  r UM9t  P OF  RESPONSE  CAT  tGOF  It  3 MINUS  i FOR 
rACH  ITEM  ADMIMNMRtO 

Pt  P CN  T I '-rP  CENTAL  OF  ITEMS  ANSWERF  0 HY  A SUBJECT  FOP  WHICH  THF 
TEST  FrSPONSt  HAS  BE- N CHOSEN 

01  A CONSTANT  US  0 TO  CHANGE  THE  METRIC  OF  THF  LOGISTIC  GRADED 
MODEL.  0=1. 7 O"1  1.03  THF  OrFAULT  VALUE  IS  1.0 


COMMON  0 

u i Mt  ns  ion  item ( i lo • . ac;og,ioi,  p(1og,i:i,  inaouooi,  irejugjt,  i 

1PESPI100I,  AOM TlCJ.lO),  80M (100,10),  NAMtCM,  IFORM(8l,  NCAT  (130) , 
2 I F OP  Ml  (81  , INAOSdOOl,  OE  SC ( ?4T 
INTi  GtH  U°Tl,C»T2,0'>T4,AOIM,BOIM 


N=NC=C 

c * 

CPE  AO  OPT  TUNS  ANO  PROGRAM  PrtP  AMt  TEF  FPOM  INPUT  FILE  DATA  IS  ON  T APEC  * 

C * 

c,,, 

I °OOL  =2 

READ  TO,  INIJP,  M,  MAYCAT,OPT  1,  OPT2,OPT4,  0,  IOMTT 
CALL  CHK1N°  (INUP,M,OPT4> 

IF  (ll.  NE.  1.71  0=1.0 

PE  AO  33,  (IF  OR  HKII  ,1  = 1,8) 

PE  AO  J 2,  (NCAT ( IT, 1=1, MI 
C***  I N PI)  I ITEMS  ANU  THEIR  PARAMETERS 
IKI  = AUIM(MAXCAT,0*“T4T 
K=BOI m(MAXCA  T, OPT.I 
00  1 I = 1 , I NUP 

1 PEAU  ( I P00L, IFORM 1 1 I T E m ( I » , ( A ( I , J J| , J J = 1 , IK  I T , (fl ( I , J » , J= 1 , <) 

C 

RE  AO  31,  ( INAD( IT , 1=1, MT 

READ  32,  (IPEJ(I),I=1,MI 
READ  33,  (IFOPMTIt, 1=1,81 
READ  34,  OESC 

PR  INI  36,  I NUP, M, 10MIT, OPT1, 0PT2,0PT4, MAXCAT ,0, IF  OR Ml , IF QPM , DFSC 


L 

C**»  IHOSt  ITEM  10  S WHICH  A-’E  IN  THE  REJECTION  VECTOR  WILL  BE  SET 


TO 

ZERO.  ZFF’O 

ITt M in  S WILL  BE  SXIPPlO  DURING  COMPUTATIONS 

00 

2 1 = 1, M 

z 

IF 

( IRE J (I ) ,tO. 

IT 

I N AO ( I ) = 0 

IF 

IOPT2.Fa.1T 

GO 

TO  10 

C 


1 

? 

x 

u 

5 

6 

7 

8 
O 

11 

11 

1? 

13 

14 
IF 
16 

17 

18 
IB 
2 a 
21 
2? 

23 

24 
2F 
26 
27 
23 
20 

30 

31 
3? 

33 

34 
36 
3F 
3? 

38 

39 
41 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 


I 
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t 


c 

67 

SEARCH  FOR  ITFMS  ADMINISTERED  IN  POOL 

bfl 

DO  9 J=1,M 

69 

IKI=AOIM(NCAT (Jl , OPT  41 

7P 

K=B0IM(NCAT(JI,0PT4» 

71 

IF  (INAD(J).NE.O)  GO  TO  5 

72 

BLANK  OR  REJtCTEO  ITEM  f NCOUNTFPEO 

73 

DO  3 L=1,I<I 

74 

3 

ADM(J,L)=0.0 

75 

DO  4 L = 1 » K 

76 

4 

HOM(J,LI =0.0 

77 

GO  TO  9 

70 

5 

CON  T I NUt 

79 

c 

HO 

00  1 T = 1 » I NU P 

81 

15  (INAO(J).Nt.IT  M(I)|  GO  TO  0 

02 

DO  (■  L = 1,IKI 

M 

b 

A DM ( J , L 1 = A ( I , L I 

84 

i)0  7 L = 1 , K 

05 

7 

fiOMIJ.LI  =B  ( I ,L  I 

86 

GO  TC  9 

87 

5 

CONTINUE 

80 

INAIM  JI  = U 

09 

9 

CONTINUE 

90 

GO  TO  13 

91 

C 

9? 

£♦** 

ALL  I 1 F Mb  IN  POOL  have  BEE  N AOMIN  I STFPF  0 SO  POOL  OOr.S 

NOT 

93 

2*  ** 

NEtO  TO  rit  SF  ARCH!  0 

94 

10 

DO  12  1=1, M 

95 

I K I - All  IM  < NCA  T ( I I , 0 3 T 4 1 

96 

K = nil  M(NCAT  (I)  ,OPI’.l 

97 

no  jj  = i , iki 

90 

11 

A0M( 1 , JJ»= A(I, JJ) 

99 

no  12  J=l.< 

10P 

BOM (1 , Jl =0 (I , J) 

101 

12 

CONTI NUF 

1C? 

c 

103 

C 

104 

c*** 

PIIND“  THE  ITk  m PAs4MEIhrS  E X I T MA  Tt.  S CORJt  S°  TNDI  NG  TO 

1 HE  TTEMS 

105 

c*  *• 

III  I Ml  TEST 

106 

13 

IF  (OPTl.NE.il  GO  TO  14 

107 

00  19  1=1, M 

100 

IF  (INADIII  .EO.OI  SO  TO  14 

109 

IK  l-Mi)I  9 (I.CAT  ( 1 1 ,0  n T4  I 

111 

K=BD1M( NCA T(I), DPI .1 

111 

PUNCH  Jb,  INATM1I .IKEJ(I), (AnM«I,JJ| ,JJ=1,IKII 

112 

PUNCH  37,  (n0M(I,JI,J=l,K) 

113 

14 

cont  tmjf 

114 

C 

115 

C 

PP1NI  0U1  THE  I T t MS  ANO  THE  IP  PAFAMEURS 

116 

15 

DO  1 h 111=1, M 

117 

IKI=AOIM (NCAT ( 111) ,0PT4l 

110 

Kr.iCIM  (NCAT  ( I II  1 .0PT41 

119 

PPINT  30,  INAOdlll  ,IREJ(IIII.  (APM(III,JJI,JJ=1,IKI1 

120 

15 

PRIM  43,  (aDM(HI,J|  , J=1,KI 

121 

PRINI  43 

122 

C*.*» 

123 

C*** 

C*** 

c**' 

C*<M 

c*< 

17 


PF4G  SUOJ'CT  FROM  TAPF.l,  C AL  CUL  4 T E THr  T A , LOOP  BACK  IP  TOO  UNTIL* 
f H t 1 >•  APE  NO  MO«t  SUBJECTS  TO  SCORE 


N=  0 

NtAr  ( i , it-  trmi  nam-:, id,  ciPtspm  ,i=i,mi 
io=umt  (ii 

IF  (10. LI. 01  (.0  TO  10 


124 

125 
125 
127 
12* 

129 

130 

131 

132 
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R I vi  1 S6,  10 

IP  • 1 0, £ U , Q ) fcO  TO  29 

CHEC*  F OK  VALID  KrSPONSFS  ANO  FOP  SPECIAL  Pl.  SDONSF  VECTORS 
CALL  CHKOSP  ( IKESP,  IOMI  T , I NA  0,  NC  AT  , M,  NAME ,1  1, TNAOS) 

CALL  CHKVtC  (IKESP, INAQS.NCAT.W, ICHK.N9?  ST, NANS) 

bO  TO  ( 19*  21?  . 21*  2-TT  , rCMK 

ALL  HEPS  OMITTED 

PRIM  45,  NAME, ID 

T=-6U,00 

P£RCM=-50.0C 

bO  TO  22 

ALL  RESPONSES  INCORRECI 
PPINI  46,  NAME, ID 
I = - lt>  . 9 

PFPONr=0 .0 
GO  TO  22 

ALL  ‘•►SPbNSES  CORRECT 
PPINI  47,  NAME, ID 
T = 1 2 . C 
PE  RON  1 = 1,0 

SFOFM=0.U 
NUMI I S = C 
TINF 0=0,0 
GO  T(1  28 
COM  1NUE 

pfronT=FLG4T (F8EST ) /FLOAT  ( NANS » 


DEI: RHINE  WHICH  MOO1 L To  USF 
IF  ( OPT  4-2  > 24,25,26 

CALL  LOGRAO  < IRESP, ADM,MOM,M,NCAT, INAOS, . C31, c3,SFORM, I^AIL ,T,NUMI 
1 TS, TINFO,SE» 

GO  TC  27 

CALL  NOGKAO  ( IKES". ADM, BOM, M.NCAT, INAOS. ,001,60  , SF JRM,  I^AR  , T.NUMI 
1TS,T1NF0,SEI 
GO  TO  27 

CALL  NOMLOG  ( IKE  SO, AOM,  1DM , M , NC A T , I N AOS , .0 51, 60, S FORM, IF  A IL  « T , NUM I 
IT  S , T INFO, SET 

OUT  PL'  I RESULTS  TO  TAPES 
IF  (IFAIL.EQ.OI  GO  TO  26 
CONVERGENCE.  NOT  OBTAINED 
PRIM  41,  NAME, 10 
N C = NC  ♦ 1 


WRITF  ( 6,421  NAME  ,10, PEPCNT,T,SFC'»M, NANS,  NilMITS,!  INFO, St' 

N = N ♦ 1 
GO  TO  17 

PPINI  44,  N,  NC 

STOP 

FORMAT  (214, II, 211, IX, II, F5. 2, 277,12) 

FORMAT  (1616) 

FORMAT  (6011 ) 

FORMAT  ( 6 A 10 ) 

FORMAT  ( 6 A 10 ) 

FORMAT  ( T 50, *LINPSSO*, 7, T50* ======*, 777/77, T2  0, ’LINEAR  POLYCMOTOMU 

IS  SCORING  WITH  TWO  PARAMETER  MOOEL S* • 7 777 , T4 3 *PSYCHOMF TR ICS  MFTHOO 
2S  PROGRAM*, 7, T40*OF PARTMENT  OF  PSYCHOLOGY *, 7, T40*UN  IVE RSITY  OF  MIN 
3NESOI  A*,  7,  T4  0*MPLS.  MINN.  55*»55 *,  7777  , T2 0,  * I NUP*  , T 27  , • = * , 1 4 ,7, T? 0 
4,*MMAX*, T27,*=*l4,/.T20.*IOMIT*.T27,*=*.  1 4, 7. T20* OP Tl* . T27*  = * 14. 7 , 


>1  2J*0PT2*,T27*  = *,  1<*,/,T  , 0* OP T4* , T 27*=* , I4.,/,T?0,*MAXCAT*,T27,*  = *,I 
64,/T2'',*0*,T27,*=*,Fi..l,/T20,*VAF  IABLE  FORMAT*,*  FOR  POOL  =*,8A10, 
7/T2P, ’VARIABLE  FORMAT  f- OR  OAT  A = * , 8 A 10  , / T 2 0 , 8 A1 0 , / T 20  , 8 A 1 0 , / , T 2 J , 8 
8 A 1 0 ) 

38  TOR  MA 1 ( IS, IX, 12, cX , 10FR.2) 

3?  FO~M/\r  llOFF.c) 

38  FORMAT  (*PAPITy  FF?OR  ON  T APE  * , 1 1 IX, 12 T 

30  FPRMAI  C/13X,*ITfcM  10  = * , 1 5 , '•  X , * PE  j- ;C  T ION  =*  , I 3, /I  LX  , * A « *,10F6. 

1 2 I 

format  ( 1 2 X,  *0 1 • , 1 QF N . 2 ) 

FORMAT  ( 1 ox,  *CUNVi  :,.GLNCr  NOT  iotainfo  FOR  SUBJECT  =*,<-410,*  ID  =*, 
’AO) 

2 FOKM.T  ( i»,cAiO , A 4,  F t,.A  ,2F  f.  14,  14,. IF  7. 2> 

FORMS)  (////) 

FORfrtl  IHT.’CAbf;  3 1 A 0 = *,!'.,*  C A S r S NOT  OOF'VcPGfcJ  =*,I5) 
fORR**T  (10X,*ALL  ITcMS  OMIlTrO  F t - SUBJECT  =*,2A10,*  IT  =*,AN> 
F"P<M  < lTX,*SU"Jt 01  =*  ,241",*'  T r =»,AM,»  MAO  AIL  ANSMr-’S  I'4*,*COP 
1 h r C I * ) 

FORMA)  ( 1QX,  *GIJHJt  OT  =*.2Ain,»  10  =*,AM,*  MAS  ALL  A„bHEr'<;  CTR*,*PT 

;ct*» 

I NO 

IN)  i.f  P FUNCTION  A01M  C NIJMCA  T , VT  '.) 

INI'.  R U':T4 

Of  I ‘Ml  ML'.  TF>l  MU  ' If  R O'  A R A - A'V  TtF  0 FOR  A GTVFN  ITfm 
A 0 I m=  i 

IF  f I’  =»T  U . - 1.3)  Al)  I M = F.'jmo  A T ♦ 1 
»•*  T"-N» 

0 ID 

IM-l-i"  FUNCTION  DIM  ( > U'MCA  I , Or-l  a) 

1 NT  -1.fr  (TATA 

U t I - ’ F I '•  r i T F r NL'“MfP  0*  ' ’«  '^rT.>:  FOR  A G T V L‘  I*  IlfM 

" 0 I "'JMCr  ) 

;f  ii  ) T4 o, a)  to i m= ml m ? at ♦ i 

F ) i:f  I 
- f • 0 

I'  D I,T!I,‘  C H ik  INF  ( I .NIIR  , M,  JOT  i. ) 

1*41  c Ou  1 4 

*;m-  ■<  _ K-  -rRO-s  !*.  TH  :\"U»  OiT  a 

i,  ..*  fr-o-  i:  fo  .-d,  a ssa<  prtnt-  o ai  o im  - r r. s a ■* 

*■  A L I j 

1-  R Du 

if  ii  »!Np.l=.  ; c G » 1 o T ' : 

“‘lM  , I N1  D 
iff  •=; 

if  I * . L . INUP)  SC  10  ? 

PR  I*  I 4,  P.IM'P 
1‘  RR=1 

IF  <:  .L-‘ .ORTr.ANU.  J PI  4.12.  T>  SO  IT  ’ 

OF  I>.  I S,  pPT  4 

n rf  = i 

IF  I I r p"  • r J.D  STOP 

+>  r.  I L»k  *» 


FORMA  I ( 1 ' X,  * I NPtjr  iRRJr:  NO.  OF  IT-  Ml  IN  I T -M  t>03L  = » , 18 , / 1 j X , • N 
1U.  MUST  t*  . Lr  . I.f  * » 

FORMA)  ( 1 P X , * I NRU  T lP’C'I  *’0.  H-MS  AOmj  NI  iT'.RFp  = * , I T , / 1 C X , * NO  . 
i»'l>r  nr  .1-..  NO.  UF  1 1 1 "3  IN  *T,  m oool  =*,TS) 

8 FORMAT  ( i?x,  *1*4  PUT  E RF  ')•••  | CO-'TTN  4 =*,I3,*  OOtS  NOT  L T~'RE  S*  , *POnO 

1 TO  ANY  OF  TF_  A V 4 1 L A ML  -’c  r PN'45f  MONELS*  • 

*.N'l 

GIJMDNTINF  NOI.pAP  ( I»*  S',  A,  “,  1 1 T t MS,  NO  A T , IN  AO  , E °S  , M AX  I T , SF  OR*  , I F A I 
1 L , I'"  I A,l,UMITS,TlNFr,Gf  T 
i»T‘fi,AL  F OF  V*  P,S"P  Vf,0 

OIMt’MtCr,  IRISPTI.T),  ATiUu.l'T,  ' I 1 Ij  U , 1 0 I , NCATCOH,  INAOCuJI 
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’ 


. 


CALL  PISt  TT  (FORVNO,  11-fc  A,  <?,NirrMS,NCAr  , TNAn,<;,GJfS5» 
call  -j E W 1 -> A *->  (F  OWVNO,  Sll-'VNO,  1“’«-SP,  A,  I,  NITF  MS,  SCAT  , INAO,  MAVI  T ,t  PS,N 
I U M I T S » GOt  S3* IHlTA, S PR  V » ’FAIL! 

SFOF’«  = -5Dy  V 

IF  (1FAIL.E0.1T  tiCJ  TO  1 

CALL  NOINFO  ( A, H,NI TEMS,NCAT , TNAP, 1*- T4,T  INrOI 

SE-=i.c/sci-’r(APS(snpv) » 

WFTU*JN 

c 


1 

t iNFc=-9H.qq 

SF  =-HH.C)tj 

Rf  TURN 

F NO 

FUNCTION  IVALUt 

()  PC  SF , VC A ! > 

CHrl'fS  Vt  SPOT  St 

FOP  SPECIAL  CASES 

z*** 

IVALHr  PfcHJVNS 

. . . 1 IF  IP.  SO  IS 

IF  ST  P*.  .,PONSl 

i.  IF  IP-  SP  IS 

WORST  PfSPONSE 

c 

IV  AH’t=  5 

T OTHFPHTSF 

IF  (IPESH.FIO.  1)  IVALUF=! 

IF  llkF.SP.eO.CNCAT*in  IVALME-? 
me  turn 

t NO 

SUShuUTINL  fcCAfl  (A, a, THETA,  T mINBO , T *IN 31 , 3, E T07G ,£ TOZ1 T 
C***  COMPOTES  VALUES  NF  CC  SC  A c Y FOR  THF  CALCULATION  OF  TH;  DEV  IV- 
C**»  ATIVFS  OF  THE  NORMAL  OGTV'T  G P A OF  0 MOJtL  F UP  THE  SPECIAL  CASF 
C***  HHtN  IMESP  IS  THF  3- ST  PESPONSF 
C 

TMI*|ho  = THlTA-H 
TMINU1=0 . r. 
r=A*lMlNU0 
P=  C OF  N ( V ) 

e T07t=EXPC-Y*Y/2.  l 
t TO7l=0. U 
RE  TUFN 
ENO 

StlOkUUTlN  A CAT  N (A,  H,  I Ht  T A,  TMIr.'Bd,!  ‘'j  N U , p.cTOZO  ,ETOZU 
c»**  computes  v alik  s NtCFSsA  y fo=>  ihl  calculation  of  th"  uep- 

C**»  IVATIVES  OF  THF  NORMAL  OGIV1-  GPAtiF_0  UC O' L cOR  THF  SPECIAL 
C**»  CASi  WHFf*  IfFjF  IS  THL  WORST  ISSF'ONST 
C 

TFUNBQ  = 0 .0 
( MIN81=TM* TA-n 
Y 1 - A*  T M I NF  1 
P=1  .C-COF  MV1  ) 
t T07G  = 0. 0 

ETO?l  = EXM-Vl*Vl/c.0> 

KE  f UWN 
ENO 

SUBROUTINE  PCATOT  ( A , MO , HI , T M*  T A , T Ml W 0 , T MI N T \ , o, E T 070 , c T02 1 T 
C***  COMPUTES  VAUKS  NFCeSAMY  FOR  THF.  CALCULATION  OF  ncRIVATIVCS  OF 
C***  THE  NORMAL  OGIVE  GPADFD  “ODFL  FOE  ALL  OTHER  CASES 
C 

TMIN«(I  = THF  7 A-RO 

f M INH1  = THF TA-01 

T = AMMINMT 

Y1=A*TMINP1 

P = COFN«Y»-COFMYH 

tTOZJ=EXP(-Y*Y/?.OT 

ETOZ1=EXP(-Y1*YI/2.0) 

RETURN 

ENO 

FUNCTION  F Ou’ V’’0  ( I ?F  SP , I NA  D,  NI TF  MS,  NOAT  , A , H , T HE  T A I 
0 1 Me  NS  I ON  IPESPllOOT,  INAOdGOI,  NCAT(IOC),  A(100,10>,  0(100,10) 
C»**  CALCIILATFP  IHt  FI'  ST  OFPIVATTVE  FOR  THE  NORMAL  OGTVF  GRAOFO  MOOe  L 


. ’ —w-  * 


f, 

7 

ft 

Q 

1C 
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1? 

IT 

14 

15 
lfe 
17 

3 

4 
r, 
f, 
7 

ft 

0 

I A 

II 

1 

n 

7 

L 

c 

F 

7 

ft 

q 

n 

n 

l? 

IT 


3 

4 
F 
f, 

T 
A 
q 
in 
11 
12 
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SUM=p.O 

on  8 i=i,mteks 

IF  (1N«P(I>  .tO.O)  r.O  TO  5 
k=  i*f  sp  ( n 
J=IV«LU£TK,NCA1 (I)> 

GO  f n ( 1 « 2 , 3 » , J 

c 

c*»»  iPt-.ii  TC  »=ST  9Eb30NSc 

1 CALL  -'CAI1  (6  ( I. II  , )(1|<)  , TH  T A,  THIN".:.  THIN'U.o.E  T0ZC.rT071» 

GO  T u 4 
C 

c»  ♦*  IPt-H  18  WJPST  9F  ;°ONSi. 

2 :;Ul  -CAIN  (A  (I.ll  ,9  C I,K-1)  . THF  TA,  TMIG33.T  *IN'I1,P,!:  10Z«J.ET02i*  . 

GO  10  4 

**»  ALL  Jl  HF  M it  SF 0 Mb-  7 

3 CAUL  -CAinr  (MI.ll  .IMI.KI  ,!l(T,K-:l  , I HE  T A , T M I N9 o , T M 1 NlU  . P. E T OZ  J . E T 

1 0 / H 

4 C 0 M T i \UC 

c 

IF  ( " . K C • ' • C I ° = 0.'j2Ul 
vim.-.- mm-  a » i,  1 1 » (- 1 ozo -h  to  1 1 1 /- 

5 CCiiTIMUr 

F ihi/'.  ) = SU  VbOF  T ( c.  0 *3. 142) 

- L T'l-  '! 

F N'T 

F IJV  I iO-(  VWNO  t Kt  b--,  !NAO,  MiTr  M->,,MCA  T,A  , 3,  THETA) 
l 1 1 r.  > n g r^KC.iSi,  i N.n  < to 1 ) « -i  cat  (too i * 4<inc»i.>. 

C » • » LAL""L4T-  ' CfCO-C  O- F>  I V A T I V‘l  rOc  TH'  NOcmAL  Or-IVt  GCAOEn  MO:l-  L 

bl'-Tr,  . n 

VC  IT  *>I  = 1 . I /sn-'  T (2.  U »d . I <E) 

I • ° - 1 = 1* Milt 

IF  ( IMA;M  I) .Cu. J)  .0  TO  ' 

< = :-■■  •>->  < i ) 

J= I M-LU-  t< ,\CA T ( 1 ) I 
go  r - ( i , 1 , 3 ) , j 

c*».  i?.  j,  >»e'4t  p-b-^MSi 

1 CALL  -CM:  ( A ( I , ’ ) . ‘M  I , ^)  , TH-  1 a,  T-IM  10  , TM1M--1  ,P,  t T07  j,  - Tn,M  ) 

)0  I 4 

Q»».  I J.  ; - XS  40’ST  ?F ' °nk|Sr 

2 "AL^  -CAIN  ( A « : ,1  ) , T ( I , v-1 ) , T'l- Tf  , TMIN'3',  TM1NP1.P,- T02  .FTC7H 

GO  T r * 

c*.»  £L!_  ,jM-h  -Tf0o0Nb-C 

S Lua  - C A T n T (A  ( I,  *.  I ,H(  t ,<)  , I ( • ,K-;)  , M- TA  , T Ifno,  TMis,  l.  o,f  r^V.T,f  T 

IC/I  I 

4 r;c  n i vic 

IF  (‘-.--0.  .”1  oiC.LCCl 

Sl|l.  = A ( U»VG0IHI*  (ETflTC-LTii'l » /-' 

b'!M;=-S'IM:*<'llMl 

||M  - ( A IT  « i)  **  31  * 00)  ° r • ( U M°r»FT  0/  Jl  - CT  11991*)  T 07  it  1 /o 

il'l-  1 1 '1  ♦ S 1 1 T1  ♦ C U M 2 
•»  CTN1I.MIL 

:=cl.'* 

-Mil.-. 

. -JM 

• it-  t;*.  '.ClMFO  (A.t,  IT-MC,NCAT,X  IA0, THETA, TIMFO) 

I ' X ■ I'N  A«:;Ctl  l,  -M100.il>.  f CAT  ( ; J o I , INAOTIOGT 
• ••  - If  INFC- MAT  1 0*.  Ffj-  r.Ffi'IO  NORMAL  OK  IV-  “OflfL  AT  r.TVKN 
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lf> 
IT 
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2C 
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2T 
2* 
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31 
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'♦ 
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18 

17 

18 

19 
2: 
21 
22 

23 

24 
2F 
2F 

27 

28 
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C»»*  VALU-  OF  TH-Tp  9 

C 5 

r INF 0=0.0  9 

K00TP1=1.  I/SUM  (2.0*3.192)  7 

c*»»  LOO'-*  TVKk  1TFNS  A 

no  5 1 = 1, MIFFS  1 

if  (lNAnm.t-c.ci  -,o  io  s io 

C 1 1 

C**»  1NIT1ALI2ATI0*  .--VALUES  CALCULATE  D FOP  FTPS!  CATEGORY  OF  I IF  M I 12 

Y = A (1 , 1) * ( THt  T A-U< l , 1 ) > II 

PO  = cnFNm  ii» 

pi=1.l  in 

t U173=EXP(-Y*Y/2.  ?)  IF 

troz;=o.o  it 

K A I GP  Y = 0 11 

C 19 

C*»*  LOO.'  OVER  ALL  THE  OAlrGARIES  OF  I r F M I 20 

1 K A I GF'  Y = K A T Gc  Y ♦ 1 21 

P = o r . p l 2 T 

if  (m.fo.c.u)  y=u.caoi  2* 

F QRVP  = A ( 1, 1) *FOOT  ’T * (lT  )2U-t  TQ71 ) 29 

1 lNFn=T  INF  0*  (FONV- *F  D f V 71 ) / p 29 

C 29 

t T 0 / 1 = f T 0 / 0 27 

P1  = PC  29 

C 2° 

IF  (*ATGPY-(f.  CAI  (II  Fill  F.3,  4 30 

C***  CURRENT  CATtGORY  I'NDEP  CONS  I Of  K9 T I ON  IS  NOT  ONE  OF  IhF  FXTRcMFS  31 

2 V = A ( T,  1 1 * ( THE  T A-fl(  I , KAT  GRY  ) ) 32 

P0=COFN(YI  33 

i TOZl'=EXP(-Y*Y/2.ii)  39- 

GO  m 1 • 39 

0***  LAST  CATEGORY  OF  AN  I1-"  IS  Hr  INF,  CONSIOcPcO  39 

3 P0=l.u  37 

£1070=0.0  31 

GO  TO  1 30 

C***  ALL  CATE  GORIE  S Fir  IT-M  I HAV1  OF  N XAmINEO  90 

9 CONTI NUF  91 

C 92 

9 CONTINUE  93 

P F T U P N 99 

L Nl)  95 

SUBP'.’UTtN-  Lor-PACI  ( IRF  s A,  3.MII  NS,  NCA  r,  M-in.tPS.N  rir  P.3FOPM,  IFAI  1 

1L,  THt  T A,  MJBUS.T  1NF0,S£1  2 

F X I •'  f NAL  FDPVLL  *S09\/Ll  3 

OINFNSION  IPFGPIU01,  A(1CU,1G),  0(100,10).  NOA 1(103),  INAO(IOO)  9 

CALL  HJ  St  C I (FDPVLL . I^LSP, A, R, Nil EMS.NCAT  , T N A 0, 5 , GU f S3 1 5 

CALL  NF  W ) p Ap  (FORVLL,SU‘ VLL, IPESP,A,M, Nirf 13, PCAT.INAD, NITER, EPS, N 5 

111*1 1 Ts,GUtSS,  THE T A, 3 OR V,  IF  ail  ) 7 

SF  JPM  = -SU»V  R 

IF  (IFAIL.tO.ll  GO  10  1 9 

CAUL  LLINFO  ( A , 9 , f, I T £ m<j  , NC  A t , INAO,  THU  T A , T l*(FO  I 10 

Sfc  = l.D/SrjPT(A«S(SO  = V)  ) 11 

if  rtjpw  12 

C 13 

1 1 I NF(i=“99  • 99  19 

SE  =-99.99  15 

RETURN  IF 

FNO  17 

SUBROUTINE  CALC°  ( I RE  SP«  INAD.NCA  I , A, 9,  NITtHS,  THt  T A,  P|  1 

COHON  0 2 

OIltNSION  IFtSP(lOP),  INAOdOOl,  NCAT(IJO),  9(100,10),  P(100,?l*  » 3 

1(100,10)  9 

C**»*  CALCULATES  UPPER  AnO  LOWER  ° FOP  EACH  ITEM  WITH  GIVEN  ANSWER  VFCT  9 

UO  9 I=1,NITF*S  9 

IF  ( INAO  (1) .£0.01  GO  TO  9 7 


J 


O CJ  C J c 
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j=ipfs«mi> 

IF  (J.EO.l)  GO  TO  1 

M(i,i)  = i.n/(i.CFtxj(-o*«u,ii*(THLTA-n(r,j-im) 
go  ro  ■? 

P(I,1)=0.0 

IF  (J.EQ.(NCAT(I)Fl)  I GO  TO  3 

P(I,<.’»=l.n/(l.O+fXP(-n*MI,l»MTNFTA-P<I,J»)n 
GO  TO  4 
° ( I • <* » = X.  Q 

CONI  1 NIIF 
PE  TIIk’.X 
r.  NO 

FUNCTION  FIVVLL  ( I SP,  I NA  U,  N 1 1 1 M I , NC A T , A , 3 , T ML  T A ) 

commcn  n 

0 I Nr.  NS  ION  ->(100,2),  AIIOO.IO),  P<i00,13l,  TNAn(lOO),  I Pf  SP ( 1 CO ) » ’ N 
1 C A I (.’JO) 

CALdHATt"  F 1 *•  S T J °ll/AIiyF  OF  LOG-LIKE  FUNCTION 

stM  = r.  .o 

CALL  C AL  C c (I*'FSP,INAO,nCAT  ,A, P,  MITE  NS, THETA, °) 

DO  ’ I = 1,MTE«S 

IF  (IGAO(I)  .EP.O)  GO  TO  1 

S.U,1  = SUM+  A (I,  1)  * ( 1.  O-FM  1 , 1)  -P  ( 1,2  )) 

CONTI SUF 

f i>  vLL  = su”*n 
{ • 1 1 . «. 

• NT 

FUNGI  ION  LJPVLL  I L SP,  INAIO,  MT'*1i,NCAT  , A ,*,  THETA) 

CO  'T"lN  D 

Oil  -CION  3(110,2),  A(IJ0,10>,  •’(100,101,  I N A n ( i o 0 ( , I p-' 3P  d : 0 > , N 
i c a i ( : : o ) 

C AL  C'.'l  A Tf  G jFCONU  D^PIVATIVF  f LOGLIKE  FUNCTION 
diM  = , . 3 

CALL  LAI  OF  ( I A>t  SF> , I A 0 , ''C  A f , A , M , I TE  MS  , I HE  I A , P) 

03  I I=1,NITEMS 

IF  (1‘ 40(1). E 0.3)  GO  TO  1 

j : = : . c - ° ( i , : ) 

0 2-1.  . -p  ( ! ,2 ) 

p 1 = ( l , ? ) - p ( 1 , 1 ) 

i’If  f ;-2.  c»o:-:  .0 
0TF'r2  = 2.  U*li.-1  . 3 
p i -p( i , l ) * u 
<•'  ?’  ' ( 1 , 2 ) * J? 

sum:,  dm,  ( a (i , : ) •*.  / on*  (-nirF:*-'i,r)iFF2*Pc  - ( c-i-f  e ) **  2 ) / °i) 

LC  ,T|  Nil1 
;OPVLL  = IT*v*Sl'M 

PF  I'lLf. 

: N 3 

SIJ  (Of’t'T  I F PLAT  (u,  M , I H - TA  , TT  M,  N'  INC  AT  , aL  0«T;  < , P'/P^E  .* ) 

C ON  Ml  N n 

HIM'  ,^IOf.  I ( 1 C U , 1 ..  ) , PL'K-Pd  T,  Pt)PPEP(lO) 

CALCl'LAIt’  FO  all  EGPOV  p CATEGORIES  of  A GIV  E ITEM 

PLOH*-P(l)sO.O 
00  I I = 1,N')MCAI 

PIJ  ppf  p(I)=1,C/(1,  Ff*P( -D*  A*  ( (M  TA-tUITM,I>  I )) 

PL  Jw>  f ( I ♦ 1 ) sPUPPF ‘ ( I ) 

CONTI  fiUF 

PIJ  op,  ( VtJ**CAT  »1)  = 1 . c 

Ft  I II  G 

E MU 

SUOM'MT  IN-  LL1NFU  (A, IT,  T I T E MG  , NC  A T , I NA  3,  T Hc  TA.TIMFO) 

C 0 M“0.\  0 

ITIM-NSION  A(100,ic),  •’(100,10),  NCAT(IOO),  IN»3(100I 
IJ  I M-  ns  I ON  “(I  PT  E P ( 1 C ) , OLQWFP(IO) 

COM'UIES  XnFOFMATION  F 0 3 TH‘  GPAnTO  LOGISTIC  MQ0EL  AT  A GIVEN 
valu*  OF  HETA 


o o n n k*  o r\j  ►-» 


r inf  )=o. o 

DO  . I = iiMH"S 

if  ii'innm.-o.ji  jo  10  2 

call  pcai  <a  ( r , n , >,  i h»- t a,  i , \m  ( i) , -vjwL  i-,  oiipp-.f  i 

C*»*  LO'J-*  IVM-’  AIL  IH1  SPONSc  C4J&GC-IES  OF  IT  M I 

ncati  [,=nca  r m n 

DO  ; J=  1 « * CA  T f G 
DUD''.  y = l-c jppLp(j| 

QLGW-  k=l-PLCWi-  R CJ) 

R=R'l-*Rt  ■>  t Jt  -PI  OWL-  ( J ) 

IF  (P.FI.r.tl  p = U . 3 u J 1 

f {!VVP  = 0*A  < I,  1)  * (H'iupE  - ( J)  *QUD-  - - PLOW-  R ( J ) » JLOW>k  I 
I INK'  = TINT  Xf  URV  D • W o ) / p 

co  n r l Nti  f 

COnT  IVJt 
R b.  1 1 N 
F NO 

CUD-  'OH  N1*  NC<L  Cl.  ( IRr  o ’,A  , 1,  (Hi  NS,  xCAT,  I JA,),t  ’S,  Ml*  I T,  .■iF0RM,I':  A! 
1L  • I ml  1 A , NUNTTL,  T INF  Ci  IF  1 
t X |r  >.  MAL  tQ RV'  L 

uiK-'  jIon  i^lsp(i..d)  , a(i:o,i;i,  Hii’o.ioi,  n''4I(ioii,  iNAn(ini) 
call  'TNt  :T  < FDRV  NL  , I - c.  °,  a.  '<  • NITr  NS,  NO  AT  , t NAP,  G.r.J.  SH  r , &!lf  SL  0) 
CALL  St  C A F.  1 (FOWL,  I t SP, A , n. V I T-  NS , F CA T , l N AO, M A X I T , r , GU t SHI , GU 
It  SLO,, '.UNITS.  THlT  A,S0RV,  IF  AIL  ) 

«jF0i-,'<  = -Sl)PV 

IF  ( IF  AIL  .tn.l I GO  10  1 

CAUL  .LlNFl  ( A,  1,M  Tt  NS,  NCAT  , f NAD,  TH'  TA,  T INF  J1 
Si  = : . O/SUPT I A *S  ( S )°  V 1 ) 

V £ 1 IK  N 


f INK’r-NS. 

St  =-  Al.qq 
RF  r It  r. 
r Ml) 

TUNC  I ION  FOkVNL  1 I *' 5°.  I NAP, NtTc  MJ, NCA T , A , T, THE TA 1 
UIN-F.sIOl.  IHFSPU.n,  NCAT(HQ),  4(100,1.11,  0(100,13) 

»**  CALCULATES  FI1  ST  IlFPIvAlIVf  OF  MINIMAL  LOG!  , T TC  KINOIIIN 
***  NOT,  ! F ON  THIj  HintL,  TH  A”  S'  Af  IMF  SLOP-  PARAMETERS 
AN'-  TH  -r*S  ARt  1 Hr  IN(F  ►>(:«- PT  PA  AMFIM’S 


SUNsi'.O 

on  i = i,mtims 

if  ( 1 1 4T( t i . t*c  . 3 1 sc  in  J 

XNUMsy.a 

of  Nr=a . o 

NUMi.l’  I = NC  A T ( 1 1 ♦ 1 

on  . j = i,numc/.T 

CONS  I = A ( 1 , I R t S P ( I M -A ( I , Jl 

AP'»=A(I,JI*TMF  1 A ( I , j| 

t Z = t *P ( AkGl 

X NIJ  A-  XNU  M»  CONS  T »t  t 

OrN  )*:  Oh  F ON*  F 7 

CONTI  NUf. 

S(/N  = StlN*  XNUN/DE  \|Qh 
CONT I NUr 
F riHVNL  = Sll* 

RF  PIRN 
F Ml 

Sim-MUlTINr  n L INFO  I A , h , SIT  1 1FS,  NC  AT  , I *TA  (1,  T Hr  T A , T I NpO ) 

LI  I N‘  NS  I ON  A(1CU,1G>,  'M103,131,  NCAHIQJI,  INAOdOO) 

D Inf  US  ION  EdOl 

COMPIIJES  1 N F C F NA1  1 n F to.  THF  NOMINAL  LOGISTIC  FUNCTION  AT  A 
GIV..N  VALUt  OF  T M !.  T A 
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C*** 


c 

c*** 

c*** 

c**» 


c 

G**» 


C 

c*** 


2 


i 

4 


* w 

,*  * * 
' * * * 


Q • * » 

1 


2 

C 

J 


- * • • 
C»” 
c»»* 


I INFO=0. 0 

LOOP  OVFk  THE  ITEMS 
110  <.  IMtMTFMS 
IF  (INAOTII .EO.OI  50  10  4 
NO  A IK.sNCST  ( II  *1 

COMPUTE  IMF  DENCMImATOK  UStO  10  CALCULATE  p ANO  THfc  FI®ST 

OEKlVttTIWr  OF  P--IT  IS  1 HF  SAMfc  VALUE  FOR  ALL  POSSIBLE  ®F- 

S°ONSES  10  ITEM  I 

OE  NOm=Q.  o 

00  1 <=1,NCATLG 

L < IO  = ■-  X o < A < r , F I * T W‘r  A ♦ M(  I , K I I 

OE  NOf*=nr.NOM*t  (K) 

CONTINUE 

SUM  OVER  ALL  r0SSI9LE  Rf S“ONS;  S 

tjo  < j=i  .ncateg 

R>  = - ( J» /OEAOM 

CALICLATt  THE  F IF"  f TVRIVATIV'  OF  P 

n-  vuuf's  o.o 
no  2 R=l,\CATcG 

) EVMMSDHVVUM  + E 00»  { A (I  . Jl  - A ( I ,K  I I 

rosi : mie 

t ORV-’  = CE  < J » » r*1-  VNI|M>  / ((ic-  0M*0-'NOM) 

T I NRO=T  TI'ROR  (Fr»ov''»F  Dc  V'1!  rv 

CONTINUE 

c.n  i r i \ur 

RE  TORfl 
- M) 

SI/n-T'  TIN.  CHFRSP  < IR;  S'.  IUm:' , INAO.NC  AT  , NITF  MS.NAM  .lfl.INAOSI 
JIM-  tsiOF  I c E S °( 1 INAOCUJI,  OCAH100).  N4MEI21.  INATSTioOl 

U:  L I : S OMITTED  1 1 - MS  F0C  A vnjfnr  FrCJM  LI'T  OF  rr  **S 
nn.llMSIl  Vf.’  (R*-  S"L  T S R-  TtIRNLi  IN  TNAilJ*  ANO  CMFr<5  10  SEE 
IF  -LL  "tSPONOi-S  ..  >fc  VALID.  Ic  SOT  AN  r -*R09  M-SSAO;  IC  WRITTEN 
ON -I  l'*c  PROGRAM  MA'.TS 

')o  ' 1 = 1*  SITE  “5 

r.MMi  a l rNCA  I (11*1 
INAr  . ( I)  =]  NA[)C  I) 

IF  ( . E‘  R<  1)  .N-  . ) >*1  I 1 '■•0  TO  i 

ana  " ( i i =: 

VJ  »*•  2 

CP-  • «<  FOR  INVRLllI  RF  SPUNST 

IF  I • i.Lt.  Ifcrr®l  III  . AND.  ClRLr.PTI  » .LF  .E.JMCATTJ  GO  TO  2 
PR  I’.l  <*  NAMt  , 10,  I -'I  S M T I , N IMf-AT 

ST  r- 

CCN'  I'jtJR 
R-  T N‘  N 


FORMAT  ( K RRtv  IN  ) ' T A I 'I  J.GT  = »,£Ai3.*  10  = *,A9.«  HAS  TL 
ILU.Al  RCSFJNSt  = ♦.  U,/n».  °ON»£S  P'JST  ur  T *<Er.  N 1 ANO  NUMO 

(>■*  OR  OAI-GCFUS  =*,I.T,*  INCL'iOIV-  *1 
L N 0 

S'lMRf  J T I R.f  CHRVt  C I I Rf  S'1 , I N A n . NC'a  T , N I T rl  MS  , I CM  <*  NPE  S T , N A NS  I 

DIM  .SION  I-’tSHIl  j>*  INADIlo'l,  ICSTC1J0I 

CM<V  C IE  >TS  IMF.  RAW  FfSPO Nrr  VCTOP  TO  DtT  1MI NT  IF 

IF  II  IS  ON-;  LF  TM-*  F 'JL LOWING  tpo  L CASES! 

-IL  OE  >T  s S»l.  .'i  S.  . . ISM*  SrT  TO  < 

AIL  WO- ST  *•  E SPP'ISt'i  • • • TCMK  >fT  TO  2 

• ' 9-  SR  ON*  • ICR*?  St  T TO  * 

1 1 MR  vw  | SL I OPK  s-  T TO  m 


7 

8 
9 

10 

11 

12 

13 

14 

15 
IF 
1* 
18 

19 

20 
21 
2? 

23 

24 

25 
2F 
?7 
2« 
29 
32 

31 

32 
3* 

34 

35 
35 
37 

1 

? 

* 

. 4 
5 
F 
7 
0 
0 
1" 
li 
1? 
1 3 

14 

15 
15 

17 

18 
19 
21 
21 
2? 
?S 

24 

25 
2» 
27 

1 

2 


* * * 


F. 

7 

8 


o o o o 


NANS=0 
NW0RS  1 = 3 
NBtSI =0 

DO  1 I=1,NITEMS 
IF  UNAO(I).EO.O)  30  10  1 
N ANS=NANSf1 

IF  (IRESPIII ,L0. 11  NHr  ST=NRFSTtl 

if  (iKESPin .fo.ncat m fii  nhops r=NMOPsr* i 

1 CONTINUE 

c 

ICHK=- 

IF  (NANS.t'l.O)  I CM  < = 1 
IF  (NWOPST.t  Q.NANjI  ICMR^ 
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